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ABSTRACT 

The  more  important  results  and  developments  achieved  as  a direct 
result  of  this  contract  are  as  follows: 

A model  for  the  calculation  of  turbulent  flows  was  constructed  which 
incorporates  effects  of  Reynolds  stress  relaxation  and  the  generation  or 
destruction  of  turbulent  energy  by  rotation.  The  model  was  verified 
by  comparison  with  experiment.  It  was  then  used  to  predict  the  properties 
of  isolated  turbulent  vortices  and  applied  to  the  turbulent  mixing  layer, 
in  particular  to  the  formation  and  interaction  of  the  organized  structure 
in  the  mixing  layer.  An  alternative  approach  to  the  mixing  layer  in  which 
it  is  modelled  by  a rolling -up  vortex  sheet  between  fluids  of  different 
density  has  also  been  formulated. 

Calculations  of  strained  flame  elements  for  use  in  the  coherent 
flame  model  of  Ha  , F combustion  showed  that  only  at  very  high  strain 
rates  was  sufficient  vibrational  nonequilibrium  produced  to  be  of  interest 
for  a chemical  laser.  This  result  is  reflected  in  the  fact  that  only  in  the 
early  portion  of  the  jet,  where  strain  rates  are  a maximum,  does  this 
vibrational  nonequilibrium  occur. 

Experimental  results  were  obtained  in  a new  facility  and  by  a 
new  technique,  for  the  amount  of  reaction  product  obtained  between  two 
chemically  reacting  aqueous  streams.  Results  at  high  and  low  Reynolds 
number  showed  the  dependence  of  the  mixing  on  Reynolds  number  and 
Schmidt  number.  A transition  Reynolds  number  region  was  found  and  the 
nature  of  the  generation  of  small  scale  motions  which  dominate  the  forma- 
tion of  reaction  product  was  explored. 

A new  facility  to  explore  energetic  reactions  with  substantial  heat 


release  has  been  designed,  is  being  built,  and  should  be  used  for  the 


first  measurements  in  the  summer  of  1979.  The  design  problems  were 
formidable  but  the  facility  is  unique  and  offers  the  prospect  of  obtaining 
new  and  very  important  data  with  which  to  compare  existing  models.  As 


a result  of  the  design,  the  range  of  problems  and  parameters  that  can  be 
studied  in  the  facility  and  the  possibilities  for  the  development  and  applica- 
tion of  new  instrumentation  and  diagnostics  are  very  great. 

State-of-the-art,  very  high  speed  and  precise  laser  Doppler  velocity 
measurements  were  completed  under  this  contract.  The  acquired  expertise 
will  find  application  in  the  turbulent  combustion  experiments.  It  has 
already  been  applied  with  considerable  success  up  to  a Mach  number  of 
2.2.  The  extension  of  two -point  velocity  measurements  has  now  been 
made  and  a multipoint  multichannel  LDV  system  has  been  designed  and 
is  nearing  completion. 

A laser  induced  fluorescence  technique  was  developed  and  success- 
fully applied  to  the  observation  of  turbulent  mixing  in  a water  jet.  It 
allowed  the  direct  observation  of  a concentration  field  to  the  smallest 
turbulent  scales.  Both  direct  photographic  techniques  and  imaging  of  a 
line  on  a 1024  detector  Reticon  array  have  been  used  to  obtain  the  data. 
Image  processing  will  give  quantitative  measurements  of  the  concentration 
field  and  potentially  an  inferred  velocity  field.  / 
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n.  INTRODUCTION 

The  general  aim  of  the  research  carried  out  under  this  contract 
was  a combined  attack,  both  experimental  and  theoretical,  on  the 
problem  of  trying  to  predict  and  understand  chemically  reacting 
turbulent  flows.  This  major  contract  made  it  possible  to  focus  a wide 
diversity  of  interest  and  expertise  at  Caltech  on  the  problem.  As  a 
result,  four  separate  but  overlapping  tasks  were  identified  and  described 
in  the  first  proposal,  two  concerned  with  predictive  models  and  two 
with  experiments. 

Section  III  is  a summary  of  the  progress  made  in  each  of  these 
tasks.  Section  IV  lists  the  publications  derived  from  this  research 
and  also  lists  the  seminars,  presented  papers  and  invited  lectures  given 
by  GALCIT  faculty,  research  fellows  and  students  on  the  work.  Section 
V is  an  appendix  which  includes  details  of  the  research  referred  to  in 
each  Task  summary  but  may  not  yet  be  available  as  a final  publication 
in  the  open  literature. 

This  research  was  sponsored  by  the  Air  Force  Office  of  Scientific 
Research  (AFSC)  and  the  Air  Force  Weapons  Laboratory  (AFSC),  United 
States  Air  Force. 
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III.  GENERAL  SUMMARY 

TASK  I.  A Study  of  Turbulence  Models  and  Applications 
Principal  Investigator 
P.  G.  Saffman 

The  problem  of  developing  a turbulent  model  that  would  allow  the 
quantitative  examination  of  the  interaction  or  organized  turbulent  structures 
was  considered.  One  of  the  main  features  of  the  large  eddies  is  the  presence 
of  mean  vorticity  or  mean  solid  body  rotation.  To  deal  with  their  interaction, 
it  is  necessary  to  have  a turbulence  model  which  has  been  demonstrated  to 
be  effective  in  the  presence  of  rotation.  Two  important  physical  effects 
need  to  be  taken  into  account.  First,  there  is  the  generation  or  destruction 
of  turbulent  energy  by  rotation  and,  second,  there  is  the  relaxation  effect 
of  the  Reynolds  stresses  in  following  the  rotation  of  the  principal  axes  of 
strain. 

A model  was  developed  which  appeared  to  contain  the  essential  physics 
of  turbulence  in  rotating  fluids  and  tested  by  comparison  with  the  experi- 
mental data  on  turbulent  flow  in  the  gap  between  rotating  cylinders.  The 
model  was  then  applied  to  determine  the  structure  of  an  isolated  turbulent 
vortex  and  the  coalescence  of  turbulent  large  scale  structures  in  the 
turbulent  mixing  layer.  In  the  first  application,  the  existence  of  a 
similarity  structure  for  the  decaying  turbulent  vortex  was  demonstrated 
with  an  inner  laminar  core  whose  size  is  inversely  proportional  to  the 
Reynolds  number.  In  the  second  application,  a turbulent  time -dependent 
spatially -homogeneous  mixing  layer  was  calculated  from  the  model 
equations.  A small  disturbance  with  wavelength  equal  to  that  of  the  most 
unstable  disturbance  of  the  associated  Orr-Sommerfeld  problem  was 
imposed  and  followed  to  finite  amplitude.  This  temporally  developing 
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mixing  layer,  instead  of  the  spatially  developing  one,  was  considered  for 
reasons  of  computational  cost.  The  growth  of  the  disturbance  into  two 
noticeable  eddies  which  then  coalesced  was  followed.  It  was  concluded 
that  turbulence  models  can  be  constructed  which  will  allow  the  description 
of  the  interaction  and  coalescence  of  organized  structures. 

Work  was  also  started  on  an  alternative  formulation  of  the  mixing 
layer  in  which  it  is  modelled  by  a rolling-up  vortex  sheet.  The  equations 
of  motion  of  the  vortex  sheet  have  been  formulated  as  a singular,  non- 
linear integro -differential  equation.  The  mathematical  problems 
association  with  this  equation  are,  however,  severe. 
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TASK  II.  The  Coherent  Flame  Sheet  Model  for  Turbulent  Chemical 
Reactions 

Principal  Investigator 
F.  E.  Marble 

In  order  to  adapt  the  coherent  flame  model  to  the  H-F  chemical 
laser,  an  extension  of  the  theory  was  demanded  that  developed  into  a 
more  serious  task  than  had  been  anticipated.  Before  beginning  the 
analysis  of  turbulent  jet  structure,  calculations  of  laminar  diffusion 
flames  under  various  rates  of  strain  were  required.  This  was  accom- 
plished through  use  of  the  Blottner  code,  obtained  from  TRW  and  the 
Defense  and  Space  Systems  Group,  and  reworked  for  the  computing 
facilities  at  Caltech.  In  order  to  achieve  the  degree  of  vibrational  inver- 
sion that  we  considered  appropriate,  it  was  necessary  to  strain  the  flame 
at  rates  far  in  excess  of  other  computations  we  have  made.  Strain  rates 
of  104  sec-1  showed  a first  indication  of  non -equilibrium  vibrational 
states  while  strain  rates  of  107  sec-1  achieved  a large  degree  of  total 
inversion.  This  strain  rate  was  accompanied  by  a distinct  change  in  the 
laminar  flame  structure. 

Where  the  strain  rate  was  generally  less  than  10®  sec-1  , the  flame 
behaved  very  accurately  as  if  the  reaction  rate  was  infinitely  fast,  at 
le<  st  so  far  as  the  reactant  consumption  rates  are  concerned.  Thus  the 
reactant  consumption  in  this  range  was  ~ V D e where  D is  the  appropriate 
binary  diffusion  coefficient  and  e is  the  strain  rate.  The  numerical 
calculation  gave  the  factor  of  proportionality  and  its  dependence  upon 
reactant  concentration.  The  strain  rate,  entering  only  in  the  manner 
shown,  demonstrated  that  the  consumption  rate  was  controlled  by 
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diffusion  and  totally  insensitive  to  chemical  rates.  On  the  other  hand, 
when  the  strain  rate  was  greater  than  106  sec’1,  the  reactant  consump- 
tion varied  as  e Examination  of  the  results  showed  that  the  flame 

was  chemically  limited  in  this  range  and  that  the  chemically  active  zone 
filled  the  entire  flame  thickness.  It  was  then  clear  that  the  proper  form 

of  the  reactant  consumption  rate  for  large  values  of  e was  VD/e 

Tc 

where  Jd/c  is  proportional  to  the  thickness  of  the  reaction  zone  and 
tc  is  a chemical  time.  For  the  entire  range  of  strain  rates,  the  reactant 


consumption  rate  is 


J De 

1 + t e 
c 


which  behaves  as  J De  for  strain  rates  e <<  1/ Tc  and  as  (1/tc)  J D/ e 
for  e » 1/  T . 

To  accomodate  this  more  complex  burning  law,  the  entire  integral 

formulation  was  reworked.  This  extended  model  was  used  in  the  study 

of  a circular  jet  of  hydrogen  issuing  into  a background  of  stationary 

atomic  fluorine.  Calculations  were  made  of  the  distributions  of  flame 

surface  area  chemical  species,  and  vibrational  states  of  HF  at  several 

values  of  z/d  along  the  jet  axis.  The  calculations  were  done  also  for 
W0  T c 

two  values  of  — g — , where  W0  is  the  initial  jet  velocity  emitted  from  a 
nozzle  of  diameter  d.  This  parameter  essentially  determines  the  quantity 
rce  and  hence,  according  to  the  above  relationship  for  the  reactant 
consumption  rate,  accentuates  the  type  of  flame  structure  that  leads  to 
non -equilibrium  of  the  vibrational  states  of  HF. 

This  work  is  covered  in  detail  in  Appendix  C which  will  be 
available  as  a GALCIT  Technical  Report. 
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TASK  III.  Dilute  and  Energetic  Chemical  Reactions  in  Turbulent  Flow 

Principal  Investigators: 

A.  Roshko  and  G.  L.  Brown  (since  October  1978) 


A.  Isothermal  Chemical  Reactions  in  Turbulent  Flow 

The  blow-down  water  channel  facility  was  utilized  fully  during  the 
present  reporting  period  for  several  investigations  of  turbulent  free  shear 
layers  and  wakes.  The  principal  experimental  techniques  made  use  of 
the  visible  (red)  reaction  product  which  is  produced  when  the  reactants 
phenolphthalein  and  sodium  hydroxide  mix  chemically  in  the  turbulent 
mixing  region.  Flow  photographs  and  movies  were  used  to  obtain 
information  about  the  large-  and  small-scale  structure  of  the  flows 
and  a light  absorption  technique  was  used  to  infer  the  instantaneous 
amount  of  visible  product. 

The  main  results  of  the  past  year's  research  are  as  follows: 

1)  At  high  Reynolds  number  (Re)  the  mixing  in  the  shear  layer  is 
independent  of  Re  and,  based  on  a comparison  with  previous  measurements 
of  Konrad  in  a gaseous  flow,  at  most  only  weakly  dependent  on  Schmidt 
number  (Sc  = ^ ). 

2)  At  low  Re  the  mixing  is  a strong  function  of  Sc. 

3)  A transition  region  exists  in  which  small  scale,  three-dimen- 
sional motions  are  introduced  and  develop  in  the  layer.  They  are  super- 
imposed on  the  large,  quasi -two -dimensional  vortex  structures.  These 
small  scale  motions  generate  a very  large  reaction  interface  area  and 
strain  rate,  thereby  permitting  the  mixing  rate  to  be  as  large  as  the 
entrainment  rate  of  fresh  reactant  into  the  shear  layer.  The  transition 
appears  to  be  determined  mainly  by  the  initial  conditions  of  the  shear 
layer  rather  than  by  the  large  scale  Reynolds  number.  The  parameters 
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which  seem  to  best  characterize  the  transition  are  the  development 
distance  x/Qj  based  on  the  initial  momentum  thickness  and  the  Reynolds 

u,  e, 

number  — ± of  the  boundary  layer  on  the  high  speed  side. 

4)  The  wakes  from  a blunt  and  from  a sharp  trailing  edge  of  the 
splitter  plate  were  also  founa  ,o  exhibit  a mixing  transition  across  which 
the  aqueous  mixing  increases  by  an  order  of  magnitude  or  more.  Again 
the  physical  explanation  is  in  the  emergence  of  small  scale  motions  in 
the  transition  region.  At  large  Re,  the  wake  mixing  is  likewise 
independent  of  Re. 

Reference 

Konrad,  J.  H.  "An  Experimental  Investigation  of  Mixing  in  Two-Dimen- 
sional Turbulent  Shear  Flows  with  Applications  to  Diffusion  Limited 
Chemical  Reactions ",  Ph.  D.  Thesis,  1977,  California  Institute  of 
Technology. 
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B.  Energetic  Chemical  Reaction  in  Turbulent  Flow 

Work  on  the  overall  design  and  fabrication  of  this  facility  continued. 
The  major  recent  work  has  been  in  the  final  design  of  the  upstream 
contraction,  test  section,  and  downstream  gas  cooling,  scrubbing  and 
collecting  region.  Construction  is  presently  underway  at  the  Caltech 
Central  Engineering  Shop.  The  major  concerns  have  been  material 
compatibility  with  Fluorine,  compatibility  with  corrosive  and  high- 
temperature  HF,  machining  costs  and  overall  safety.  Safety  remains 
a prime  concern  as  both  Fluorine  (Fa)  and  the  reactant  product 
Hydrogen  Fluroide  (HF)  are  highly  corrosive  and  extremely  toxic. 

The  main  results  of  the  past  year's  work  are  as  follows: 

1)  The  apparatus  has  been  moved  to  what  was  formerly  the 
"Clean  Air  Car  Laboratory"  from  its  former  location  in  the  Guggenheim 
penthouse.  The  major  reasons  for  this  move  are  (i)  availability  and 
(ii)  preferable,  particularly  from  the  safety  standpoint  as  this  location 
is  much  more  accessible,  affords  more  room  for  working  and  allows 
more  isolation  in  terms  of  the  exhaust  gas  containment.  Since  the 
experiment  is  to  be  performed  here  in  the  heart  of  Caltech,  containment 
is  considered  one  of  our  strongest  safety  requirements  because  of  the 
toxicity  of  the  Fa  and  HF  and  their  low  threshold  limits  (0.  1 ppm  for 
Fa;  3 ppm  for  HF).  The  pipe  link  that  joins  the  high  pressure  supply  to 
the  reactant  holding  tanks  has  also  been  completed. 

2)  The  two  adjustable  metering  valves  have  been  built.  Each  is 
made  of  a 316  stainless  steel  housing  into  which  the  reactant  gas  enters 
and  then  chokes  through  a series  of  8 rectangular  slots  which  are  milled 
into  a Monel  tube.  The  choked  area  is  variable  by  means  of  a movable 
stainless  steel  piston.  With  100  psi  (6.9  x 10B  N/ma)  inlet  pressure 
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schlieren,  L.DV)  of  the  mixing  layer  in  the  sideview.  The  top  and  bottom 
walls  are  also  of  glass  and  allow  for  similar  study  of  the  planform  of  the 
mixing  layer.  A traverse  mechanism  platform  can  be  substitute  for 
the  top  and  bottom  walls,  while  accessibility  remains  through  the  side- 
walls.  Some  parts  of  the  traverse  mechanism  will  be  subjected  to  hot 
HF  attack  and  those  parts  have  been  made  replaceable.  Provision  has 
also  been  made  for  film  cooling  of  the  sidewalls  by  injecting  cool 
Nitrogen  at  the  beginning  of  the  mixing  layer.  Etching  of  the  glass  by 
the  HF  and  thermal  shock  due  to  heat  release  are  considered  the  major 
problems  and  film  cooling  may  be  required  only  in  the  maximum  heat 
release  cases  (temperature  rise  of  1000°  C). 

Drawings  for  this  section  are  complete  and  construction  will  begin 
immediately  upon  completion  of  the  contraction  section. 

(c)  The  diffuser/afterburner  section:  This  divergent  section 
allows  the  reaction  to  go  to  completion  before  the  exhaust  gases  are 
discharged  downstream  into  the  exhaust  gas  section. 

(d)  The  exhaust  gas  section:  The  design  problems  associated 
with  collecting  or  processing  the  hot  and  highly  toxic  exhaust  gas  in  a 
built-up  environment  like  Caltech  are  formidable.  Because  of  the 
large  flow  rates  associated  with  this  facility,  processing  of  the  hot 
exhaust  gases  on  the  fly  and  discharge  to  the  atmosphere  were  initially 
considered  as  impractical  because  of  the  extremely  low  concentration 

of  exhausted  HF  or  Fa  which  could  be  considered.  In  particular,  packed 
absorber  beds  commonly  used  in  chemical  laser  or  rocket  engine  work 
are  totally  unsatisfactory  because  of  size,  pressure  drop  and  time  to 
reach  steady  state  conditions.  Some  typical  magnitudes  for  the  present 
facility  are  as  follows: 
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mass  flow  rate  on  one  side  of  layer 

maximum  molar  concentration  of  Fa 

maximum  temperature  rise 

initial  volume  (2  reactant  tanks)  at 
100  psi  (6.9  x 106  N/ma 

typical  heat  release 

final  volume  at  ~ 1000°  C and 
1 atmosphere  pressure 


~ 1 kg/sec 
~ 5i 

~ 1000° C 
~ 30  ft*  (.  85m*) 

~ 10*  J/s 
~ 800  ft*  (22.  6 m*) 


Collecting  800  ft*  of  corrosive  and  toxic  gas  at  1000°  C,  and  then 
subsequently  scrubbing  presents  a substantial  problem.  A heat  exchanger 
and  gas  containment  gas  were  designed  for  the  facility  and  appeared  to 
be  a possible  solution.  The  design  is  briefly  described  below.  More 
recently,  however,  after  discussions  with  J.  Ortwerth  and  using  some 
results  on  HF/F3  scrubbing  obtained  at  the  Boeing  Research  Laboratory, 
it  has  been  decided  to  replace  the  heat  exchanger  and  gas  containment 
bag  by  a system  in  which  5#  sodium  hydroxide  solution  is  sprayed  in  a 
fogjet  into  the  exhaust  gas  to  cool  it  to  about  100°  C and  partially 
neutralize  the  HF  and  Fa.  The  final  volume  to  be  handled  will  now 
be  significantly  reduced  and  probably  no  more  than  300  ft*  (8.  5 m® ). 

The  exhaust  gas  and  carryover  spray  are  then  caught  in  a large  teflon 
bag  and  subsequently  scrubbed  and  discharged  to  the  atmosphere,  this 
latter  process  occurring  at  a much  slower  rate.  Some  additional  work 
remains  to  be  done  on  the  caustic  spray  system  before  the  design  is 
complete. 

(1)  Heat  exchanger  and  gas  containment  bag:  This  design 
was  pursued  at  some  length  because  of  its  many  intrinsic  advantages. 

It  will  not  now  be  built,  however,  and  has  been  superseded  by  the 
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scrubbing  system  outlined  above.  The  exchanger  was  designed  to 
remove  approximately  10®  J/sec  for  four  seconds  up  Lug  a thermal  mass 
concept.  Fifteen  laminar  flow  elements,  each  consisting  of  many 
closely  spaced  stainless  steel  plates  precooled  with  liquid  nitrogen, 
would  reduce  the  exhaust  gas  temperature  to  typically  100°  C.  This 
exhaust  gas  is  then  trapped  in  a large  stainless  steel,  very  thin  sheet, 
bag  of  a metal  bellows  design  (10  ft  by  10  ft),  A teflon  bag  was  not  used 
so  that  the  system  would  be  failsafe.  As  before,  the  trapped  gas  would 
then  be  neutralized  at  a slower  rate. 

(2)  Gas  feed  system:  Finally,  the  supply  lines  and  valves 
for  transferring  Fa  from  the  high  pressure  supply  bottles  to  the 
reactant  tanks  remains  to  be  built.  We  have  had  useful  inputs  from 
JPL  on  this  design  and  forsee  no  major  problems  here. 
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TASK  IV.  Advanced  Instrumentation  and  Measurement  Techniques 
Development 

Principal  Investigators: 

P.  E.  Dimotakis  and  H.  W.  Liepmann 

A.  Laser  Doppler  Velocimetry 

1)  High  speed  flow  measurement.  The  experiments,  data 
reduction  and  analysis  on  turbulent  boundary  layer  measurements  of 
mean  u,  v,  u',  v',  and  u'v'  in  high  speed,  high  Reynolds  number  turbulent 
boundary  layers  were  completed  during  the  contract  period.  Preliminary 
results  were  presented  at  the  Third  International  Workshop  on  Laser 
Velocimetry,  at  Purdue  University  in  July  1978  (see  Appendix  D).  The 
final  results  will  be  published  as  an  AEDC  report,  under  whose  sponsor- 
ship the  work  began.  A pre-print  of  the  report,  by  Dimotakis,  Collins, 
and  Lang  is  appended  (Appendix  E).  We  believe  these  measurements 
constitute  the  state-of-the-art  in  precision  high  speed  laser  Doppler 
velocimetry,  which  also  specifically  address  the  difficulties  in  making 
measurements  in  the  vicinity  of  solid  boundaries  (walls). 

In  these  measurements,  a flow  Mach  number  range  of  0.  1 to  2.  2 
was  covered  in  two  differont  facilities  (GALCIT  Merrill  Wind  Tunnel 
at  M ~ 0.  1 and  the  20"  JPL  Supersonic  Wind  Tunnel  for  the  higher  speed 
flows).  As  a result  of  our  measurements,  we  were  able  to  show  con- 
clusively that  the  discrepancies  that  had  arisen  in  the  past  in  high 
speed  flow  measurements  of  the  Reynolds  stress  in  the  vicinity  of  walls 
were  unrelated  to  compressibility  effects,  as  opposed  to  what  had 
originally  been  claimed.  We  were  also  in  a position  to  make  investiga- 
tions into  the  nature  of  turbulent  transport  in  the  vicinity  of  the  viscous 
sublayer  which,  even  though  it  is  usually  only  a small  fraction  of  the 
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turbulent  boundary  layer  thickness  (typically  ~ 0.001  boundary  layer 
thicknesses  at  high  Reynolds  number),  can  have  significant  effects  that 
extend  over  the  lower  20$  of  the  turbulent  boundary  layer. 

2)  Multipoint,  multichannel  laser  Doppler  velocimetry.  The 
design  and  fabrication  of  the  electronic  components  of  the  multipoint, 
multichannel  laser  Doppler  velocimetry  system  is  proceeding  satis- 
factorily even  if  a little  behind  the  schedule  we  had  hoped  would  have 
been  possible.  The  design  is  by  now  complete  and  the  effort  is  presently 
expended  in  the  fabrication  and  assembly  in  a manner  that  will  render 
the  manufacture  of  additional  channels  as  inexpensive  as  possible. 

To  this  end,  computer  aided  design  and  computer  aided  manufac- 
turing (CAD/ CAM)  techniques  have  been  used  wherever  practical.  In 
particular  the  special  purpose  wire  wrap/printed  circuit  boards,  on 
which  the  system  is  being  assembled,  were  designed  on  the  computer, 
which  also  computed  and  stored  the  coordinates  for  the  7,  000  holes 
that  had  to  be  drilled  and  generated  the  artwork  on  a digital  plotter. 

The  resulting  plots  were  then  used  to  generate  the  negatives  used 
to  etch  the  printed  circuit  boards.  Reduced  copies  of  the  resulting 
board  etch  pattern  are  included  in  Appendix  F.  The  drilling  was 
subsequently  done  by  using  the  computer  to  feed  the  stored  coor- 
dinates of  the  holes  to  a digital  3 -axis  positioning  table  on  which  a 
precision  45,  000  RPM  drill  motor  (fitted  with  tungsten  carbide  drills) 
was  mounted. 

As  of  this  writing,  the  high  speed  analog  front  end  signal 
processing  circuitry  for  four  channels  has  been  fabricated.  We  anti- 
cipate that  the  first  measurements  using  the  new  system  will  be  realized 
during  the  current  calendar  year  and  that  it  will  be  available  in  time 
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to  be  used  in  the  high  heat  release  combustion  facility  presently  under 
construction. 

Concurrently  with  the  design  and  fabrication  of  the  multichannel 
system,  we  were  able  to  measure  at  the  same  time  the  velocity  at  two 
points  in  a high  Reynolds  number,  two-dimensional  shear  layer.  This 
was  possible  through  a modification  of  the  existing  processor  that 
allowed  the  signals  from  two  optical  channels  to  be  time  multiplexed  into 
a single  digital  processor,  thereby  essentially  measuring  the  selected 
velocity  component  at  two  points  in  the  flow  simultaneously.  The  results 
are  summarized  in  a paper  submitted  to  the  ALAA  Journal  for  publication. 
A pre-print  of  this  paper  by  Koochesfahani,  Catherasoo,  Dimotakis, 
Gharib  and  Lang  is  appended.  The  measurements  and  resulting  new 
conclusions,  that  were  made  possible  by  our  ability  to  monitor  velocity 
at  one  additional  point  in  the  flow,  suggest  to  us  that  our  efforts  at  truly 
multipoint  measurements  should  yield  a wealth  of  new  information  about 
turbulence. 


, - 
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B.  Laser  Induced  Fluorescence. 

One  of  the  most  rewarding  parts  of  the  research  under  this  task 
has  been  the  development  of  the  Laser  Induced  Fluorescence  (LIF) 

■ technique.  This  has  allowed  the  positive  identification  of  the  concentration 

field  of  one  of  two  aqueous,  mixing  species  that  is  labeled  in  terms  of  a 
dye  that  can  be  selectively  excited  to  fluoresce  by  a laser  illuminated 
field  (beam,  sheet  or  volume).  The  technique  not  only  allows  flow 
visualization  under  conditions  that  are  inaccessible  by  conventional  flow 
visualization  techniques,  but  also  allows  quantitative  measurements  to 
be  made  since  the  fluorescence  intensity  emitted  from  a point  in  the 
flow  is  linearly  dependent  on  the  concentration  of  the  fluorescing  dye 
and  the  local  intensity  of  the  laser  field. 

As  of  this  writing,  two  methods  of  acquiring  data  that  result  from 
this  technique  have  been  used.  The  first  method  utilizes  the  laser  as  a 
sheet  illumination  source  and  records  the  data  photographically  by 
focusing  the  plane  illuminated  by  the  laser  on  the  image  plane  of  a 35  mm 
camera.  Examples  of  data  recorded  in  this  fashion  appear  on  Figures 
2a,  6,  7,  8,  and  9 of  Appendix  H.  The  second  method  utilizes  the  laser 
as  a line  illumination  source  (beam)  and  images  the  illuminated  line  on 
a linear  charged  coupled  device  (CCD)  detector  array  of  128  or  1024 
detectors  that  are  electronically  scanned  at  high  speed  to  produce  a 
sequence  of  synchronous  voltages  that  are  proportional  to  the  light 
intensity  incident  on  the  corresponding  detector.  These  voltages  are  then 
used  to  modulate  the  intensity  of  a cathode  ray  tube  beam  which  is  driven 

! «• 

by  a raster  scan  pattern  synchronously  with  the  electronic  scanning  of  the 
, detector  array.  The  resulting  pattern  on  the  CRT  screen  is  then  photo- 

graphed to  produce  an  image  corresponding  to  the  concentration  of  the 
labeled  mixing  species  on  the  illuminated  line  as  a function  of  time. 


* 
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Examples  of  data  recorded  in  this  manner  appear  on  Figures  4 and  5 
of  Appendix  H. 

Several  important  discoveries  were  made  possible  by  this  new 
flow  diagnostic  method.  It  was  possible  to  determine,  for  example,  that 
entrainment  of  irrotational  fluid  into  the  turbulent  region  of  a jet  was 
by  means  of  a mechanism  not  significantly  different  from  the  mechanism 
that  was  found  in  the  two-dimensional  shear  layer.  In  particular,  it  was 
found  that  entrainment  does  not  proceed  as  the  result  of  a diffusive 
process  of  a turbulent  - non-turbulent  interface,  but  rather  as  a result 
of  the  induced  velocity  as  a result  of  the  large  structures  that  are  also 
found  in  the  turbulent  jet.  Secondly,  it  was  established  that,  at  least  in 
the  case  of  water  (high  Schmidt  number^it  was  possible  to  find  unmixed 
reservoir  fluid  all  the  way  to  the  axis  of  the  jet,  contrary  to  the  conven- 
tionally accepted  ideas  in  this  matter.  Thirdly,  it  was  also  established 
that  the  notions  of  intermittency  and  a turbulent  region  that  is  separated 
from  the  non-turbulent  region  by  an  interface  that  is  topologically  relatively 
simple  are  of  questionable  validity.  Lastly,  evidence  was  found  that  the 
structures  in  the  turbulent  jet  grow  in  size  by  a mechanism  akin  in  many 
ways  to  the  coalescence  process  that  is  dominant  in  the  two-dimensional 
shear  layer.  These  are  important  observations  because  they  suggest  that 
the  organization  of  turbulence  that  was  found  to  exist  in  the  two-dimensional 
shear  layer  may  also  be  found  in  other  turbulent  shear  flows.  This  leads 
us  to  hypothesize  that  all  turbulent  shear  flows  are  organized  to  a lesser 
or  greater  extent,  depending  on  the  topology  of  the  geometry,  in  a 
hierarchy  of  flows  that  places  the  organized  wake  of  the  circular  cylinder 
on  top. 
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It  should  be  emphasized  that  all  these  conclusions  were  made 
possible  by  the  fact  that  the  data  that  result  from  the  Laser  Induced 
Fluorescence  measurements  are  field  measurements  as  opposed  to 
point  measurements.  The  information  regarding  the  structure  of  the 
various  flow  patterns  that  were  observed  would  simply  not  have  been 
available  by  measurements  on  fewer  dimensions.  A preliminary  report 
of  some  of  the  early  conclusions  is  included  in  the  form  of  Appendix  H. 
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C.  Acoustical  Imaging 

Transmission  of  sound  waves  through  a given  flow  field  can  be 
utilized  to  extract  information  about  the  flow.  Since  a temperature 
change  is  the  principal  cause  in  changing  the  speed  of  sound,  trans- 
mission imaging  was  demonstrated  by  using  a steady  plume  as  the  object. 
A heater  was  mounted  at  the  bottom  of  a water  tank,  while  two  piezo- 
electric probes  were  rigidly  connected  and  traversed  through  the  tank 
(Fig.  1).  A pulse  generator  triggered  a function  generator  to  generate 
3 MHz  sine-wave  bursts  of  about  50  msec  duration,  which  were  fed  into 
the  transmitter  probe.  The  receiving  probe  picked  up  the  sound  bursts, 
which  were  then  amplified,  bandpass  filtered  and  averaged.  The  plotter 
recorded  this  signal,  which  is  a measure  of  the  transmitted  sound 
amplitude,  as  a function  of  the  traverse  position.  Figure  2 shows  such 
a plot  while  Figure  3 shows  the  actual  temperature  variation  through  the 
plume  as  measured  with  a thermocouple.  As  expected  from  ordinary 
shadowgraphy,  the  acoustic  trace  shows  the  second  derivative  of  the 
thermocouple  trace,  where  the  amplitude  decrease  in  the  center  is  due 
to  the  deflection  of  the  center  rays  out  from  the  center. 

In  order  to  investigate  unsteady  flows,  an  acoustical  sensing  array, 
comprised  of  100  piezoelectric  elements  in  a line,  has  been  designed  and 
fabricated.  Each  element  measures  .035"  x 0.  5"  to  give  a total  length 
of  3.  5".  Since  utilizing  sound  waves  give  the  opportunity  to  measure  also 
the  phase  of  the  signal  directly,  processing  electronics  are  under  design 
to  extract  amplitude  and  phase  from  the  incoming  signal,  which  will  then 
be  recorded  in  digital  form  ontape  as  a function  of  time.  The  recording 
speed  is  designed  to  be  500  frames/sec.  The  unique  feature  to  have 


amplitude  and  phase  available  at  the  same  time  allows  one  to  obtain 
various  images  from  the  same  recording,  e.g.,  amplitude  only  images 
and  phase  only  images.  It  is  also  possible  to  apply  various  numerical 
reconstruction  techniques  like  Fresnel  Transforming,  which  require 
the  knowledge  of  both  amplitude  and  phase. 

It  is  planned  to  apply  this  technique  in  a gas  flow  using  sound 
frequencies  of  about  350  kHz.  Since  the  speed  of  sound  is  a strong 
function  of  the  temperature  only,  an  acoustical  image  can  give  com- 
plementary information  to  the  optical  image.  Another  application  is  a 
flow  situation,  where  light  is  too  sensitive,  like  in  the  HF  mixing  layer 
facility. 


Figure  1.  SCHEMATIC 
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Appendix  A 

NUMERICAL  INVESTIGATION  OF  THE  COALESCENCE  OF  LARGE  SCALE 
STRUCTURES  IN  THE  TURBULENT  MIXING  LAYER 
D.  D.  Knight 

I.  Introduction 

During  the  past  twenty-five  years,  considerable  attention  has  been  devoted 

to  the  phenomenon  of  the  large-scale  structures  ("eddies")  in  turbulent  shear 

1 2 

flows.  The  early  work  of  Townsend  ’ and  his  colleagues  emphasized  the  importance 

of  the  large  eddies  in  the  dynamics  of  turbulent  flows.  Recently,  Brown  and 
3 

Roshko  discovered  the  remarkable  coherence  of  the  large  eddies  in  a two- 
dimensional  turbulent  mixing  layer.  Additional  evidence  of  coherent  structures 
has  been  found  in  a wide  variety  of  turbulent  shear  flows,  including  the  two- 
dimensional  wake,^  boundary  layer, ^ and  axi-symmetric  jet.**  These  results  have 

7 8 

led  to  new  theories  as  to  the  role  of  the  eddies  in  the  phenomena  of  entrainment,  ’ 
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mixing,  combustion  and  pollutant  formation. 

The  purpose  of  this  research  is  to  investigate  the  ability  of  turbulence 

11  12 

models  (in  particular,  the  model  equations  of  Saffman  * ) to  describe  the  features 

of  the  large-scale  structures  in  the  Incompressible  homogeneous  turbulent  mixing 
layer.  Specific  attention  is  focused  on  the  process  of  amalgamation  ("pairing") 
of  large  scale  eddies. 

II.  Definition  of  Problem 

The  problem  under  consideration  is  the  two-dimensional  temporally-developing 
turbulent  mixing  layer.  Indicated  in  Figure  la.  The  mixing  layer  is  formed  by  two 
parallel  streams  of  infinite  extent  in  the  x-dlrection.  The  flow  is  statistically 
unsteady,  since  the  width  of  the  mixing  region  grows  in  time.  In  contrast,  the 
typical  experimental  configuration  is  the  spatially  developing  mixing  layer  formed 
by  two  semi-infinite  parallel  streams,  as  Indicated  in  Figure  lb.  The  flow  is 
statistically  steady,  and  the  mixing  width  grows  linearly  downstream. 

The  temporally-developing  mixing  layer  is  considered  for  reasons  of 
computational  cost.  In  this  configuration,  the  mean  convective  velocity  of  the 
eddies  is  approximately  zero,  thereby  allowing  the  use  of  a fixed  computational 
domain  whose  streamwise  extent  is  twice  the  eddy  size  (for  the  case  of  a single 
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pairing"),  in  the  spatially  developing  mixing  layer,  the  mean  convective 
velocity  la  approximately  (U^  + U2)/2,  thereby  requiring  a computational  domain 
whose  atreamwlae  extent  is  several  times  the  typical  eddy  size  in  order  to 
consider  a similar  amalgamation. 

There  la  no  precise  transformation  between  the  temporally-  and  spatlally- 
<*®vc^°P*n8  cases.  Qualitatively,  however,  the  configurations  are  similar  if 
related  by  the  Galilean  transformation 

*’  - + U2)t  (!) 

The  freestream  velocity  in  the  temporally  developing  case  is  thus  defined  as 

«-  “ *<Ul  - u2>  (2) 

The  turbulent  motion  is  assumed  to  consist  of  a combination  of  large-scale 
coherent  motion  and  fine-scale  random  fluctuations.  Considering  the  instantaneous 
streamwlse  velocity,  U(x,y ,z,t),  a conditional  average  <(/>,  denoted  by  u for 
simplicity,  is  defined  by 


where  T is  defined  by: 
o 

tl*a  ecale  of  random  motion  <<T^<<  time  scale  of  coherent  motion. 

2 

Based  on  the  experimental  data,  the  conditionally-averaged  ("coherent") 
motion  la  assumed  to  be  two-dimensional  and  unsteady.  The  instantaneous  streamwlse 
velocity  is  thus 

U * u + u" 

where  u"  represents  the  fine-scale  random  fluctuations.  Similarly,  the  Instantaneous 
velocities  in  the  y and  z directions  are  defined  by 


Mb*' 
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V - v ♦ v" 

W - w + w" 

where  in  the  present  case  w ■ 0. 

The  conditionally-averaged  conservation  equations  are 
3u,  3u.  * 


Momentum 


a * *2yi\ 

■■  tt  ♦ “j  - - is,  * (-  <»r‘j->> + vs72 


Mass: 


3x . 


(A) 


(5) 


Xj  ■ x,  *2  • y and  v is  the  kinematic  viscosity.  The 


where  u^  • u,  Uj  * 

11  12 

turbulence  model  equations  of  Saffman  * provide  additional  rate  equations  for 

the  "fine-scale  Reynolds  stresses"  <ui"uj">,  thereby  yielding  a closure  of  the 

equations  of  motion.  For  purposes  of  brevity  the  model  equations  are  not 
“v. 

presented  here;  complete  details  are  given  in  references  11  and  12. 

The  initial  profiles  for  the  velocity,  pressure  and  fine-scale  Reynolds 
stresses  are  taken  to  be  the  similarity  solution  of  the  equations  of  motion, 
wherein  there  is  no  streamwlse  (i.e.,  x)  variation  (e.g.,  u - U^f(n)  , where 
n " y/VatQ  is  the  similarity  independent  variable  and  t#  is  an  arbitrary  initial 
time).  In  addition,  a perturbation  is  added  to  the  velocity  and  pressure  field. 

The  perturbation  is  the  sum  of  two  contributions:  1)  the  (approximately)  most 
unstable  eigenfunction  (harmonic)  corresponding  to  an  lnvlscld  (Raleigh)  stability 
analysis  of  the  initial  velocity  profile,  and  2)  the  first  subharmonic  [i.e.,  the 
eigenfunction  whose  wavelength  is  twice  that  of  the  disturbance  in  1)].  The 
relative  phase  of  the  two  contributions  is  chosen  to  yield  an  initial  vorticlty 
field  which  is  symmetric  upon  reflection  through  an  origin  located  at  the  center 
of  the  computational  region.  The  pair  of  eddies  which  initially  develop,  therefore, 
have  the  same  "strength,"  as  discussed  later.  The  kinetic  energy  of  the  initial 
perturbation^equi-partltloned  between  harmonic  and  subharmonic,  is  taken  to  be  4Z 
of  the  initial  kinetic  energy  defect  of  the  mixing  layer.  The  initial  Reynolds 
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based  on  the  velocity  difference  (2Uw),  velocity-profile  maximum  slope  thickness 
end  kinematic  viscosity  is  480. 

The  boundary  conditions  assume  periodicity  In  x over  the  length  of  the 
computet ional  domain  which  is  equal  to  the  wavelength  of  the  first  subharmonic. 
Freestream  boundary  conditions  are  taken  at  sufficiently  large  |y|  where 
u • 0^  and  all  fine-scale  Reynolds  stresses  vanish. 


HI.  Results 

Mean  Flow  Properties 

The  mean  velocity  u in  the  streamwlae  direction,  denoted  by  U,  is  defined  as 


U • u = ^ ^udx 


(6) 


where  L ■ length  of  computational  domain.  The  velocity-profile  maximum  slope 
thickness  6^(0  ('Vorticity  thickness")  Is  defined  by 


2U 


6 - 

01 


3U 

3y 


(7) 


Using  equations  (1)  and  (2),  the  computed  results  yield  a growth  rate  given  by 


dfi 


<»i  - v 


— - - 0 24 

dx'  (Uj  ♦ U2) 


The  experimental  best-fit  indicated  by  Brown  and  RoshkoJ  is 


d<5 


(ux  - u2) 


dx*  " °*10(U1  + U2) 


although  the  numerical  coefficient  is  dependent  upon  the  state  of  the  boundary 

13 

layer  at  the  end  of  the  splitter  plate  with  values  varying  from  0.16  in  the 

14 

experiment  of  Liepmann  and  Laufer  to  0.24  in  the  results  of  Wygnanski  and 


Fiedler . 


15 


In  view  of  the  only  approximate  similarity  between  the  temporally  and 


mm 


— ^ 
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spatially  developing  mixing  layers,  the  agreement  with  experiment  is  reasonable. 

In  Figure  2,  the  mean  velocity  profile  U(y,t)  is  shown  as  a function  of  y/8, 

g 

where  6 is  the  momentum  thickness  defined  in  the  usual  fashion.  In  this  figure, 
the  time  t has  been  non-dlmenslonallzed  using  U and  the  initial  half-width*  of 
the  mixing  layer,  with  t - 0 implying  the  initial  conditions.  It  is  apparent  that 
the  mean  velocity  profile  has  remained  self-similar  despite  substantial  growth  in 
the  coherent  perturbations.  In  particular,  define  the  coherent  perturbation 
kinetic  energy  by 


“coh  L 


f \ \ i(u*2  + v,2)dxdy 


where  u'  and  v'  are  given  by 
u'  • u - U,  v'  • v 

and  define  the  mean  kinetic  energy  defect  by 


- j HU.2  - U2)dy 


As  indicated  above,  the  initial  coherent  perturbation  kinetic  energy  is  4* 


of  the  initial  E 


At  t“  10,  however,  the  coherent  perturbation  kinetic 


energy  is  23X  of  the  current  mean  kinetic  energy  defect.  The  self-similarity 
of  the  mean  velocity  is  in  agreement  with  experiment.2 


The  initial  half-width  is  defined  as  the  value  of  y for  which  U ■ .99U 

a 

for  the  initial  velocity  profile. 
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Urw  Seal*  Structure 

The  evolution  and  interaction  of  a pair  of  large-scale  eddies  is  indicated 
in  figures  3a  through  31.  Bach  figure  is  a contour  plot*  at  a particular  tlae 
of  the  conditionally-averaged  vorticity  C where  c ■ ^ • Selected 

numerical  values  of  the  vorticity  are  shorn,  where  the  vorticity  has  been  non- 
dlaanslonallxed  using  and  the  initial  half-width  of  the  mixing  layer.  As 
Indicated  earlier,  the  particular  choice  of  the  relative  phase  of  the  harmonic 
and  subharmonlc  yields  initially  two  eddies  whose  maximum  absolute  vorticity  is 
identical  (Figure  3a). 

The  vorticity  contour  plots  clearly  indicated  the  growth.  Interaction  and 
eventual  amalgamation  or  "pairing"  of  the  two  eddies.  The  features  are 

8 

qualitatively  similar  to  the  experimental  dye  photographs  of  Wlnant  and  Browand, 

and  the  experimental  vorticity  contour  plots  of  the  large-scale  structures  of 

Browand  and  Weidman. ^ In  particular,  the  correlation  coefficient  of  the  large- 
st* 

scale  perturbations  at  t • 8 (see  Figure  3d)  is  greater  than  0.9  for 

■1.3  < y/0  < 1.3.  This  is  in  agreement  with  the  experimental  result  of  Browand 

and  Weldman^  who  found  the  correlation  coefficient  of  the  large  structures  to  be 

greater  than  0.9  for  -2  1 y/0  s 2 at  a stage  in  the  vortex  pairing  where  the 

*** 

vorticity  contours  display  a strong  resemblance  to  Figure  3d. 


*The  contour  plots  were  produced  on  a line  printer.  The  algorithm  employed 
divided  the  vorticity  at  a particular  time  into  a fixed  number  of  equal 
Increments  between  its  minimum  and  max'! mum  values.  The  projection  of  these 
levels  onto  the  x-y  plane  were  printed  alternately  as  characters  and  blanks. 

The  result  is  more  precisely  a "level"  plot,  rather  than  a contour  plot;  for 
example,  the  dark  region  at  the  center  of  the  eddies  in  Figure  3a  corresponds 
to  -1.68  < C < -1.58,  and  the  immediately  surrounding  blank  annular  region 
corresponds  to  -1.58  < c < -1.47.  For  purposes  of  simplicity,  the  numerical 
values  indicated  in  Figures  3a  to  31  are  the  mid-range  values  (l.e.,  -1.63, 

-1.52,  for  the  above  case). 

•♦Defined  by  T • - uTvr/(ur2r  v1"2)^ . Note  definitions  of  equations  (6)  and  (9). 

***This  is  Stage  II  of  Browand  and  Weidman1^  (Figure  6b).  In  comparing  the  figures, 
due  account  must  be  given  to  the  differences  in  definition  of  direction  of 
mean  flow. 


The  qualitative  features  of  fluid  entrainment  by  the  large-scale  structures 
have  been  examined  by  inclusion  of  a passive  scalar  (i.e.,  diffusion)  equation.  The 
initial  passive  scalar  (e.g.,  fluid  "dye")  profile  was  located  slightly  above  the 
mixing  layer  with  a uniform  concentration  in  x and  a Gaussian  distribution  in  y. 

The  "dye"  concentration  at  suceeding  times  is  indicated  in  Figures  4a  through  4i , 
where  three  arbitrary  fixed  levels  of  concentration  have  been  plotted  with 
increasing  darkness  in  order  to  give  a qualitative  visual  impression  of  the  dye 
concentration.  At  early  times,  both  the  diffusion  of  the  dye  by  fine  scale 
turbulence  and  the  convection  by  the  large-scale  structures  is  apparent.  As  the 
eddies  amalgamate,  a large  fraction  of  the  dye-bearing  fluid  is  entrained  into 
the  mixing  layer,  and  the  region  of  maximum  concentration  eventually  coincides 
with  the  final  amalgamated  eddy.  It  is  interesting  that  the  dyed  fluid  appears  to 
be  swept  across  the  mixing  region  before  eventually  appearing  as  a concentrated 
eddy-like  structure. 

IV.  Conclusions 

The  dynamical  evolution  of  the  large-scale  structures  in  the  temporally- 
developing  mixing  layer  have  been  studied  numerically  using  the  turbulence  model 
equations  of  Saffman.  Preliminary  results  Indicate  quantitative  agreement  with 
mean  growth  rate  and  self-similarity  of  the  mean  velocity  profile.  The  results 
display  the  observed  phenomena  of  "pairing"  of  large-scale  structures  in  qualitative 
agreement  with  experiment.  Further  work  is  currently  being  pursued  to  provide 
more  detailed  comparison  of  the  features  of  the  large-scale  interactions. 
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Figure  2:  MEAN  VELOCITY  PROFILE  AT  SEVERAL 
STAGES  DURING  THE  EDDY  PAIRING 


Figure  3a:  VORTICITY  CONTOURS:  t - 2 

Values  of  vortlcity  are  shown.  The  top  and 
bottom  row  of  asterisks  indicate  vertical 
extent  of  the  computational  domain. 
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Figure  3c:  VORTICITY  CONTOUR 


Figure  3d:  VORTICITY  CONTOUR: 
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VORTICITY  CONTOUR 


Figure  3g:  VORTICITY  CONTOUR 
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Figure  31:  V0RTIC1TY  CONTOUR 
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Figure  4a:  DYE  CONCENTRATION:  t - 2 

The  edges  of  the  mixing  layer  y « + 6 
(where  U - + .99U  ) are  shown  by  arrows. 
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Figure  4c:  DYE  CONCENTRATION: 
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Figure  4d:  DYE  CONCENTRATION 
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Appendix  B 


On  the  Approximate  Motion  of  Vorticity- Bearing  Eddies 
in  a Perfect  Fluid  in  Three  Dimensions 


by 

Doyle  D.  Knight 

Applied  Mathematics  Department 
California  Institute  of  Technology 
Pasadena,  California  91125 


The  approximate  equations  of  motion  of  a dilute  collection 
of  vorticity- bearing  eddies  in  an  unbounded  perfect  fluid  in  three 
dimensions  are  derived.  The  planar  interaction  of  a pair  of 
eddies  is  considered  numerically  for  two  initial  configurations. 
Application  to  a deterministic  model  of  the  large  eddies  in  a 
turbulent  shear  layer  is  discussed. 
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I.  Introduction 

The  past  several  years  have  seen  a resurgence  of  interest 
in  the  role  of  the  "large  eddies"  in  turbulent  flows.  First  utilized 
extensively  by  Townsend  in  his  research  on  turbulent  shear  flows 
(Townsend  1956),  the  concept  of  the  large  eddies  has  recently  been 
the  subject  of  considerable  research  (for  a brief  review,  see 
Roshko  (1976)). 

The  present  results  are  a first  step  in  an  approximate 
deterministic  model  of  the  large  eddies  in  a turbulent  flow. 

The  "large  eddies"  are  viewed  as  simply- connected  regions 
or  "blobs"  of  fluid  with  non- zero  vorticity  in  an  otherwise  perfect 
fluid  in  three-dimensions.  They  form  a dynamical  system,  whose 
approximate  equations  of  motion  are  described  below.  The  effects 
of  viscosity  (i.e.,  entrainment)  and  generation  are  currently  under 
study. 

II.  Formulation  of  Problem 

We  consider  a collection  of  simply- connected  regions  of 
non- zero  vorticity  (the  "large  eddies")  in  an  otherwise  perfect 
fluid  in  three-dimensions.  For  simplicity,  we  assume  that  the 
flow  is  unbounded,  at  rest  at  infinity,  and  contains  no  sources, 
sinks  or  body  forces.  (For  further  discussion,  see  IV). 

As  the  fluid  is  assumed  perfect,  Kelvin's  theorem  allows 

that  vortex  tubes  move  with  the  flow,  and  thus  the  volume  of 

rotational  fluid  comprising  each  eddy  remains  constant.  We 

. tlx 

define  the  center  of  mass  R of  the  i n eddy  with  volume 
V.  in  the  usual  fashion, 

l ’ 
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where  x is  the  position  vector.  It  is  straightforward  to  show 
that 


V 


dRl 
i dt 


= I*  + f (p  n dS 
9V. 


(2) 


where 


- Impulse  of  i4^  eddy  = ■—  x A u dV 


V. 
x 

w = vorticity  = curl  u 

u - velocity 

q>  = velocity  potential 

n = outward  normal  on  surface  of  i^  eddy 
As  the  fluid  is  assumed  unbounded  and  at  rest  at  infinity, 
the  velocity  field  can  be  expressed  as 


*) = - h 5 

All  space 
.th 


(x  - x')  A o(xr,  t)dV(x') 


X - X 


I I s 


(3) 


At  some  point  within  the  i eddy,  it  is  convenient  to  represent 
the  velocity  as  a sum  of  two  terms, 

= 2e  + Si 

where  u_(x,t)  is  the  contribution  due  to  the  presence  of  the 
eddies  exterior  to  the  i**1  eddy  (i.e.,  the  volume  integral  in  (3) 
is  divided  into  the  region  exterior  and  the  region  interior  to  the 
i**1  eddy).  It  is  then  possible  to  derive  the  dynamical  result 
(Saffman,  1971)* 


dl* 

•dT 


= /.  “e  a £ dv* 


(4) 


t Professor  Saffman  has  noted  that  the  essential  idea  of  equation 
(4)  is  evident  in  Betz  (1933). 
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The  system  of  equations  (2)  and  (4)  provide  the  frame- 
work for  an  approximate  description  of  the  motion  of  the  eddies. 


III.  Approximate  Equations  of  Motion 

We  first  consider  the  equation  for  the  impulse  (4).  By 
utilizing  the  fact  that  in  the  vicinity  of  the  i**1  eddy  u^.  may 
be  expressed  as  the  gradient  of  a scalar,  and  expanding  u_  in 
a Taylor’s  series  about  the  center  of  mass  of  the  eddy,  we 

can  obtain  after  some  manipulation 

dl 1 

-ar  s - A(t)Il  (5) 

where  A(t)  is  a second  order  tensor  given  by 


A(t)  = grad  u (R  (t),t) 

i*. 

and  contributions  from  higher  order  terms  in  the  Taylor  series 
expansion  are  ignored.  The  velocity  field  at  the  center  of  mass 
of  the  i**1  eddy  due  to  the  presence  of  all  other  eddies  (i.e., 
Uj.(R1(t),  t))  can  be  expressed  using  the  asymptotic  form 
(Batchelor,  1970), 


'"Tt F.X  grad{  I^(t)  . grad  ( )} 

£1 


x = R (t) 


where  the  summation  omits  the  i n eddy. 

It  is  possible  to  give  a physical  interpretation  of  the  term 
on  the  right  side  of  equation  (5),  analogous  to  the  "vortex  stretching" 
effect  in  the  vorticity  equation  (Batchelor,  1970).  The  quantity 
- A(t)I^  represents  the  negative  of  the  product  of  two  terms: 

1)  The  magnitude  of  the  impulse  (i.e.,  |ll|  ),  and  2)  The 
relative  rate  of  change  of  an  infinitesimal  material  line  element 
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vector  at  R (connecting  two  points  along  the  direction  of  I1)  due  to 
the  irrotational  straining  field  induced  by  the  other  eddies 

We  next  consider  the  equation  for  the  center  of  mass  (2).  We 
divide  the  velocity  potential  in  the  surface  integral  over  the  i**1  eddy 
into  two  contributions,  arising  from  the  vorticity  exterior  and  interior 
to  the  i*^  eddy: 

V ~ <?£.  + 

where  u^,  = grad  We  then  obtain 

/ ^ n ds  = / uEdV  = u (Rl(t),t)Vl 


where  higher  order  (quadratic)  terms  arising  from  the  Taylor  series 
expansion  of  u_  are  omitted.  Utilizing  the  asymptotic  expansion  for 


we  obtain 


C 1 ,i 

J e nds  = - 3 I 


where  higher  order  terms  arising  from  the  non- spherical  nature  of  the 
eddy  surface  are  omitted.  The  result  is 


vi^*Tl  + vi5E<5  <7> 

where  u^,  is  given  by  (6). 

For  a collection  of  N eddies  in  an  unbounded  quiescent  fluid, 
equations  (5)  and  (7)  represent  a closed  set  of  6N  equations  for  the 
vectoral  quantities  R1,  I1,  i = 1,  N.  There  is  in  addition  the  global 
requirement  that  the  total  impulse  be  conserved.  This  condition  can 
be  easily  shown  to  follow  from  equations  (5),  (6)  and  (7)  for  N = 1 or  2, 
and  it  is  expected  that  it  can  be  shown  for  arbitrary  N. 

It  is  trivial  to  verify  that  equations  (5)  and  (7)  reduce  to  the 
exact  description  of  a single  Hill  spherical  vortex  in  an  unbounded 
flow,  i.e.. 
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The  basic  assumption  inherent  in  the  above  approximations 
is  that  the  collection  of  eddies  is  "dilute",  i.e.,  that  their  average 
separation  is  much  larger  than  their  typical  size. 


III.  Interaction  of  Two  Eddies 

For  the  purpose  of  providing  a simple  qualitative  picture, 
the  two-dimensional  interaction  of  two  eddies  of  equal  volume  was 
considered.  The  governing  equations  were  solved  numerically  using 
a fourth-order  Runge-Kutta  scheme  (Isaacson  and  Keller,'  1966). 

The  accuracy  of  the  method  was  verified  by  computing  the  behavior 
with  the  impulse  of  one  of  the  eddies  set  identically  to  zero.  The 
motion  of  the  "eddy"  with  zero  impulse  is  equivalent  to  that  of  a 
fluid  parcel  in  irrotational  flow  about  a sphere,  thus  providing  a 
check  on  the  numerics. 

Figure  1 shows  the  motion  of  a pair  of  eddies  whose  impulses 
at  t = 0 were  equal  in  magnitude  and  opposite  in  sign,  directed 
along  the  x-axis  (i.e.,  horixontal).  The  positions  of  the  eddies  are 
shown  at  equal  time  increments,  beginning  with  t = 0.  A notable 
feature  is  the  development  of  a component  of  impulse  in  the  vertical 
direction  by  each  eddy,  as  evidenced  by  the  inclined  trajectories  at 
large  time.  As  the  total  impulse  is  conserved,  the  final  vertical 
velocities  are  of  equal  magnitude  and  opposite  sign. 

An  altogether  different  behavior  is  indicated  in  Figure  2, 


which  shows  the  motion  of  two  eddies  whose  initial  impulses  were 
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+ 1 and  +3,  directed  along  the  Ix-axis.  Again,  the  effect  of  the 
straining  field  produced  by  one  eddy  on  the  other  is  evident  in  the 
development  of  a component  of  impulse  in  the  vertical  direction. 

By  use  of  a Galilean  transformation,  the  initial  velocities  of  the 
eddies  (though  not  the  impulses)  in  Figure  2 can  be  made  approx- 
imately equivalent  to  those  of  Figure  1,  and  the  straining  field  ex- 
perienced by,  for  example,  eddy  #1  is  qualitatively  similar  in  both 
figures.  The  observed  difference  between  the  two  cases  in  direction 
of  vertical  motion  of  eddy  =1  at  large  time  can  thus  be  seen  to  result 
from  the  difference  in  orientation  of  its  impulse  (see  equation  (5)). 

IV.  Extensions 

The  eddy  model  as  described  above  is  as  yet  incomplete. 
Some  of  the  additional  effects  which  require  incorporation  are  in- 
dicated below. 

a.  Entrainment.  The  model,  being  inviscid,  provides 
no  means  for  entrainment  of  irrotational  fluid  by  the 
large  eddies.  The  turbulent  entrainment  process 
must  be  modeled,  as  it  is  the  basic  growth  mechanism 
of  a turbulent  shear  layer. 

b.  Generation.  For  turbulent  boundary  layers  on  smooth 
surfaces,  the  mechanism  by  which  the  vorticity  gener- 
ated at  the  surface  is  'ejected'  into  the  boundary  layer 
is  a complicated  phenomena  (see;  for  example,  Kline 
et  al,  1967).  This  generation  cycle  for  the  eddies  re- 
quires additional  modeling. 

There  are  additional  questions  of  interest,  in  particular, 


the  possibility  of  extending  the  approximate  equations  of  motion 
by  considering  higher  order  terms,  thereby  introducing  tensor 
descriptions  of  the  vorticity  distribution  within  each  eddy. 
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Figure  1:  Motion  of  pair  of  eddies  whose  impulses  are  of  opposite  sign. 

Positions  of  each  eddy  are  shown  at  equal  time  increments. 
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THE  TURBULENT  CHEMICAL  REACTION 
OF  A HYDROGEN  JET  DISCHARGED  INTO  FLUORINE 
Frank  E.  Marble,  J.  E.  Broadwell, 

O.  P.  Norton,  M.  V.  Subbaiah 
ABSTRACT 

The  coherent  flame  model  for  turbulent  chemical  reactions  des- 
cribes the  flame  structure  in  terms  of  a distribution  of  laminar  flame 
elements,  each  of  which  retains  a flame-like  behavior  while  it  is 
being  stretched  in  its  own  plane  and  transported  by  the  turbulent 
field.  In  its  simplest  form,  which  is  employed  here,  it  should 
describe  fast  chemical  reactions  in  the  early  stages  of  mixing  with 
some  accuracy. 

In  the  present  work  this  model  is  applied  to  the  turbulent  hy- 
drogen-fluorine reaction  with  reference  to  the  H-F  laser.  Because 
the  non- equilibrium  distribution  of  vibrational  states  of  the  HF  mole- 
cule is  central  to  the  laser  problem,  flame  strain  rates  £ must  be 
considered  such  that  £ ^ ^ / where  is  the  relevant 

chemical  time.  This  requirement  necessitates  modification  to  the 
consumption  rate  of  reactants  that  is  utilized  in  the  model.  For 
these  high  strain  rates,  the  flame  is  no  longer  diffusion  controlled  with 
a relatively  small  zone  of  chemical  reaction.  Rather,  the  chemistry 
occupies  the  entire  diffusion  zone.  The  model  has  been  modified  to 
accommodate  this  new  description  of  reactant  consumption  and  the 
mechanics  of  calculating  the  distribution  of  vibrationally  existed  states 
has  been  developed. 

Because  the  description  of  the  turbulent  jet  and  the  chemistry 
of  the  hydrogen-fluorine  reaction  are  fixed,  the  only  parameters 


i 


required  to  define  the  problem  are  the  initial  states  of  the  hydrogen 
and  fluorine  streams  and  the  value  of  V\/*  3c  / <L  (sometimes 

referred  to  as  a Damkohler  number)  where  VJo  is  the  initial  velocity 
of  the  jet  issuing  from  a port  of  diameter  d.  , and  is  a character- 
istic chemical  time  of  the  reaction.  Detailed  distribution  of  the  HF 
product  species  and  the  corresponding  populations  of  vibrational  states 
have  been  calculated  for  two  values  of  V/»  Tc  / cL  . These 

results  are  given  at  distances  of  2,  7 and  15  jet  diameters  downstream 
of  the  port.  At  similar  locations  within  the  jet,  the  larger  value  of 
U/c  Te  / d,  tends  to  accentuate  the  non- equilibrium  population 

of  vibrational  states.  The  effect  of  distance  along  the  jet  is  more 
complex.  It  is  shown  that  the  combined  effects  of  flame  surface 
density  and  strain  rate  produce  a maximum  in  the  vibrational  popu- 
lation non-equilibirum  that  occurs  several  diameters  downstream  of 
the  jet  origin. 
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1.  INTRODUCTION 

The  coherent  flame  front  model*1*  is  a description  of  fast 
chemical  reactions  in  turbulent  flow  in  which  the  reactions  are 
assumed  confined  to  thin  flame  surfaces.  The  turbulent  flame 
structure  then  consists  of  a distribution  of  these  surface  elements. 

The  model  describes  the  manner  in  which  the  flame  elements  are 
stretched  and  the  distribution  is  dispersed  by  the  turbulent  motion, 
as  well  as  the  mechanism  by  which  neighboring  flame  surfaces  con- 
sume the  intervening  reactant  and  annihilate  each  other.  A consider- 
able advantage  of  the  coherent  flame  model  is  that  it  effectively  sep- 
arates the  detailed  structure  of  the  laminar  flames  from  the  fluid 
mechanics  so  that  systems  with  complex  chemical  reactions  may  be 
treated  nearly  as  easily  as  simple  ones. 

The  work  described  here  treats  a turbulent  combustion  process 
associated  with  the  hydrogen-fluorine  CW  laser.  The  geometry  of 
the  problem  consists  of  a circular  jet  of  hydrogen  gas  discharging 
into  a stationary  field  of  atomic  fluorine.  Both  gases  are  diluted 
with  helium.  Although  this  configuration  differs  somewhat  from  that 
which  occurs  in  practical  lasers,  it  serves  to  illustrate  the  mechanics 
of  applying  the  coherent  flame  model  in  a situation  where  the  assump- 
tions of  the  model  are  quite  closely  satisfied. 

The  central  issue  is  the  determination  of  the  degree  to  which  the 
vibrational  degrees  of  freedom  for  the  HF  molecule  differ,  during  the 
turbulent  combustion  process,  from  their  normal  equilibrium  distri- 
bution as  well  as  the  distribution  of  these  states  over  the  jet  flame. 
The  complexity  of  the  chemistry  required  to  describe  even  a some- 
what simplified  model  of  the  process  demonstrates  quite  clearly  the 
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ad vantages  that  accrue  to  the  coherent  flame  description. 

The  present  application  of  the  flame  sheet  theory  differs  in 
one  essential  way  from  the  previous  developments,  references  1,  Z. 

In  earlier  applications  we  have  employed  the  quasi-steady  strained 
diffusion  flame  result  that,  for  a diffusion-controlled  flame,  the  re- 
actant consumption  rate  per  unit  flame  area  is  proportional  to 
l/X>£  where  & is  the  coefficient  of  molecular  diffusion  and  £ is 
the  strain  rate  (sec  *)  in  the  plane  of  the  flame  surface.  In  this 
circumstance  the  actual  region  of  chemical  reaction  is  thin  in  compar- 
ison with  the  thickness  of  the  diffusion  zones.  In  other  words,  the 
characteristic  chemical  time  3*.  is  short  in  comparison  with  the 
value  of  t _1;  hence  £ Tc  <■  <■  / . But  it  is  almost  evident  that 

in  order  to  disturb  the  normal  distribution  of  HF  vibrational  states, 
the  fluid  dynamical  time  must  be  short  compared  with  the  relevant 
chemical  time,  so  that£Tc>>l  . When  this  circumstance  prevails, 
the  reactant  consumption  is  no  longer  diffusion  controlled  and,  indeed, 
the  entire  diffusion  zone  is  chemically  active.  Then  the  reactant 
consumption  rate  is  proportional  to  the  thickness  of  the  diffusion  zone, 
v^T  divided  by  the  chemical  time,  or  Jj>/£  . Thus  a 

local  reactant  consumption  law  has  been  introduced  that  encompasses 
both  the  diffusion  controlled  and  chemically  controlled  flame.  The 
same  law  is  appropriate  to  describing  the  distribution  of  vibrationally 
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excited  states. 
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2.  THE  FLAME  STRUCTURE  MODEL 

To  make  this  turbulent  flame  picture  quantitative,  define 

the  flame  surface  density  £ to  be  the  flame  surface  area  per  unit 

volume;  the  dimensions  of  the  flame  surface  are  very  large  in 

comparison  with  its  thickness.  The  local  flame  surface  density  is 

altered  by:  1)  production  resulting  from  fluid  straining  motions, 

2)  turbulent  transport,  and  3)  flame  annihilation,  in  which 

neighboring  flame  surfaces  consume  the  intervening  fuel  or  oxidizer. 

The  structure  of  the  flame  surface  elements  is  dominated 

by  the  straining  motion  in  the  plane  of  the  flame.  Then  if 

2 

£ is  the  straining  rate  and  D - (length)  /time  is 

the  molecular  diffusion  coefficient,  the  flame  thickness  is 
and  consumes  reactants  at  a rate  "-\/P£  per  unit  flame  surface  area. 

Now  when  the  straining  rate  is  very  high  so  that  £ T«. 
where  is  a characteristic  chemical  reaction  time,  the  reactant 

consumption  rate  is  proportional  to  the  thickness  of  the  diffusion 
zones  of  the  flame,  because  it  is  within  this  layer  that  the  reactants 
are  molecularly  mixed.  Then  the  reactant  consumption  is***  ■^'\/p''£ 
The  representation  of  the  law  which  encompasses  both  diffusion 
controlled  and  chemically  controlled  reactant  consumption  is 

(it 

It  £?<. 
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When  the  flame  surface  density  Iwand  the  strain  rate  € are  known 


locally,  the  reactant  consumption  rates  within  this  region  are  known 


also  because  they  are  fixed  by  the  flame  structure.  The  actual  flame 


structure  may  be  determined  to  any  degree  of  approximation  between 


the  simplest,  where  the  chemical  reaction  rates  are  very  rapid  so 


that  the  reactant  consumption  is  controlled  by  molecular  diffusion,  to 
the  detailed  structure  involving  the  complete  chemistry  with  appro- 


priate rates.  The  flame  structure  with  finite  rates  requires  numerical 
calculation  but,  because  it  is  a steady  one-dimensional  problem,  the 


actual  calculations  may  be  accomplished  with  relative  ease. 


The  appropriate  form  of  the  equation  describing  the  flame  surface 


density  may  be  deduced  from  first  principles  by  considering  the  dis- 


tortion and  migration  of  surface  elements,  fixed  to  the  fluid,  in  a 


(3) 

turbulent  medium.  Here,  it  will  be  sufficient  to  motivate  the  form 


by  physical  reasoning  and  then  to  suggest  the  manner  in  which  the 
various  terms  scale  with  features  of  the  flowfield  and  the  chemistry 


of  the  flame  structure.  Now  in  a fluid  with  mean  velocity  components 


, the  expression 


M , IS,  „ 


gives  the  change  in  flame  density  following  a mean  fluid  element. 


According  to  our  model  described  earlier,  this  change  may  be  written 


in  the  following  form 


- turbulent  diffusion  of  flame  elements  into  the  region 


+ increase  in  flame  surface  element  area  by  turbulent 
straining  motions 
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- reduction  in  flame  surface  resulting  from  local  consump- 
tion of  the  reactants. 

In  these  calculations,  the  turbulent  diffusion  of  flame  surface  will  be 
described  using  a turbulent  diffusivity  and  the  assumption  that  the 
turbulent  fluctuation  velocities  are  large  in  comparison  with  diffusion 
velocities,  familiar  from  usual  treatments  of  turbulent  heat  and  mass 
transport.  Thus,  if  we  denote  J5"  the  diffusion  coefficient  that  arises 
in  the  description  of  momentum  exchange,  the  turbulent  diffusion  term 
applicable  to  the  treatment  of  the  turbulent  circular  jet  is 


! 


i 
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where  the  appropriate  boundary  layer  approximation  has  been  made. 

The  rate  of  increase  of  an  element  of  flame  surface  area  is  proportional 
to  the  strain  rate  in  the  turbulent  fluctuations.  Under  the  assumption 
that  this  strain  rate  is  proportional  to  the  rate  of  strain  in  the  mean 
motions,  the  rate  of  increase  in  flame  surface  density  is  thus  pro- 
portional to  the  product  of  the  strain  rate  of  mean  motion  and  the  local 
flame  density.  Thus  we  take  the  second  term  on  the  right  equal  to 

c £ 5 * /fr/  <2 

where,  again,  we  have  written  the  term  in  the  form  appropriate  to 
the  circular  jet,  W being  the  velocity  component  of  the  mean  motion 
in  the  direction  of  the  symmetry  axis. 

The  process  by  which  flame  surface  is  removed  from  the  field 
is  best  pictured  by  considering  two  neighboring  laminar  diffusion 
flame  fronts  parallel  to  each  other  and  containing  one  constituent, 
say,  fuel,  between  them.  As  the  motion  progresses,  the  intervening 
fuel  is  consumed  and  both  elements  of  flame  surface  are  extinguished. 


t 


P 
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j 
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A similar  process  takes  place  if  the  intervening  constituent  is  the 
oxidizer.  To  make  the  mechanism  quantitative,  we  note  that  fraction 
of  local  volume  occupied  by  fuel  is  proportional  to  the  concentration 
K,  of  fuel.  Moreover,  the  rates  of  reactant  consumption  by  an 
element  of  laminar  flame  are  presumed  known  from  detailed  calcula- 
tions of  that  flame  structure.  If  we  call  V,  the  volume  rate  of 
consumption  of  fuel  per  unit  flame  area,  the  rate  at  which  volume 
occupied  by  fuel  is  being  consumed  is  l/(  £ so  that  the  fractional 
rate  of  fuel  consumption  in  a unit  volume  of  space  is  V,  . Thus, 

if  the  flame  surface  area  is  nearly  uniformly  distributed  over  the  region, 
this  expression  gives  also  the  fraction  of  flame  surface  that  is  annihi- 
lated so  that  the  flame  surface  area  reduction  due  to  fuel  consumption 

(¥)«  ' 

In  a completely  similar  manner,  we  may  reason  that  the  flame  surface 
area  reduction  due  to  consumption  of  oxidizer  is  proportional  to 


and  because  in  a given  region,  the  surface  is  being  reduced  by  either 
the  exhaustion  of  fuel  or  oxidizer,  but  not  both,  the  two  expressions 
will  be  considered  additive. 


If  now  we  collect  the  various  terms  that  have  been  discussed, 
it  is  possible  to  write  in  detail  the  conservation  equation  for  flame 
surface  density 
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The  form  of  equation  has  been  chosen  to  be  that  appropriate  for  the 
analysis  of  the  circular  jet. 

The  conservation  equations  for  the  individual  reactants  are  con- 
ventional except  for  the  terms  describing  the  reactant  consumption  by 
chemical  reaction.  These  will  be  given  in  terms  of  the  flame  surface 
density  and  the  reactant  consumption  for  a unit  flame  surface  supplied 
by  one-dimensional  flame  calculations.  For  the  fast  overall  chemical 
reaction,  we  need  consider  only  three  constituents,  the  fuel  k,  , 
the  oxidizer  K*  , and  the  product  . related  by  the  fact  that  the 
sum  of  mass  fractions  is  unity 


K,  *»  K,  ■ » 


(4) 


so  that  only  two  of  these  K,  and  k»  , need  be  calculated.  The 
detailed  chemistry,  in  contrast  to  the  overall  reaction  between  fuel 
and  oxidizer,  is  contained  within  the  one-dimensional  flame  structure. 

The  consumption  rate  of  fuel  per  unit  volume  is  simply  the  pro- 
duct of  the  effective  influx  velocity  V,  , the  flame  surface  density  £>  , 
and  the  mass  fraction  K,"*  of  fuel  in  the  fuel- containing  component 
in  the  unmixed  state.  The  conservation  equation  for  fuel  may 

be  written  down  directly  as 


_-  -l  A.  £ 

y Jr  vaj  r dr  V 3*"  1 


(5) 


Similarly,  the  conservation  equation  for  the  oxidizer  component  is 


(6) 


It  remains  to  define  the  consumption  rates  V,  and  V*.  for  the 
fuel  and  oxidizer  components.  Regardless  of  whether  these  quantities 
are  defined  through  use  of  infinitely  fast  reaction  rates  or  by  detailed 
calculation  of  the  one-dimensional  flame  structure,  they  depend  upon 
the  reactant  concentrations  on  each  side  of  the  flame,  and  , 

and  the  local  strain  rate  £ of  the  mean  flow.  It  is  this  latter  item 
which  couples  the  local  diffusion  flame  structure  to  the  gasdynamic 
structure. 


-11- 


3.  THE  TURBULENT  CIRCULAR  FUEL  JET 
The  detailed  solution  for  a turbulent  flame  structure  described 
in  the  manner  outlined  in  the  previous  section  requires  a turbulence 
model  for  the  mean  flow  and,  because  there  is  a choice  of  turbulence 
model  to  be  used,  it  is  preferable  to  go  directly  to  the  problem  of 
interest  rather  than  describe  the  procedure  in  general.  For  the  pro- 
blem of  the  circular  fuel  jet,  we  shall  choose  the  elementary  model 
utilizing  a scalar  turbulent  diffusivity.  Furthermore,  we  shall,  in 
the  interests  of  simplicity,  neglect  the  change  of  mean  gas  density 
associated  with  the  chemical  reactions.  This  restriction  in  no  way 
implies  that  the  density  change  is  a negligible  factor.  Rather  the 
use  of  a turbulence  model  for  flows  of  non-uniform  density  introduces 
a degree  of  uncertainty  of  its  own  which  makes  it  additionally  diffi- 
cult to  judge  the  merits  of  the  flame  model. 

Under  these  restrictions,  the  gasdynamic  field  is  described  by 
the  equations 


JL  A 
r dr 


(rU)*2^ 


-r  O 


(7) 


j / + i /S>y.  -J-JLfr  _£r*  ) 

u dr  ^d*  r dr  l f J 


(8) 


where  the  pressure  of  the  far  field  has  been  assumed  constant.  Using 
the  eddy-diffusivity  model  for  the  turbulent  shear  stress  , we  write 


(9) 


where  i®*  is  the  turbulent  analogue  to  the  kinematic  viscosity.  The 
customary  difficulty  arises  here  in  the  choice,  rather  than  the  deter- 
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mination,  of  the  eddy  diffusivity.  For  the  circular  jet,  however,  the 
choice  is  made  simpler  by  the  fact  that  the  jet  momentum,  which  is 

2 

conserved  along  the  direction  of  flow,  has  the  dimension  of  (velocity) 
(length)  . And  because^  has  the  dimension  of  (length)  x (velocity), 
there  being  no  other  characteristic  dimensions  in  the  problem,  it  is 
evident  that  the  turbulent  diffusivity  is  proportional  to  the  square  root 
of  the  jet  momentum,  and  hence  is  a constant. 

It  is  possible  to  find  a similarity  solution  for  the  turbulent  jet, 
or  more  properly  for  a point  momentum  source.  To  carry  this  out, 
it  proves  convenient  to  introduce  a modified  radial  dimension 


1.  * 


(10) 


where is  the  constant  jet  momentum  flux 


•/ 


l*/x  Pair 


The  stream  function  , with  the  usual  properties  that 
W * , may  then  be  written  in  the  form 


(11) 


u*-i 


24 

3K 


p ->  fry)  (12) 

Substitution  into  the  momentum  equation,  Equation  8,  yields  the 
ordinary  differential  equation  for  F<y\ 

which  may  be  integrated  once  to  give 

. _L  c r ’ - f " L F ' 

(13) 

the  constant  of  integration  being  zero  because  yF  and  f vanish  for 
large  “f  . Fortunately,  the  solution  for  Equation  13  may  be  written 
simply  as 


x 
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F(*n  • 


V 


which  satisfies  the  requirement  that  the  stream  function  vanish  on  the 
symmetry  axis,  * o , and  satisfies  identically  the  momentum  con- 

servation relation.  Equation  11,  as 


/(V)* 


7^7 


The  velocity  components  of  the  mean  flowfield  are  then 

. , . _L.  4*  -A-  / Jl ) 

£ J_ t 

* 4 » * 0* 


- tK-p‘) 

sr  j. 

vr/*  * (,-;>>)*• 


Now  explicit  solution  of  the  flame  density  and  reactant  conservation 
relations  require  a quantitative  form  for  the  reactant  consumption  rates. 

If  the  chemistry  were  actually  of  infinitely  fast  rate  and  if  the  mutual 
reactants  were  uncontaminated  by  products  of  combustion,  the  values 
of  vt  and  Vi,  appearing  in  Equations  5 and  6 could  be  written 


K/V,  * p,  \/l>£ 


; = p*. 


where  V is  the  coefficient  of  molecular  diffusion,  proportional 

to  the  turbulent  strain  rate,  and  (3,  , ^9*.  are  constants  for  a given 

chemical  system,  depending  upon  the  reaction  stoichiometry  and  the 
amount  of  diluent  accompanying  the  reacting  species  in  the  unmixed 
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ga*  stream*.  Numerical  calculations  of  strained  laminar  diffusion 

flames,  which  will  be  described  subsequently,  showed  that  strain  rates 
3-1  7-1 

between  10  sec  and  10  sec  were  required  to  achieve  a large  alteration 

in  the  distribution  of  vibrational  states  of  H-F.  These  high  strain  rates 

were  accompanied  by  a distinct  change  in  the  laminar  flame  structure. 

f>  - 1 

When  the  strain  rate  wa b less  than  10  sec  the  reactant 
consumption  rate  behaved  as  if  the  chemistry  were  infinitely  fast. 

In  this  range  the  reactant  consumption  rate  is  ~v/D/  where  P is  the 

appropriate  binary  diffusion  coefficient.  When  the  strain  rate  was 

*>  - 1 - Zx 

greater  than  10  sec  , the  reactant  consumption  varied  as  £ , 

the  results  showing  that  the  reactant  consumption  rate  was  chemically 

limited  and  that  the  chemically  active  zone  filled  the  entire  flame 

thickness.  It  was  clear  then  that  the  proper  expression  for  the 

reactant  consumption  rate  at  large  strain  rates  is 


(19) 


where  v/p/f  is  proportional  to  the  flame  thickness  and  is  a chemical 

time. 

For  the  entile  range  of  strain  rates  the  reactant  consumption 
per  unit  flame  area  is 

x/pT 

(Z0) 


x*  Vi  * /*; 


IS  7*  £ 


where  p.  * (i,  , depending  upon  whether  we  are  considering  the 

consumption  of  hydrogen  or  fluorine  atoms.  This  gives  .i  rate 
a ' for  strain  rates  £3,**/  and  (ft  / for 

larger  strain  rates  £ > > / . The  equations  for  reactant  eoniump- 


«nd  flame  density  may  then  be  written 
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\A>  C 


}r  w 2 i f >»“'  9r  J l'|-»3t£  ^ 


.,?Kx.t,)Kt  . JL  JL.  ( r f;  ] - fl  ---  - £ 

U T7  wTit  ' rarl  r*  .*  J*  t *“* 


^<vySC.  . -I  j-(r>2iL)- 

) »,  3 i c 5*  ’ Sr''  * i + 7+  l 


K,  Kt 


K,  -t  l<v  ♦ * v - < (24) 

When  the  chemistry  is  fast,  as  it  is  in  this  case,  it  is  convenient 

to  use  a linear  combination  of  K.  and  Kj  that  satisfies  a homogeneous 

i 

diffusion  equation.  Such  a combination  is  K,  “ (*Z7p,  S ^ 

and  because  a diffusion  flame  with  fast  chemistry  utilizes  fuel  and 
oxidizer  in  the  stoichiometric  ratio  ^ , the  quantity  and 


K,  -* 


Thus  J(r,  z)  satisfies 


5i  ^ . A . LJ  .-  -L  2.  / r fi  J ) 


with  the  conditions  *®and  hence,  comparing  Equations  26,  8,  and 

9,  it  is  clear  that  - W/t*,*)  . The  constant  of  proportionality 

between  J and  W is  found  by  noting  that  if  the  actual  flow  of  fuel  in- 
jected into  the  stream  is  X , then  the  flow  of  the  product  formed  is 
(i-frvflir  for  large  values  of  z where  the  reaction  is  complete  and 
Ki  = ° . Then  the  integrated  flux  of  the  product  is,  as  ? 


AT,  U/  rdLr  * X ( lf*f ) 


9 [~ 

and  because  integration  of  Equation  26  shows  that  jJWr'ct*  is  oon- 
served  along  the  z direction,  it  follows  that 
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o» 

f JW  rdr  * r 

Jo  (28) 

Because  J is  proportional  to  W,  the  integral  is  the  same  as  evaluated 
to  ensure  conservation  of  momentum  with  z and  hence  it  follows,  com- 
paring Equations  11  and  28,  that 


J • 


- A.  JL 

« tr 


(29) 


We  have  then  the  algebraic  relation 


in  addition  to  the  relation  Equation  24  and,  as  a consequence,  we  need 
only  work  with  a single  species  conservation  Equation  21. 

In  order  to  make  this  identification  between  J(r,  z)  and  W(r,  z),  it 
has  been  necessary  to  introduce  the  volume  flow  rate  X as  well  as  the 
momentum y*.  of  the  jet.  These  two  quantities  may  be  used  to  char- 
acterize the  jet  velocity 


velocity  — 

and  the  characteristic  dimension  of  the  jet. 


(31) 


lenqih  ~ 


(32) 


where,  in  fact,  the  length  — *\Zrt‘R,  where  R is  the  initial  radius 

Y/*- 

of  a jet  with  uniform  velocity. 
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4.  THE  STRAINED  LAMINAR 
HYDROGEN-FLUORINE  FLAME 

To  complete  the  detailed  description  of  the  model,  we  require  the 
numerical  solution  for  the  strained  hydrogen-fluorine  flame  which  de- 
termine the  flame  structure,  including  the  distribution  of  vibrationally 
excited  states,  and  the  reactant  consumption  rates.  Results  of  this 
type  can  be  obtained  analytically  when  the  chemistry  is  simple  but  in 
the  present  case  in  which  the  vibrational  populations  result  from  a 
complex  interaction  of  chemical  kinetics  and  molecular  diffusion,  a 
numerical  calculation  is  necessary.  These  were  made  with  the  Blottner 

boundary  layer  computer  program  modified  so  as  to  apply  to  the  laminar 
(4) 

diffusion  flame'  The  boundary  layer  equations  and  the  numerical 
scheme  employed  in  their  solution  are  described  in  references  5 and 
6.  The  program  modifications  together  with  a comparison  of  results 
with  experimentally  determined  HF  vibrational  population  profiles  are 
given  in  reference  7.  The  following  brief  description  of  the  computer 
program  is  taken  from  the  latter  reference. 

The  physical  situation  is  sketched  in  Figure  1.  Two-dimensional 
streams  of  hydrogen  and  fluorine  (diluted  with  helium)  flowing  from 
opposite  sides  of  a splitter  plate  interdiffuse  and  react  to  form  the 
laminar  flame  sheet,  in  which  the  various  vibrational  levels  of  hy- 
drogen fluoride,  HF(v),  are  produced.  When  the  free  stream  reactant 
concentration  are  uniform  and  the  stream  velocities  are  constant  and 
equal  to  Vj<o  , the  flame  is  equivalent  to  the  "standard"  time  dependent 
diffusion  flame  (with  replaced  by  t );  the  strained  flame  is 

generated  when  is  caused  to  increase  linearly  with  z. 

If  the  chemistry  is  simple  and  the  reaction  rate  fast,  it  is 


1 J 


i 


>. 
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1 


straight-forward  to  prove^  that  the  structure  and  reactant  con- 
sumption rates  of  a flame  experiencing  a strain  rate  € are  identical 
to  those  of  an  unstrained  time  dependent  flame  at  the  time  t * ‘/>te  , 

A few  comparisons  of  our  numerical  results  for  H*.'  F flames  in  the 
two  cases  indicated  that  the  two  solutions  are  similar  when  the  same 
transformation  of  independent  variable  is  made.  Consequently,  the 
shorter  time  dependent  computation  was  used  to  obtain  the  results 
described  in  the  following. 

The  Blottner  program,  which  solves  the  standard  boundary 
layer  conservation  equations,  was  designed  especially  to  treat  many 
chemical  species  and  chemical  reactions  and  to  allow  multi- component 
diffusion.  In  the  present  case,  the  kinetic  system  consisted  of  the 
"pumping"  reactions: 

Ht  * F — *■  Hf(\r)  + H : \r*  / .2. 3 

1 (33) 

together  with  equations  describing  the  important  deactivational  effects 
of  H,  F as  well  as  HF  itself.  The  set  of  reaction  equations  to- 
gether with  the  associated  rate  coefficients  are  given  subsequently. 

The  temperature  dependence  of  the  species  viscosities  are  program 
inputs.  The  other  transport  processes  were  treated  with  the  simplify- 
ing assumption  that  the  Prandtl  and  Lewis  numbers  were  constant  at 
unity. 

Both  reactant  streams  were  at  300°K  and  5 torr.  The  fuel  stream 

8 3 7 

composition  was:  hydrogen,  2.67*10  g mole/cm  , helium  2.47-10 

3 

g mole /cm  . The  composition  of  the  oxidizer  was  atomic  fluorine, 

*8  3 - 7 3 

2.67*10  g mole/cm  , helium  2.47*10  g mole/cm5. 


The  program  outputs  consist  of  various  property  profiles  of 
temperature  and  concentration  for  instance,  and  of  H,  and  F consumption 

u 


rates.  The  dependence  of  the  latter  on  strain  rate,  £ , is  shown  in 

Figure  2.  As  the  arguments  presented  in  earlier  sections  lead  us  to 

expect,  the  consumption  rates  at  large  values  of  the  strain  vary  like 
-Ji  14 

£ and  at  low  values  like  £.  * . Equation  20,  given  previously. 


K i * 


P'  /+  7c  £ 


is  a fit  to  the  curve  in  Figure  2.  The  numerical  constants  required 
for  our  work  are 


2c  - S.IJT*  to  iec  (34) 

The  vibrationally  excited  and  ground  state  populations  within 

the  flame  structure  depend,  as  do  the  reactant  consumption  rates, 

upon  the  rate  of  strain.  These  populations,  in  the  form  of  integrated 

mass  (g  moles)  per  unit  flame  area  are  shown  in  Figure  3.  We  note 

the  large  departures  from  equilibrium  at  the  high  strain  rates,  indeed 

r -i 

total  inversions  exist  for  £>loiec  . Empirical  approximation  to 
the  results  in  Figure  3 that  are  convenient  for  the  turbulent  jet 
analysis  are  given  in  Section  6. 

The  "complete"  chemical  reaction  system  for  the  H-,  - re- 
action as  presented  in  reference  8 consists  of  some  100  chemical 
reactions.  For  given  ranges  of  temperature  and  concentration,  many 
of  these  reactions  have  negligible  influence  on  the  HF  vibrational 
populations  and  need  not  be  considered.  For  purposes  of  the  present 
investigation  it  was  decided  that  the  important  excited  state  pumping 
and  deactivational  processes  were  contained  in  the  eighteen  equations 
given  in  Table  1.  Also  shown  are  the  coefficients  in  the  expression 
for  the  forward  and  backward  reaction  rate  coefficients  which  are  of 


V.- 


the  form: 


Af  « Ca  7"c‘  <f*p  /oooC./t) 
Ak>  - d,tt>x  0KJ3  ('■too © Z>,/t) 
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5.  INTEGRAL  FORMULATION  OF  HYDROGEN  JET 

IN  FLUORINE 

To  treat  the  jet  of  hydrogen  issuing  into  a stationary  background 
of  fluorine  atoms,  consider  equations  21,  23  and  30.  If  there  were 
no  reactions  and  if  the  flame  straining  were  negligible  ( ot *o  ),  then 
the  problem  would  possess  similarity  solutions  of  the  form 
and  £ . There  is  a certain  advantage  in  writing  these 

quantities  in  that  form,  as 


A2,  - 'e'  & ( “i  • & ) 


(35) 


and 


£ 


, \/jL  i 

i-&) 


(3,2 


«■  C 7 .£  ) 


(36) 


Now  because  we  shall  use  an  integral  technique  to  obtain  our 
final  result,  it  is  convenient  to  rewrite  Equations  21  and  23  in 
similarity  variables  and  integrate  them  over  the  flame  cross  section. 
Introducing  the  similarity  coordinates,  *>£  and  z,  we  have 

7>K, 


& 


-±i>  (,»§*)-  t'VZZh 


l*¥l 


(3  7) 


and 


£ 


/3wi  , /_hJi \ 


(38) 


Integrate  these  from  7*0  to  ^ and  take  account  of  the  behavior  of 
A2,  and  ^5  for  large  radii.  We  then  have  from  Equation  3 7 


— I ( a \ * f ) - 


<r(7,07  ^7 


(39) 
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where  we  have  used  a new  variable  for  distance  along  the  axis, 

/ = jS.  % 

S 3 r 

and  introduced  a Damkohler  number 

r * t/$~  * -Tr  7‘ 

Treating  the  flame  density  relation  in  the  same  manner 


(40) 


(41) 


«U 


»(i) 


& i 


ydiy 

l '+  ■&£(-&  *' 


(42) 


where  in  this  notation 


*x  - 1 ~ ~j~  ( 


(43) 


To  complete  the  solution  using  the  integral  technique,  it  is  re- 
quired to  select  profile  shapes  for  the  mass  fraction  distribution  and 
for  the  flame  surface  density.  These  are  chosen  specifically  in  the 
following  form: 


■4,(7, n * (t  ) ( 1 * 
kxC7.U  --  ' ‘ (y  )f  (‘* 

U1,H  * 


(44) 


(45) 


(46) 
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where  -fti)  and  ini  f ) , unknown  functions  of  distance  along  the  axis, 
will  be  determined  through  the  integral  relations.  Equations  39  and  42. 
This  choice  posesses  several  features  in  its  favor: 

a)  They  are  exact  when  there  is  no  chemical  reaction. 

b)  They  satisfy  exactly  the  conditions  given  by  Equation 
30  among  species  mass  fractions. 

c)  They  satisfy  the  appropriate  boundary  conditions  at  the 
symmetry  axis  and  at  distances  remote  from  the  axis. 

d)  The  flame  surface  density  vanishes,  as  it  must  physically, 
where  either  the  fuel  or  oxidizer  reactants  vanish. 

Now  the  value  of  Ftiy) « £ 7\)  » Equation  14,  is  the  exact 

solution  for  the  stream  function  in  the  jet  problem  utilizing  the 

turbulent  diffusivity  formulation.  Utilizing  this  and  the  representations 

for  it,  <”7,1)  » $ ) given  in  Equations  44,  45,  and  46,  the 

integrations  indicated  in  the  integral  relations,  Equations  3 9 and  42, 

may  be  carried  out.  The  integrals  are  of  three  classes;  (i)  those 

which  may  be  evaluated  directly,  (ii)  those  which  may  be  reduced  to 

an  integral  of  and  have  the  value 


/ ^ 

/ o**rr 

o 

and  (iii)  those  of  the  form 


IT  ( "n'Vzl ) 

1 * \ l(nq)  ' 


. [ til?  Air.£J 


This  last  class  may,  after  a transformation,  be  written  in  the  form 


-- 


xn*V»  . »/♦ 

jt L »-*) 

±JL  (,-*)' 


which  is  suitable  for  numerical  integration. 

Now  if  we  utilize  these  results  in  the  integral  relations. 


1 


I 
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Equations  39  and  42  , for  the  hydrogen  concentration  and  the 

flame  density  respectively,  then  these  become  respectively 


o-n 


f/+ 

- (-£-)  A 


(\\o-n 


( t ; t ) - 


and 

• Upfi-wX-tHiuO-ol'-ii  i-LitL'>] 

V*  r ) 

“ 1 (t)  <«;»>} 


(51) 

These  are  to  be  integrated  from  f * 2.  , the  virtual  origin  of  the  jet 

being  at  i * o , to  an  arbitrary  value  of  I ; the  initial  values  of 
■f  l - •>  and  III.)  ■ A.  , where  gives  a measure  of  the 

amount  of  flame  surface  that  is  present  at  the  start  of  calculation. 
Physically,  the  initial  value  of  (\U.)  is  related  to  an  ignition  process 
which  the  turbulent  flame  model  does  not,  of  course,  describe.  It  is 
to  be  expected  that  varying  the  initial  value  will  have  a significant 

effect  upon  the  structure  of  the  flame  i * l0  but  a much  smaller  effect 
away  from  the  origin.  The  values  of  and  — • are  fixed  by  the 
reactants  and  by  the  turbulence  model,  respectively.  The  values  of  tx. 
and  A , on  the  other  hand,  are  numerical  constants  which  are  universal 
in  the  sense  that  they  do  not  depend  upon  stoichiometry,  momentum  flux, 
or  geometrical  size  of  the  hydrogen  jet. 
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6.  DISTRIBUTION  OF  REACTION  PRODUCTS 
The  results  that  are  desired  from  our  calculations  are 
the  distribution  of  the  product  HF  molecule  within  the  jet  and  the 
distribution  of  this  product  over  the  vibrational  states  that  are  im- 

i 

portant  in  the  laser  operation. 

First,  the  distribution  of  HF  mass  fraction  follows  from  the 
hydrogen  mass  fraction,  which  we  have  just  determined 


where  is  mass  of  HF^  per  unit  area  of  flame  surface,  is 

the  maximum  value  of  this  quantity.  T ^ is  the  characteristic  time 
(reciprocal  of  strain  rate)  at  which  this  maximum  occurs  and  6^  is 
i constant  giving  the  breadth  of  the  distribution.  As  discussed  in 
section  4,  these  numbers  follow  from  the  calculations  of  figure  3, 


o 

II 

> 

v = 1 

v = 2 

v = 3 

c(i) 

2 x 10'9 

o 

H 

1 

o 

r-H 

X 

r- 

5 x IQ’10 

IQ’10 

7.  46xl0-6 

4.38  x 10'5 

2.  46  x 10*4 

1.  38  x 10'3 

25 

10 

4.  17 

1.  67 

Table  2.  Constants  for  Representation  of 
Distribution  of  Vibrational  States,  Eq.  54 

Now  to  evaluate  the  distribution  of  states,  explicit  calculations 

of  strain  rate  and  the  actual  flame  surface  density  are  required.  The 

strain  rate  then  becomes 

*IttI  * ) 7*  (</->  ) (55) 

and  recalling  the  definition  of  T from  Equation  41, 

fT‘  * 1 "e  o^m’  (56) 

and  hence  the  quantity  £ 1 required  for  the  evaluation  of  is 


€ 7* 


where  the  second  factor  is  given  in  Table  1. 

— <i) 

The  distribution  of  vibrational  states,  C*.  , is  then 


where  , the  flame  surface  density,  may  be  conveniently  written 


£ ••  - ( - zfz  \ )J-  y<u 

4 P.  („DJ  - ( j ( > > ' Y ( 

The  value  of  Y<.t\  has  already  been  calculated  and  hence  it  is  a 
simple  matter  to  compute  the  normalized  densities  of  vibrational 


states. 


where 


£ /j.  (_*£.)*-£. 

(i.  ' ' * 


It  should  be  noted  that  although  the  normalized  population  densities  are 
convenient,  in  that  they  remove  some  of  the  system  parameters  from 
the  results,  they  do  not  provide  a direct  scaling  analysis  for  the 
hydrogen  fluorine  system.  The  reason  for  this  is  simply  that  the 
system  parameters  are  involved  also  in  the  parameter  T,  which 
enters  Equation  56  for  the  representation  of  £ 7C  and  hence  in  the 
values  of  S^. 


• - U HI 
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7.  RESULTS  OF  SPECIFIC  CALCULATIONS 

We  shall  show  two  specific  cases  of  a hydrogen  jet  reacting 
with  a stationary  atmosphere  of  fluorine  molecules.  These  were 
made  with  the  specific  purpose  to  explore  the  behavior  of  the  model 
for  a particular  kinetic  scheme  rather  than  to  describe  the  distri- 
bution of  chemical  species  in  the  jet  for  a wide  range  of  jet  momenta 
and  geometric  scales. 

The  calculations  divide  naturally  into  two  steps:  (i)  Determination 
of  the  centerline  reactant  concentration  and  flame  surface  density, 
and  and  (ii)  Determination  of  the  distribution  of  chemical 

species  within  the  flame  element  and,  hence,  over  the  jet  cross 
section. 

The  solution  for  the  first  step  follows  from  equations  50  and  51 
which  are  rewritten  here  as 


oU 

tfC  t 


A ( i ; 7 ) 


(62) 


dLV 

c*.  1 


/«V>T\  B(l)  y _ T can)  yL 
{ & ) \ >-< 


(63) 


i ; 

• ■ 

i 


where  now 
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To  carry  out  this  calculation  we  require  values  for  the  several  con- 
stants that  arise  in  the  equations.  The  dimensionless  group 
is  determined  empirically  by  the  known  structure  of  the  non- reacting 
jet  and  eC  is  universal  constant  which  relates  the  local  turbulent 
strain  rate  to  the  strain  rate  j ^ J of  the  mean  flow.  In  fixing 
the  value  of  e(  from  jet  measurements  it  is  as  easy  to  fix 
directly,  and  this  appears  to  be  approximately 


= r 


The  value  of  X , which  occurs  in  connection  with  the  flame 
shortening  terms  of  Equation  63,  has  been  chosen  on  the  basis  of 
the  experiments  of  Hawthorne  et  al^  to  be 

A - Z (69) 

The  integrals./^  and  J,  , that  occur  in  the  coefficient  A(ii  7)  and 
C ( i;  3)  contain  the  ratio  of  chemical  time  Jc  » to  the  time  defined 
by  the  jet  momentum  and  flow,^t  and  X . Referring  to  Equation  41, 
the  value  of  7 is  proportional  to 


•v  o( 


(w.y) 

( w. 
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where  W*  and  are  the  initial  jet  velocity  and  jet  diameter.  As  the 
value  of  y changes  from  small  values  ) to  large  values,  the 
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local  diffusion  flame  elements  change  from  diffusion  controlled  to 
chemistry  controlled.  And  from  the  standpoint  of  non- equilibrium 
populations  of  vibrational  states  of  HF,  the  larger  values  of  T favor 
non- equilibrium. 

The  solutions  for  (1-f)  and  Y , proportional  to  the  centerline 
values  of  hydrogen  mass  fractions  and  flame  density  respectively,  are 
shown  in  Figure  4 for  7-  f and  in  Figure  10  for  7 - a ® . The 
general  character  of  the  flame  density  is  that  the  straining  process 
dominates  in  the  early  portion  of  the  jet,  increasing  the  flame  den- 
sity until  the  flame  annihilation  terms  become  important,  stopping  the 
rise  in  flame  density.  After  a distance  of  10-12  port  diameters,  the 
flame  straining  and  flame  annihilation  processes  are  nearly  in  balance, 
but  lead  to  the  gradual  decrease  in  flame  surface  area  as  the  hydrogen 
concentration  diminishes.  The  effect  of  is  evident  and 

shows  that  for  large  T values.  Figure  10,  the  flame  annihilation 
mechanism  is  somewhat  slower  in  retarding  the  growth  in  flame  sur- 
face. This  results  from  the  higher  strain  rates  involved  which  lead 
to  a reactant  consumption  rate  ~ -L.  \/v/£  f giving  lower 

values  than  for  smaller  strain  rates.  Thu6  the  active  flame  occupies 
a smaller  fraction  of  the  volume  for  Ui  Jc  / large  and  hence  post- 
pones the  peak  in  flame  density  until  it  has  reached  larger  values  of 

i ~ */<L 

Figure  5-9  summarize  the  distribution  of  reaction  products  over 
the  jet  for  three  distances  downstream  from  the  jet  port  corresponding 
to  values  of  = 2 , 7,  15.  Figure  5 gives  the  radial  distribution  of 

the  ground  state.  Figure  5a,  and  three  vibrationally  excited  states 


■ 
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of  the  HF  molecule  at  a value  of  $ = 2.  The  general  features  of 
the  distribution  are  similar  for  all  vibrational  states  but  the  concentra- 


tion,  as  indicated  on  the  vertical  scale,  shows  the  rapid  decrease  for 


higher  vibrational  levels.  It  should  be  noted  that  the  vertical  scales 


of  each  figure  have  been  adjusted  so  as  to  give  a reasonable  geometric 
representation  of  the  distribution;  the  actual  magnitudes  of  the  pop- 
ulations must  be  obtained  with  reference  to  the  scales.  Figures  6 


and  7 show  a similar  set  of  results  at  ( = 7.  While  the  relative 


populations  of  the  HF  vibrational  levels  is  changed  only  slightly,  the 


absolute  values  of  the  Cv  are  increased  by  a factor  of  nearly  3.  The 
corresponding  results  for  { = 15,  given  in  8 and  9,  show  again  a 
similar  relative  population  of  vibrational  states  but  the  absolute 
values  have  declined  to  about  half  their  magnitude  at  t = 2.  This  be- 


havior, which  was  described  earlier,  results  from  the  fact  that  both 


flame  density  and  strain  rate  enter  the  final  result.  The  strain  rate 


decreases  along  the  jet  length  while  the  flame  density  increases  at 
first  and  then  decreases.  In  the  present  example  for  7=5,  Figures 
4-9,  the  flame  density  variation  is  the  more  significant  factor. 


Two  other  features  should  be  kept  in  mind  when  looking  at 


these  results.  First,  the  actual  volume  that  contains  molecules  in 


any  state  along  the  radius  is  proportional  to  the  radius  and  hence  the 
high  values  near  the  jet  axis  carry  little  weight,  the  volume  weighted 
maximum  occurring  between  values  of  ^ = 1 and  2.  Second,  the 
physical  radius  r is  proportional  to  y * and  hence  increases 
linearly  with  distance  along  the  jet  which,  again,  influences  the 
distribution  along  the  jet  discussed  above. 


Figures  10  through  15  give  a corresponding  description  of  the 
distribution  of  hydrogen,  flame  density  and  HF  vibrational  states  for 
a value  y = 20,  in  contrast  with  the  value  of  7 = 5 for  the  previous 
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series.  Inasmuch  as  the  chemical  parameters  remain  the  same,  this 
increase  in  the  value  of  J may  be  considered,  referring  to  Equation 
70,  as  an  increase  in  jet  velocity  and/or  a decrease  in  jet  radius. 
Because  it  is  the  comparison  with  the  series  for  7*5"  that  is  of 
interest,  the  specific  details  of  the  curves  will  not  be  discussed 
extensively. 

One  feature,  however,  that  is  particularly  apparent  in  Figures 
11a  - lid,  is  the  alteration  of  radial  distribution.  This  feature 
results  from  the  dependence  of  strain  rate  upon  radius  and  is  accen- 
tuated for  larger  values  of  7/Sl  . Thus  it  appears  in  Figure  11  for 
i • Z but  not  in  the  distributions  at  larger  distances  along  the  jet 
axis. 

For  the  larger  value  of  7 , the  flame  surface  annihilation  rate 
is  considerably  reduced  over  that  at  7 * JT  and  hence  the  net  pro- 
duction of  flame  surface  is  somewhat  larger.  But,  on  the  other  hand, 
the  actual  consumption  of  reactants  is  reduced.  This  fact  is  reflected 
in  Figure  12  where  the  mass  fraction  of  HF  is  shown  to  be  somewhat 
less  than  that  for  the  slower  (7-  S’)  jet.  The  integrated  species 
in  various  degrees  of  vibrational  excitation,  however,  is  somewhat 
larger.  This  demonstrates  that  the  considerable  flame  surface  manu- 
factured in  the  range  1 t If  7 must  wait  until  the  strain  rate  has 
diminished  to  the  point  where  the  integrated  density  of  vibrational 
states  (see  Figure  3)  is  high.  This  tendency  to  delay  the  diffusional 
thickening  of  the  flame  sheet,  and  hence  the  appearance  of  reaction 
products,  is  one  of  the  unusual  effects  of  high  velocity,  high  y jets. 
By  the  time  the  value  of  $ * / »“  has  been  reached,  the  populations  of 
the  excited  vibrational  levels  appear  to  have  returned  reasonably 
closely  to  the  values  for  7*3"  . This  relative  independence  of  7 
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Figure  1.  Configuration  for  Calculation 
of  Strained  Diffusion  Flames 
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Figure  4.  Hydrogen  Concentration  (1-f)  and 
Flame  Density  (Y)  Along  Axis  of  Hydrogen 
Jet  in  Fluorine.  cp  =2.0,  0 ij  y/fi  =5, 

X =2,  T = 5.  0. 


Figure  5.  Population  of  HF  Vibrational  Level; 
t = 5,  C = 2. 
a)  Ground  State 
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Figure  5.  Population  of  HF  Vibrational  Level ; 
r = 5,  C = 2. 

c)  Second  Vibrational  Level 
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Figure  7.  Population  of  HF  Vibrational  Levels ; 

t = 5,  C = 7. 

c)  Second  Vibrational  Level 


Figure  8.  Distribution  of  HF  Mass  Fraction;  t = 5, 
C = 15 
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Figure  9.  Population  of  HF  Vibrational  Levels 
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Figure  9.  Population  of  HF  Vibrational  Levels 
r = 5,  C = 15 

b)  First  Vibrational  Level 
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Figure  10.  Hydrogen  Concentration  (1-f)  and 
Flame  Density  (Y)  Along  Axis  of  Hydrogen 
Jet  in  Fluorine,  cp  = 2,  a V |J  = 5,  \ = 2, 
T = 20. 
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Figure  11.  Distribution  of  HF  Vibrational 
Levels ; t = 20,  £ = 2. 
a)  Ground  State 
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Figure  13.  Distribution  of  HF  Vibrational  Levels; 
T = 20,  C = 7. 
a)  Ground  State 
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; 13.  Distribution  of  HF  Vibrational  Levels 
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Figure  15.  Distribution  of  HF  Vibrational 
Levels;  T = 20,  C = 15. 
b)  First  Vibrational  Level 
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1.  Introduction 


Recent  measurements  by  Johnson  and  Rose  (1973,  1975),  Yanta  and  Lee  (1974, 
1976),  and  Abblss  (1976)  have  used  laser  Doppler  velocimetry  techniques  to  make 
direct  measurements  of  the  Reynolds  stress  in  turbulent  boundary  layers  in  the 
Mach  number  range  of  1.5  to  3.0.  A serious  anomaly,  however,  is  exhibited  by 
these  measurements  in  that  the  maximum  of  -u'v'  occurs  much  further  from  the 
wall  than  is  reasonable  for  flow  at  constant  pressure.  The  anomaly  has  been 
discussed  by  Sandbom  (1974),  who  supports  the  conjecture  by  some  of  the  au- 
thors cited  that  density  fluctuations  may  contribute  substantially  to  the  tur- 
bulent stresses  near  the  wall.  This  conjecture  is  in  direct  opposition  to  the 
conclusion  by  Morkovin  (1961)  that  effects  of  density  fluctuations  should  be 
small  compared  to  effects  of  variations  in  mean  density  for  Mach  numbers  up  to 
4 or  5. 

The  purpose  of  the  present  experiments  is  to  obtain  redundant  data  over 
a substantial  range  of  Mach  numbers  (M  **  0.1  to  2.2),  in  an  effort  to  resolve 
the  anomaly  in  turbulent  shearing  stress.  The  flow  was  documented  by  conven- 
tional means  using  a Pitot  tube,  which  was  traversed  through  the  boundary  lay- 
er, to  measure  the  mean  flow.  In  addition,  surface-friction  measurements  were 
made  using  both  a floating-element  balance  and  a Preston  tube.  The  mean-flow 
scaling  suggested  by  Van  Driest  was  applied  to  the  data,  to  test  the  adequacy 
of  a single  similarity  formulation  for  both  compressible  and  Incompressible 
flow.  The  shearing  stress  was  then  computed  from  the  mean  flow  as  part  of  the 
analysis.  These  results  were  documented  by  Collins,  Coles  and  Hicks  (1978). 

In  the  work  presented  here,  the  instrumentation,  data  acquisition  and  ana- 
lysis for  laser  Doppler  velocity  measurements  in  the  same  flow  are  discussed. 
Measurements  of  u,  v,  u' z , v' z and  the  Reynolds  stress  u'v'  were  made  over  the 
full  Mach  number  range . A 

The  low  speed  (M  ^ 0.1)  measurements  were  carried  out  in  the  GALCIT 
Merrill  wind  tunnel.  The  high  speed  measurements  (0.6  M -£  2.2)  were  con- 
ducted in  the  JPL  20-inch  wind  tunnel. 

The  present  paper  is  a preliminary  report  of  a more  extensive  discussion 
of  this  work  to  appear  by  Dimotakis,  Collins  and  Lang  (1978). 
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2.  High  speed  laser  Doppler  velocimetry 

The  first  obvious  problem  that  arises  In  making  high  speed  laser 

^Graduate  Aeronautical  Laboratories  of  the  California  Institute  of  Technology 
^Jet  Propulsion  Laboratory 


Doppler  velocity  measurements  Is  that  the  resulting  signal  frequencies  are 
very  high.  In  principle.  It  would  appear  that  this  problem  could  be  controlled 
by  Increasing  the  fringe  spacing.  In  practice,  however,  conflicting  require- 
ments of  a minimum  number  of  Doppler  cycles  needed  by  the  processing  electron- 
ics, coupled  with  the  spatial  resolution  considerations  Inescapably  dictate 
the  higher  frequencies. 

The  restrictions  on  signal  frequencies,  above  and  beyond  the  problem  of 
handling  the  high  frequencies  per  se,  arise  for  several  reasons.  The  relative 
accuracy  6u/u,  with  which  the  velocity  of  a single  particle  can  be  measured, 

Is  never  less  than  the  ratio  of  the  processor  clock  period  Tc,  to  the  total 
flight  time  At  that  is  used  for  the  measurement,  i.e. 

6u/u  Xc/At  . (1) 

In  addition  to  these  considerations  which  limit  the  accuracy  even  if  the 
measurement  environment  was  noise  free,  other  factors  become  important  at  high 
velocities  which  decrease  the  signal-to-noise  ratio  of  a single  reading.  As 
the  flow  velocities  become  higher,  we  are  forced  to  use  smaller  particles  to 
minimize  problems  of  particle  lag  with  rapidly  decreasing  scattering  cross- 
sections  and  resulting  signal  intensities.  In  addition,  as  the  velocity  in- 
creases, the  scattering  particle  spends  less  time  in  the  focal  volume.  Even 
though  the  number  of  scattered  photons  per  second  remains  constant  (propor- 
tional to  scattering  cross-section) , fewer  total  photons  are  collected  by  tie 
detecting  optics. 

The  extent  to  which  the  measurement  of  particle  velocity  represents  the 
fluid  velocity  is,  of  course,  a separate  issue.  Small  particles  can  generally 
be  expected  to  track  the  flow  quite  well  if  the  (Lagrangian)  frequency  of  fluc- 
tuations in  the  fluid  velocity  is  less  than  the  reciprocal  of  a characteristic 
particle  response  time  Tp.  If  the  flow  relative  to  the  particle  can  be  de- 
scribed by  Stokes  flow  we  have 

Tp  - (dp/18) (pp/u)  - (dp/18v) (Pp/pf ) , (2) 

where  dp  is  the  particle  diameter,  pp  and  Pf  are  the  particle  and  fluid  densi- 
ties, respectively,  and  u is  the  (absolute)  viscosity.  It  can  be  argued,  and 
substantiated  from  direct  spectral  measurements  (Klebanoff  1955,  Perry  and 
Abell  1975,  1977),  that  the  expected  fluctuation  frequency  in  the  interior  of 
the  boundary  layer  should  be  given  by 

<w>  = const[U(y)/y] , (3) 

where  U(y)  is  the  local  streamwise  mean  velocity  and  where  the  constant  is  of 
the  order  of  unity.  Consequently,  the  requirement  for  good  particle  tracking 
in  a turbulent  boundary  layer  becomes 

Tp[U(y)/y]  4 1.  (4) 

By  way  of  example,  using  a dibutyl  phthalate  aerosol,  a flow  of  M ^ 1 requires 

particles  4 1.5  pm  diameter  to  track  the  fluctuations  at  y/6  £ 0.1.  For  a more 
detailed  discussion,  see  Dimotakis,  Collins  and  Lang  (1978). 


3.  Scattering  particles 

The  particle  generator  for  the  present  experiments  utilized  a Laskin 
nozzle  type  construction  to  generate  a polydisperse  aerosol  of  dibutyl 
phthalate.  Particles  greater  than  1 ym  in  diameter  were  effectively  removed 
by  an  Impact  plate  Incorporated  in  the  design.  The  resulting  particle  size 


distribution,  measured  using  both 

a cascade  impactor  and  a multichan-  4.0 - T ^ Co.code  Impactor 

nel  particle  analyzer,  is  shown  in  i 

figure  1.  35-  M I Multi *Chonnel Analyzer 

Two  different  methods  were  j j 

used  to  introduce  particles  into  , i'-o- 

the  flow  depending  on  the  flow  fa-  , 

cility  that  was  used.  For  the  high  « ■' 

speed  measurements  (0.6  M ^ 2.2),  ,§2.5-' 

which  were  conducted  in  the  JPL  20"  w I 1 

wind  tunnel,  it  was  necessary  to  £20-1  \| 

seed  the  flow  by  Introducing  the  I i | 

particles  into  the  settling  chamber  # t , 

between  the  last  turbulence  screen  “ 15  I Pl 1 1 J 

and  the  contraction  section.  The  5 mJI 

particles  were  introduced  through  a 1.0  * ^ 

tube  which  protruded  5 cm  into  the  M 

settling  chamber.  Holes  were  05_  — JjL 

drilled  into  the  tube  along  the  llj 

stagnation  line  as  well  as  the  rear  39^ 

side  at  ±30°  with  respect  to  the 

flow  vector.  0 05  10  15  20  2 5 50 

For  the  low  speed  measure-  Porticle  Oiometer  ( /*m) 

ments  in  the  Merrill  wind  tunnel 

(M  = 0.1),  although  some  naturally  Fig*  1.  Particle  size  distribution 

occurring  particles  were  present, 

the  flow  was  seeded  using  the  same  aerosol  generator  in  order  to  increase  the 
data  rate  and  to  better  control  the  scattering  particle  size  distribution. 

For  these  measurements,  the  particles  were  injected  through  a tube  spanning 
the  test  section  located  downstream  of  the  boundary  layer  plate.  The  particles 
were  convected  around  the  wind  tunnel  circuit  and  over  a period  of  a few  min- 
utes filled  the  whole  tunnel.  In  practice,  balancing  the  rate  at  which  parti- 
cles were  introduced  with  the  rate  that  they  were  lost,  in  order  to  keep  the 
particle  number  density  approximately  constant  in  time,  proved  very  difficult 
because  the  characteristic  time  lag  was  of  the  order  of  several  tens  of  min- 
utes. As  a consequence,  for  the  Merrill  tunnel  measurements  the  particle  num- 
ber density  varied  during  the  time  required  for  the  traverse  through  the  bound- 
ary layer. 


Porticla  Oiomatar(^m) 

Fig.  1.  Particle  size  distribution 


Optics  and  measurement  Reometr\ 


The  present  experiments  utilized  the  laser  Doppler  velocimeter  in  the 
single  particle,  dual  (forward)  scatter  mode.  For  measurements  in  the  JPL  20" 
wind  tunnel,  the  instrument  was  mounted  on  a vibration  isolated  traversing 
mechanism  on  top  of  the  wind  tunnel  test  section,  as  depicted  in  figure  2. 

For  measurements  in  the  GALCIT  Merrill  wind  tunnel,  the  traverse  was  suspended 
from  the  ceiling  of  the  laboratory. 

The  light  source  for  this  instrument  was  a Coherent  Radiation  Model  52B, 

4 watt  Argon  ion  laser.  The  laser  was  operated  single  line  at  0.5145  ym  and 
etalon  stabilized  to  provide  a single  mode  beam.  It  was  mounted  on  the  tubular 
support  structure  for  the  optics  as  9hown  in  figure  2. 

The  laser  beam  is  directed  into  the  transmitting  optics  cell  where  it  is 
split  into  three  parallel  beams,  of  approximately  equal  intensity,  which  form 
a right  isosceles  triangle,  whose  base  is  nominally  parallel  to  the  wall.  The 
three  output  beams  are  then  focused  in  the  center  of  the  tunnel  by  a 1 m focal 
length  lens,  to  a common  intersection  volume  ^ 0.8  mm  in  diameter.  The  optical 
axis  intercepts  the  test  plate  from  below  at  an  angle  of  1:100  in  order  to 


permit  measurement  of  the 
boundary  layer  down  to 
the  surface  of  the  pol- 
ished test  plate. 

The  resulting  Inter- 
section volume  contains 
three  independent  sets  of 
(virtual)  fringe  planes 
correspondingly  perpen- 
dicular to  the  u, 

(u  + v)//2  and  (u  - v)//2 
velocity  vectors  (see 
figure  3) . By  selective- 
ly blocking  one  of  the 
three  beams,  any  one  of 
the  three  velocity  compo- 
nents can  be  examined 
without  a change  In  the 
focal  volume  geometry. 

The  optical  axis  for 
the  receiving  optics  Is 
aligned  at  an  angle  of  ^ 

7.5*  with  respect  to  the 
test  plate  surface.  This 
allows  viewing  the  focal  volume  down 
that  the  transmitting  optics  axis  is 
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Fig.  2.  Measurement  assembly 


to  the  test  plate  surface  and,  recalling 
in  turn  inclined  with  respect  to  the  test 
plate,  reducing  the  spanwise  extent  of  the  focal  volume.  The  collected  light 
cone  is  focused  by  a 1:1  imaging  system  into  a 0.05  cm  diameter  pinhole  which 
spatially  filters  the  collected  light  accepted  by  the  photomultiplier  receiving 
optical  assembly.  The  effective  measurement  volume  diameter  is  defined  by  the 
pinhole  in  the  receiving  optics  and  therefore  is  equal  to  0.05  cm. 


Fig.  3.  Measurement  fringes 


5.  Signal  processing  and  data  acquisition 

The  output  of  the  photomultiplier  was  amplified  by  a direct  coupled  am- 
plifier installed  at  the  base  of  the  photomultiplier  housing.  That  output  was 
in  turn  amplified  by  a low-noise,  wide  bandwidth,  AC  coupled  preamplifier. 

The  resulting  signal  was  low-pass  filtered  to  remove  the  substantial  shot 
noise  above  the  Doppler  frequencies  in  each  case.  Note  that  the  pedestal,  re- 
sulting from  the  dual  scatter  optics , is  not  removed  by  this  processing  scheme. 
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The  filtered  signal  Is  fed  Into  a counter  type  processor  (Dlaotakis  and  Lang 
1974,  1977,  and  Dlootakla,  Collins  and  Lang  1978). 

The  processor  uses  a four  level  comparison  scheme  on  every  cycle  of  the 
Doppler  burst.  A TTL  signal  DIN  Is  defined  whose  negative  slope  Is  driven  by 
the  crossing  of  the  reference  level  VQ  from  above  by  the  analog  Input  (see 
figure  4a).  DIN  can  be  prescaled  by  m - 1,2, 4, 8 to  yield  DIN'  which  Is  used 
by  the  timing  circuitry.  DIN'  Is  subsequently  checked  to  ensure  that  every 
period  Is  between  two  limits,  Tj  and  Tj  + T2 . If  T was  the  period  of  the 
Doppler  signal,  we  must  have,  as  a consequence, 

Tj  < mx  < Tj  + t2  . (5) 


The  time  Intervals  and  t2  are  front  panel  selectable  and  cover  the  range 
60  nsec  < Tj,T2  < 3 msec.  The  processor  updates  the  flight  time  register  with 
the  time  between  the  first  nega- 
tive slope  of  DIN'  and  every 
successive  acceptable  negative 
slope  of  DIN'  using  a 100  MHz 
(Tc  ■ 10  ns)  crystal.  Every 
negative  slope  Is  also  recorded 
In  the  cycle  counter. 

The  data  were  recorded  In 
two  different  modes.  In  the 
first  mode  (fixed  fringe  mode), 
the  output  corresponds  to  the 
(truncated)  flight  time  for  a 
fixed  number  of  fringes.  All 
the  data  for  M £ 0.6  (JPL  20" 
wind  tunnel)  were  recorded  in 
this  mode.  Eight  cycles  were 
counted  for  the  beam  pair  at  0s 
(ui  data)  with  the  prescaler  m 
set  to  2.  Four  cycles  were 
counted  for  the  ±45°  (u2,  u3 
data)  beam  pairs  with  the  pre- 
scaler set  to  unity. 

In  the  second  mode  (free 
fringe  mode),  the  processor  only 
requires  that  a fixed  number  of 
cycles  be  reached  or  exceeded  by 
the  buret  but  records  the  flight 
time  and  number  of  cycles 
(fringes)  of  the  whole  burst. 

Both  n and  (the  truncated)  At 
were  recorded  in  this  mode  and 
the  velocity  is 
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Fig.  4.  Signal  burst  processing 


uA  ■ ns/At,  (6) 

where  s Is  the  fringe  spacing.  The  data  recorded  at  M ^ 0.1  (GALCIT  Merrill 
wind  tunnel)  were  acquired  in  this  mode  requiring  at  least  four  fringes  for  the 
145*  data  (u2,  u2)  and  at  least  eight  for  the  0*  data  (uj)  with  the  prescaler 
m set  to  2.  Thus,  In  both  modes,  the  same  (minimum)  number  of  fringes  had  to 
be  crossed. 

There  are  several  Important  differences  between  the  two  modes.  First, 
using  the  facility  of  the  processor  to  recognize  the  end  of  the  burst  while 
retaining  the  previous  valid  measurement,  the  total  flight  time  can  be  used. 
Instead  of  an  arbitrary  fraction.  In  determining  the  velocity  component  of  in- 
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terest.  This  Improves  the  measurement  accuracy  In  an  obvious  way.  Second, 
using  the  whole  burst  cancels  the  phase  errors  that  are  made  in  assigning 
equal  phase  to  the  crossings  of  the  reference  level  VQ  by  the  signal.  Third, 
the  sampling  statistics  are  different  between  the  two  inodes.  In  the  first 
mode  (fixed  number  of  fringes),  all  we  know  Is  that  the  particle  crossed  at 
leaet  that  many.  In  the  second  mode  (number  of  fringes  unrestricted  provided 
it  exceeds  a certain  minimum)  we  know  that  the  particle  crossed  exactly  that 
many.  That  alters  the  sampling  bias,  as  will  be  discussed  later  on. 

The  data  were  formatted  into  records  of  1,024  particle  measurements  each 
and  written  on  a Kennedy  9100  digital  tape  deck.  Several  records  were  re- 
corded at  each  measurement  location.  Data  rates  up  to  50  KHz  were  observed  in 
the  high  speed  flow. 


6.  Data  processing 


In  order  to  process  the  high  speed  data  (0.6  < M 2.2),  which  were  re- 
corded in  the  fixed  fringe  mode,  a histogram  of  each  record  (1,024  measure- 
ments) was  formed  as  a function  of  the  flight  time  At.  The  histograms  were 
pruned  in  two  passes.  First,  the  data  in  any  bin  that  contained  only  one 
count  and  did  not  have  neighboring  bins  with  more  than  one  count  were  discard- 
ed. In  the  second  pass,  any  data  Isolated  from  the  main  body  of  the  histogram 
by  at  least  one  zero  were  also  discarded.  If,  as  a result  of  these  two  oper- 
ations, more  than  24  measurements  out  of  the  total  of  1,024  had  to  be  reject- 
ed, the  whole  record  was  rejected. 

For  the  low  speed  data  (M  - 0.1),  which  were  recorded  in  the  free  fringe 
mode  and  for  which  the  actual  number  of  flight  time  clock  counts  was  large,  an 
alternative  scheme  was  devised.  The  data  were  sorted  in  velocity  bins  as 
follows.  The  velocity  of  each  particle  was  computed  as  an  integer  percentage 
of  the  maximum  velocity  that  the  processor  settings  admitted,  i.e. 


lOOtjnj/Atj, 


(7) 


where  Tl  is  the  minimum  Doppler  period  setting,  n^  is  the  number  of  fringes 
crossed  by  the  i**1  particle  and  Atj  is  the  time  of  flight  corresponding  to  the 
n^  fringes.  The  histogram  was  formed  based  on  the  velocity  index  Ij  and  each 
record  was  pruned  in  the  fashion  described  above . 

The  complexity  of  processing  single  particle  laser  Doppler  velocity  data 
is  compounded  by  the  fact  that  in  such  measurements  the  fluid  velocity  is 
sampled  in  a biased  manner.  This  problem,  first  pointed  out  by  McLaughlin 
and  Tiederman  (1973),  necessitates  that  the  unbiased  expectation  value  of  any 
function  of  the  velocity  u,  be  computed  consistently  with  this  bias,  i.e. 


<f(u)  > 


Ef^)®'1^)  /£s-1(u i> 
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(8) 


where  B(u)  is  proportional  to  the  probability  per  unit  time  of  making  the  par- 
ticular measurement.  The  dependence  of  the  sampling  bias  function  B on  the 
local  velocity  vector  u,  as  well  as  the  shape  of  the  measurement  volume  and 
the  minimum  number  of  fringe  crossings  required  by  the  processing  electronics, 
has  been  derived  elsewhere  (Dimotakis  (976).  If  the  angle  9 between  the  in- 
tersecting beams  is  small  then,  for  the  fixed  fringe  mode  measurements,  the 
expression  for  B( u)  simplifies  to 


sfixed^;E)  “ ui(J  + «5/u»)*ll  - e*(l  + u^/u*)1. 


(9) 


where  u,  is  the  velocity  component  perpendicular  to  the  fringe  planes,  ul(  is 


the  velocity  component  parallel  to  the  fringe  planes  and  perpendicular  to  the 
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long  axis  of  the  measurement  volume  (beam  bisector) , and  e la  the  ratio 
nmin^nT  of  the  minimum  number  of  fringe  crossings  (Doppler  cycles)  required  by 
the  processor  to  the  total  number  (maximum)  In  the  measurement  volume.  If 
(u„ /uA)2  « 1 equation  9 reduces  to  the  familiar  form 

®fixed<H)  « «;*  « (k  + k),  (10) 

where  k Is  the  (Integer)  output  of  the  flight  counter. 

In  the  free  fringe  mode,  the  bias  function  Is  given  by 

®free(— ;n>  “ + u*)3^2(n  + %),  (11) 

or,  equivalently, 

Bfree(k;n)  * (k  + ^)n2/(n  + *>•  (12) 

See  Dimotakis  (1976)  and  Dlmotakla,  Collins  and  Lang  (1978)  for  derivations 
and  a more  extensive  discussion.  Note  that  equation  12  Is  In  disagreement  with 
the  predictions  of  Hoesel  and  Rodl  (1977)  who  suggest  that  Bf *ee « At « (k  + %) . 


7.  Results  and  discussion 

The  mean  streamwlse  velocity  profiles  are  compared  to  the  corresponding 
Pitot  tube  measurements  in  figures  5 and  6.  The  agreement  can  be  seen  to  be 
quite  good.  The  momentum  thickness  6 for  the  profiles  was  computed  from  the 
laser  Doppler  measurements.  In  figure  5 four  profiles  in  the  Mach  number 
range  0.6  4 M 4 2.2  are  plotted  as  Indicated.  The  agreement  for  the  lower  speed 
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profiles  Is  very  good.  We  do  not  have  an  explanation  for  the  small  discrep- 
ancies that  can  be  seen  In  the  higher  Mach  number  profiles.  The  top  data  point6 
In  figure  6 correspond  to  a computation  of  the  mean  velocities  using  the  har- 
monic mean  (bias  function  of  equation  10),  the  bottom  points  correspond  to  the 
more  exact  bias  compensation  of  equation  12.  The  differences  are  small  as  are 
also  the  differences  between  these  and  mean  values  computed  with  no  compensation 
for  the  bias  (not  plotted).  In  figures  7 and  8 the  root  mean  square  fluctu- 
ations, normalized  by  the  edge  velocity,  are  plotted  for  the  same  flow.  The 
upper  set  of  points  In  figure  8 are  computed  using  equation  10  for  the  bias. 
Again  the  differences  between  the  two  are  small.  The  filled  squares  in  figures 
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Fig.  7.  0.6  4 M < 2.2  u’/Ue  profiles  Fig.  8.  M ^ 0.1  u'/Ue  profile 

7 and  8 are  taken  from  Klebanoff  (1955).  The  apparent  fluctuation  level  in 
the  free  stream  for  the  higher  Mach  numbers  is  largely  the  result  of  the 
finite  clock  period  (see  equation  1),  an  effect  that  was  not  removed  from  the 
data.  The  error  bars  are  computed  on  the  basis  of  estimates  from  different 
records.  The  differences  between  the  present  data  and  the  Klebanoff  (1955) 
data  are  within  our  confidence  limits.  The  normalized  Reynolds  stress  profile, 
-Pwu'v'/tw,  measurements  are  presented  in  figures  9 and  10.  The  solid  lines 
represent  the  Reynolds  stress  computed  from  the  corresponding  mean  velocity 
profiles  by  integrating  the  boundary  layer  equations  In  the  manner  described 
in  Collins,  Coles  and  Hicks  (1978).  The  filled  squares  are  taken  from 
Klebanoff  (1955).  The  upper  points  in  figure  10  are  computed  on  the  basis  of 
equation  10  for  the  bias. 

Several  questions  can  be  answered  immediately  by  the  present  data. 

Firet,  the  anomalies  reported  In  previous  supersonic  boundary  layer  measure- 
ments of  the  Reynolds  stress  are  not  a consequence  of  compressibility  effects. 
This  can  be  ascertained  directly  from  figure  9 in  which  no  systematic  varia- 
tion with  Mach  number  of  the  Reynolds  stress  measurements  in  the  vicinity  of 
the  wall  can  be  seen.  Second , the  difficulty  which  does  exist  for  y/8  ■€  2, 
cannot  be  associated  with  local  particle  lag,  in  the  context  of  equation  4, 
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Fig.  9.  0.6  M ■£  2.2  -u'v'  profiles 


Fig.  10.  M 'V  0.1  -u'v’  profile 


since  the  deviation  from  the  expected  behavior  would  have  occurred  at  increas- 
ing distances  from  the  wall  as  the  Mach  number  increased  (the  boundary  layer 
thickness  was  roughly  2.5  cm  for  all  the  measured  profiles).  Third 3 the 
previously  reported  anomalies  cannot  be  attributed  to  the  sampling  bias.  The 
differences  between  the  correctly  computed  Reynolds  stress  data  (lower  data 
points  of  figure  10)  and  the  simpler  compensation  by  means  of  equation  10  are 
small,  even  though  the  bias  function  of  equation  12  yields  results  that  are 
closer  to  the  expected  values. 

Except  for  the  M ^ 0.1  stress  measurements,  the  agreement  between  the 
computed  stress  profiles  and  the  directly  measured  values  is  quite  good,  for 
y/0  Z 2.  The  discrepancy  near  the  wall,  however,  remains  an  open  question  as 
of  this  writing.  The  problem,  we  feel,  is  related  to  the  behavior  of  the 
stress  producing  flow  in  the  neighborhood  of  the  wall.  It  is  a common  obser- 
vation that  the  walls  of  the  tunnel  become  coated  with  a thin  film  of  the 
aerosol  used  for  seeding.  This  Indicates  that  the  flow  inside  the  viscous 
sublayer  is  such  that  a particle  that  enters  the  sublayer  has  a very  low  prob- 
ability of  leaving.  The  scaling  law  for  particle  behavior  in  a turbulent 
boundary  layer,  as  expressed  by  equation  4,  is  not  applicable  in  the  viscous 
sublayer.  A particle  entering  the  viscous  sublayer  becomes  caught  in  the  mo- 
tion of  the  longitudinal  vortices,  whose  transverse  extent  is  approximately 
20  v/uT,  and  ends  up  on  the  wall.  Consequently,  a fluid  element  coming  from 
the  wall  is  less  likely  to  carry  particles  than  a fluid  element  moving  touarda 
the  wall.  This  causes  several  problems  that  show  up  In  different  measurements. 
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The  scattering  particle  number  density  profile  for  the  JPL  data  (0.6  4 M 
< 2.2)  should  look  approximately  like  a half  Gaussian,  resembling  the  concen- 
tration of  a passive  contaminant  Introduced  far  upstream  at  the  wall.  The  re- 
corded data  rate  was  divided  by  the  local  streaawlse  velocity  and  plotted  in 
figure  11.  To  the  extent  that  correlations  between  the  scattering  particle 
number  density  fluctuations 
and  the  etreesnsise  velocity  ' I 

can  be  neglected,  the  ratio  Is  • 

proportional  to  the  local  par-  |«*  *• 

tide  matber  density.  The  de-  _ • * 

crease  as  the  wall  Is  ap-  ( • * 

proached,  extending  to  a large  Kl  ^ ^••****  * * • • 

fraction  of  the  boundary  layer  ,•*■*  * * 

thickness.  Is  quite  consplcu-  y***  • 

ous.  The  velocity  component  >-  • 

perpendicular  to  the  wall  Is  — ^ _ • 

msasured  to  be  negative  in  the  * ^ * • * 

vicinity  of  the  wall,  a conse-  Q j>V  * • _ •_ 

quence  of  the  fact  that  the  £ •*  *•  * m 

positive  half  of  the  probabll-  f • • 

Ity  density  of  that  quantity  i ■.•  ■»«,T[  • * . 

Is  underrepresented  In  the  " *• 

laser  Doppler  measurements.  y <*•  a *.  w.i.w 

The  underrepresentation  of  £ * • 

that  motion  also  results  in  a 

lower  measured  value  for  the  • , •< 

-u'v'  correlation  since  the  ••  * 

upswelling  carries  fluid 

with  lever  streamwlse  veloc-  t>  ®e 

lty.  This  motion  (bursting),  se 

which  has  been  observed  to  0 * ® 

be  intermittent  and  quite  01 

violent,  is  held  responsible  °r  r it r .rT~ 

for  a very  large  fraction  of  y/thetr 

the  total  stress  near  the 

wall  (Blackwelder  and  Kaplan  Fig.  11.  Particle  number  density  profiles 

1976). 

It  appears  that  the  particle  transport  to  the  wall  via  the  viscous  layer 
is  so  effective,  even  at  the  lower  velocities,  as  to  make  the  wall  look  like 
an  infinite  sink  for  particles.  The  conclusion  is  then  that  the  discrepancy 
in  the  measured  Reynolds  stress  near  the  wall  arises  because  the  scattering 
particle  nwiber  density  is  very  strongly  correlated  with  the  instantaneous 
stress  and  not  as  a result  of  any  local  failure  of  the  particles  to  track  the 
flow. 
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Summary 


Experiments  have  been  carried  out  using  laser-Doppler  velocimetry  in 
conjunction  with  Pitot  and  surface-friction  instrumentation  for  the  measure- 
ment of  the  mean  flow,  the  velocity  fluctuations,  and  the  Reynolds  stresses 
in  turbulent  boundary  layers  at  constant  pressure  in  subsonic  and  supersonic 
flows.  TWo  facilities  were  used  for  these  studies.  In  one  facility,  data 
were  obtained  on  a flat  plate  at  M£  = 0.1,  with  Reg  = 6600.  In  the  other, 
data  were  obtained  on  an  adiabatic  nozzle  wall  at  *»  0.6,  0.8,  1.0,  1.3, 

2.2  with  Reg  - 23,000  and  40,000. 

Part  I of  this  report  is  a detailed  analysis  of  the  mean  flow  as  observed 
using  Pitot-tube,  Preston-tube,  and  floating-element  balance  instrumentation. 
The  emphasis  in  Part  I is  on  the  use  of  similarity  laws  with  Van  Driest  scaling 
and  on  the  inference  of  the  distribution  of  shearing  stress  and  of  the  normal 
velocity  component  from  the  equations  of  mean  motion.  These  data  serve  as  a 
background  for  the  analysis  of  the  laser-Doppler  data  reported  herein. 

The  present  Part  II  of  this  report  is  a description  of  both  the  mean  and 
the  fluctuating  components  of  the  flow,  and  of  the  Reynolds  stress  as  observed 
using  a dual  forward-scattering  laser-Doppler  velocimeter.  A detailed  descrip- 
tion of  the  instrument,  and  of  the  data  analysis  techniques,  have  been  included 
in  order  to  fully  document  the  data.  A detailed  comparison  is  made  between 
the  laser-Doppler  results  and  those  presented  in  Part  I,  and  an  assessment  is 
made  of  the  ability  of  the  laser-Doppler  velocimeter  to  measure  the  details 


of  the  flows  involved. 
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I*  Introduction 
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The  turbulent  boundary  layer  at  constant  pressure  has  been  the  subject 
of  experimental  and  theoretical  investigations  for  many  years,  and  provides 
a well  documented  flow  for  the  assessment  of  experimental  techniques.  These 
investigations  have  shown  that  for  the  incompressible  turbulent  boundary  layer, 
the  turbulent  shearing  stress  can  be  measured  directly  or  can  be  calculated 
from  the  distribution  of  mean  velocity  with  the  aid  of  well-established  simi- 
larity laws.  For  compressible  flows,  however,  measurements  of  the  Reynolds 
stresses  are  rare,  and  interpretation  of  the  results  is  difficult.  Recent 
measurements  by  Johnson  and  Rose  (1973),  Yanta  and  Lee  (1974),  and  by  Abbiss 
(1976)  have  used  laser-Doppler  velocimetry  techniques  to  make  direct  measure- 
ments of  the  Reynolds  stress  in  turbulent  boundary  layers  with  free-stream 
Mach  numbers  in  the  range  1.5  to  3.0.  However,  a serious  anomaly  is  exhibited 
by  these  measurements  in  that  the  maximum  value  of  -pu  v iccurs  much  further 
from  the  wall  than  is  reasonable  for  flew  at  constant  pressure.  This  anomaly 
has  been  discussed  by  Sandborn  (1974),  who  supports  the  conjecture  by  some  of 
the  authors  cited  that  density  fluctuations  may  contribute  substantially  to 
the  turbulent  stresses  near  the  wall.  This  conjecture  is  in  direct  opposition 
to  the  conclusion  by  Morkovin  (1961)  that  effects  of  density  fluctuations 
should  be  small  compared  to  effects  of  variations  in  mean  density  for  Mach 
numbers  up  to  4 or  5. 

The  purpose  of  the  present  experiments  is  to  obtain  redundant  data  over 
a substantial  range  of  Mach  numbers  (Mg  ■ 0.1  to  2.2),  in  an  effort  to  resolve 
the  anomaly  in  turbulent  shearing  stress.  The  low-speed  experiments  were 

•ff 

performed  in  the  boundary  layer  on  a flat  plate  model  in  the  GALCIT  Merrill 

Graduate  Aeronautical  Laboratories  of  the  California  Institute  of  Technology. 
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wind  tunnel  at  a free-stream  Mach  number  0.1.  The  high-speed  experi- 
ments were  performed  in  the  ceiling  boundary  layer  of  the  J.P.L.*  20-inch 
wind  tunnel  at  Mach  numbers  0.6  S Me^  2.2. 

The  flow  was  documented  by  conventional  means  using  a Pitot  tube, 
which  was  traversed  through  the  boundary  layer,  to  measure  the  mean  flow. 

In  addition,  surface-friction  measurements  were  made  using  both  a floating- 
element  balance  and  a Preston  tube.  The  mean-flow  scaling  suggested  by 
Van  Driest  (1955)  was  applied  to  the  data,  to  test  the  adequacy  of  a single 
similarity  formulation  for  both  compressible  and  incompressible  flow.  The 
shearing  stress  was  computed  from  the  mean  flow  as  part  of  the  analysis. 

The  details  of  these  experiments  have  been  discussed  in  Part  I of  this 
report  by  Collins,  Coles  and  Hicks  (1978). 

In  the  present  report,  the  instrumentation,  data  acquisition  and 

analysis  for  laser-Doppler  velocity  measurements  in  the  same  flows  are 

/2  ~t2  — 7 — 7 

I described.  Measurements  of  u,  v,  u , v and  the  Reynolds  stress  - pu  v 

were  made  over  the  full  Mach  number  range.  These  measurements  are  discussed 
and  compared  to  the  results  for  the  mean  flow  presented  in  Part  1 of  this 
work. 

II.  Measurement  Considerations  at  High  Speeds 

The  use  of  laser-Doppler  velocimetry  techniques  in  high-speed  flows 
unavoidably  results  in  very  high  frequency  signals.  Figure  1 is  a plot  of 
the  Doppler  frequency 

v_  = 2 sin8/2  u , (1) 

D \ x 
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as  a function  of  the  flow  velocity  u with  the  beam  separation  angle  6 


as  a parameter.  For  these  computations,  the  wavelength  has  been  assumed 


to  be  \ m 0^5145  pm,  corresponding  to  an  argon- ion  laser 


reduction  in  the  signal  frequency,  results  in  a reduction  in  the  number 


of  fringes  within  the  measurement  ellipsoid,  and  a consequent  reduction 


virtual  fringes  whose  spacing  is  given  by 


Because  the  (TEM^  Gaussian  beam  radius  at  the  focus  is  given  by 


where  R is  the  beam  wavefront  radius  of  curvature  after  the  focusing  lens 


(R  ~ f = lens  focal  length),  and  w is  the  1/e  intensity  envelope  radius 


at  the  focusing  lens,  the  resulting  number  of  fringes  is  then  given  by 


(e.g.,  Dimotakis  (1976)) 


tan6/2 


The  number  of  fringes,  Nf,  is  plotted  in  Figure  2 as  a function  of  the 


focal  volume  diameter  2 w with  the  angle  6 as  a parameter  for  a wavelength 


X * 0.5145  pm*  As  can  be  seen  from  Figure  2,  spatial  resolution  consid 


erations  place  a lower  bound  on  the  angle  6,  presenting  a criterion  which 


conflicts  with  the  requirement  to  keep  the  resulting  Doppler  frequencies 
small. 


Other  difficulties  arise  if  the  Doppler  frequencies  are  permitted 
to  be  large.  In  addition  to  the  problems  of  handling  high  frequencies 
per  se,  the  relative  accuracy  6u/u  with  which  the  velocity  of  a single 
particle  can  be  measured  is  never  less  than  the  ratio  of  the  processor 
clock  period,  T^,  to  the  total  flight  time  At  that  is  used  for  the  mea- 
surement, i.e.. 


6u  > 
u At 


(5) 


By  way  of  example,  a processor  with  a clock  frequency  of  100  MHz 
-8 

(Tc  ■ 10  sec)  can  measure  the  velocity  of  a Mach  ~ 2 particle  in  air 
(u^  ~ 530  m/sec)  in  a 800  pm  diameter  focal  volume  to  about  ~ 0.7%. 

In  addition  to  these  considerations,  which  would  limit  the  accuracy 
even  if  the  measurement  environment  was  noise  free,  other  factors  become 
important  at  high  velocities  by  decreasing  the  signal-to-noise  ratio  of 
a single  reading.  As  the  flow  velocities  become  higher,  smaller  particles 
are  required  to  minimize  problems  of  particle  lag.  This  results  in  rapidly 
decreasing  scattering  cross-sections  (Born  and  Wolf  (1975))  and  consequent 
signal  intensities.  In  addition,  as  the  velocity  increases,  the  scattering 
particle  spends  less  time  in  the  focal  volume.  As  a result,  even  though 
the  number  of  scattered  photons  per  second  remains  constant,  proportional 
to  the  scattering  cross-section,  fewer  total  photons  are  collected  by  the 
detecting  optics.  An  equivalent  way  of  analyzing  this  problem  is  to  consider 
the  noise  bandwidth  that  must  be  utilized  as  the  velocities,  and  resulting 


i ; 


■ 


4 


I 


8 


Doppler  frequencies,  increase.  For  a shot-noise  limited  detection  process, 
the  noise  power  increases  proportionally  to  the  bandwidth  that  is  accepted 
by  the  processing  electronics.  Consequently,  using  the  same  particles, 
the  signal-to-noise  ratio  will  decrease  as  the  velocity  increases. 


III.  Scattering  Particles 
A.  Particle  Lag 

The  extent  to  which  the  measurement  of  particle  velocity  represents 
the  fluid  velocity  is,  of  course,  a separate  issue.  If  the  Lagrangian 
fluctuation  frequencies  in  the  fluid  are  small  compared  to  the  reciprocal 


of  a characteristic  particle  response  time,  t , the  particles  are  expected 


to  follow  those  fluctuations. 

For  small  solid  particles  in  gas  flows,  the  equation  of  motion  in 
a Lagrangian  frame  can  be  approximated  by 


d .1  1 

— u + — u — u„ 
dt  p T p T f 

p p 


(6) 


If  the  particle  Reynolds  number 


Re  e . 

P v P 


(7) 


is  small,  and  the  flow  can  be  described  by  the  continuum  equations,  we 
have  the  classical  Stokes  flow  for  a spherical  particle  for  which  the 


drag  coefficient  is  Cp  “ 24/Re^.  The  particle  response  time  is  then  given 


by 


, . i»2(i) 

p 18  ' (i  / 


2 

d z /p 


(- E.) 

18v  ' pf / * 


(8) 
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which  is  independent  of  the  fluid  density.  A plot  of  the  particle  response 
time  as  a function  of  the  particle  diameter  is  given  in  Figure  3 for  di- 
butyl phthalate  aerosol  particles  in  air. 

In  order  to  estimate  the  fluctuation  frequencies  to  which  a particle 
will  be  subjected,  it  is  useful  to  examine  the  Eulerian  fluid  motions  in 
the  turbulent  boundary  layer.  From  the  measurements  of  Klebanoff  (1954) 
in  a boundary  layer  at  constant  pressure,  and  from  the  measurements  of  Perry 
and  Abell  (1975;  1977)  in  pipes,  the  behavior  of  the  Eulerian  velocity 
fluctuation  spectrum  can  be  inferred.  These  results  give  spectra  which  can 
be  described  by 


Suu(<t,»y)  TT6  f f (u6 

u(y)  °Lu(y). 

u (y) 


for  the  wake  region,  and 


S (u>,y) 
uu 

u,2(y) 


u(y) 


(9) 


(10) 


in  the  inner  portion  of  the  boundary  layer,  but  outside  the  sublayer. 
These  spectra  then  result  in  expected  frequencies  given  by 


(u))g  ■ const. 


(ID 


in  the  wake  region,  and 

(to)  ■ const.  u(y). 

y y 


(12) 
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in  the  boundary  layer  interior.  Consequently,  particles  above  the  sublayer 
should  track  the  flow  well,  provided  that 


y(y)  <1  , 


< J_  r y/6  1 

~ Me  Lu(y)/ul 


Substituting  Equation  (8)  for  t , we  obtain  an  upper  bound  on  the  particle 
diameter  for  the  particles  to  track  the  flow, 


dP  < dPmax  ^ . 

Me 


for  di-butyl  phthalate  aerosol  in  air  at  25  C. 

A prediction  of  Equation  (12)  is  that  the  turbulent  wavenumbers  k ~ oj/u 
increase  as  the  wall  is  approached. 


const . 


This  result  is  confirmed  qualitatively  by  the  photograph  in  Figure  4 in 
which  the  turbulent  boundary  layer  on  the  surface  of  a water  channel  has 
been  made  visible  using  a thick  suspension  of  aluminum  flakes  which  pref- 
erentially align  with  the  local  strain  field.  The  photograph  was  taken  in 
the  GALCIT  low-speed  water  channel  by  Brian  Cantwell  and  is  reproduced  here 
with  his  permission.  The  Reynolds  number  for  this  photograph,  based  on  x. 


HflMMI 


was  10  . The  monotonic  decrease  in  scale  size  as  the  wall  is  approached 
is  quite  evident.  These  observations,  and  the  results  of  Equation  (12)  are 
consistent  with  the  usual  mixing-length  hypotheses  made  for  the  outer  layer 
by  Prandtl  (1942),  and  are  consistent  with  the  hypothesis  for  the  region 
closer  to  the  wall,  but  outside  the  viscous  layer  by  von  Kannan  (1930). 

In  the  interior  of  the  boundary  layer,  the  quantity  in  the  brackets 
in  Equation  (14)  is  bounded  by  the  relation 


iLL. 

u(y)/ue 


u /u 
e t 

6u  /v 
T w 


(17) 


where  u is  the  friction  velocity  defined  as 
T 

\ s ■ <18> 

The  ratio  ue+/6+  has  been  shown  by  Coles  (1968)  to  be  a function  of 

Reg  for  incompressible  flows  with  zero  pressure  gradient.  The  square  root 

of  this  quantity  is  plotted  in  Figure  5,  for  the  range  of  Reg  covered  by 

3 4 

the  present  measurements  (8  x 10  < Reg  < 4 x 10  ) . The  data  for  this  figure 

are  taken  from  Table  4 of  Collins,  Coles  and  Hicks  (1978).  These  data  are 
compared  in  Figure  5 to  a straight  line  fit  given  by  the  emperical  relation 

yu>+  - f (Re0)  = 5.01  Reg'0,435  . (19) 

In  writing  the  lower  bound  in  Equation  (17),  effects  of  compressibility 
have  been  ignored.  Using  the  numerical  result  of  Equations  (15)  and  (19),  we 
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then  have 


r 


d 

p 


< d 


max 


4^3  <*,(, 
Me^  6 (pm) 


24.2 

"eX 


0.435 


far  from  the  sublayer 


near  but  outside  the  sublayer  . 


(20) 


By  way  of  example,  a flow  of  ~ 1 requires  particles  less  than  1.5  [if n 
in  diameter  to  track  the  fluctuations  at  y/6  ~ 0.1. 

B.  Techniques  for  Seeding  the  Flow 

The  particle  generator  for  the  present  experiments  used  a Laskin 
nozzle  type  construction  to  generate  a poly-dispersed  aerosol  of  di-butyl 
phthalate.  Particles  greater  than  1 pm  in  diameter  were  effectively 
removed  by  an  impact  plate  incorporated  in  the  design.  The  resulting 
particle  size  distribution,  measured  using  both  a cascade  impactor  and 
a multi-channel  particle  analyzer,  is  shown  in  Figure  6.  These  results 
indicate  little  contribution  to  the  number  density  for  sizes  greater  than 
1 pm.  While  little  is  known  about  the  distribution  below  0.4  pm,  it  was 
possible  to  ascertain  that  no  particles  with  a diameter  less  than  0.1  pm 
were  present. 

IVo  different  methods  were  used  to  introduce  particles  into  the  flow 
depending  on  the  flow  facility  that  was  used.  For  the  high-speed  measurements 
(0.6  ~ Mg~  2.2),  which  were  conducted  in  the  JPL  20- inch  wind  tunnel,  it  was 
necessary  to  seed  the  flow  by  introducing  the  particles  into  the  settling 
chamber  between  the  last  turbulence  screen  and  the  contraction  section. 

The  particles  were  introduced  through  a tube  which  protruded  5 cm  vertically 
into  the  settling  chamber  from  above.  Holes  were  drilled  into  the  tube  along 
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the  stagnation  line  ns  well  ns  the  rear  side  nt  + 30  with  lospect  to  the 


fl*iw  vector.  Using  this  technique,  data  rates  in  excess  of  SO, 000  samples 
per  second  were  achieved.  Preston-tuhe  measurements  obtained  w i t h and 
without  pare  ides  indicate  no  detectable  influence  of  the  Introduction  of 
particles  on  the  surface  friction  and  therefore,  probably,  on  the  mean  flow. 

For  the  low-speed  measurements  in  the  Merrill  wind  tunnel  (f\.  • 0.1), 
although  some  naturally  occurring  particles  were  present,  the  tluw  was  seeded 
using  the  same  aerosol  generator  in  order  to  increase  the.  data  tate  and  to 
control  the  scattering  particle  sire  distribution.  For  these  measurements 
the  particles  were  injected  through  a tube  spanning  the  test  section  located 
dow  list  ream  of  the  boundary- layet  plate,  providing  a uniform  distribution  ot 
particles  throughout  the  test  section.  In  practice,  the  adjustment  of  the 
rate  at  which  particles  were  introduced  proved  difficult  because  the  charac- 
teristic lime  was  ot  the  order  ot  several  minutes  for  equilibrium  with 
particle  losses.  As  a consequence,  tlx*  part Icln  number  density  changed 
during  the  tins'  required  lor  a traverse  through  the  boundary  layer.  Using 
this  seeding  technique,  data  rates  iu  excess  ot  2,000  samples/sec  we  re 
achieved  at  low  speeds. 

C.  The  Dist  ribut ion  ot  Particles 

Tlie  particle  number  density  was  measured  by  counting  the  average  fre- 
quency of  occurrence  of  valid  data  from  the  processing  electronics,  lltc so 
data  are  a measure  ol  the  mean  particle  t lux  (mi) . Assuming  that  correlations 
between  the  scattering  particle  number  density  and  the  stioaiuwtse  velocity 
are  small,  then 


( 
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gives  the  mean  value  of  the  particle  number  density. 

It  is  an  unfortunate  consequence  of  this  measurement  technique  that 
the  measured  value  of  the  particle  flux  (nu)  decreases  with  the  signal-to- 
noise  ratio  of  the  experiment.  The  most  common  cause  of  a decrease  in  the 
signal-to-noise  ratio  in  the  high-speed  flow  experiments  has  been  attributed 
to  the  random  occurrence  of  oil  filaments  on  the  windows  through  which  the 
observations  were  made.  Imperfections  in  the  lucite  windows  in  the  Merrill 
wind  tunnel  had  a similar  effect  on  the  measured  particle  flux.  In  addition, 
a more  serious  problem  arose  in  the  low-speed  measurements,  where  the  particle 
flux  was  a slowly-varying  function  of  time. 

The  resulting  distribution  of  particle  number  density  through  the 
boundary  layer  is  shewn  in  Figures  7,  8 and  9,  in  which  the  particle  number 
density  has  been  normalized  by  the  local  fluid  density  in  order  to  account 
for  the  effects  of  compressibility.  These  data  are  further  normalized  by 
their  maximum  value  in  the  boundary  layer  to  eliminate  the  effects  of  differing 
injection  rates. 

The  loss  mechanisms  of  the  scattering  particles  from  the  boundary  layer 
are  important  in  understanding  the  resulting  laser-Doppler  observations.  For 
the  high-speed  experiments,  the  particle  distribution  was  initially  both 
spatially  and  temporally  approximately  viniform  over  a restricted  volume  of 
fluid  in  the  settling  chamber.  The  particle  distribution  in  the  outer  portion 
of  the  boundary  layer  is  then  determined  by  the  diffusion  of  the  passive 
scalar  into  the  outer  flow.  In  addition  to  diffusion,  the  entrainment  of 
fluid  into  the  boundary  layer  from  the  free  stream  by  the  mechanism  of  large- 
scale  eddy  engulfment  contributes  substantially  to  the  structure  of  the  particle 
number  density  profiles.  Because  the  entrained  fluid  contains  no  particles. 
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the  measured  distributions  shown  in  Figures  7 and  8 exhibit  a rapid  decay 
away  from  the  wall,  with  the  particle  number  density  falling  to  zero  in 
the  outer  portion  of  the  boundary  layer.  This  entrained  fluid  is  under- 
represented in  the  Doppler  measurements  because  of  the  loss  of  statistical 
we ight . 

As  the  wall  is  approached,  the  particle  number  density,  normalized 
by  the  local  fluid  density,  decays  more  rapidly  than  the  velocity  as  a result 
of  the  loss  of  particles  to  the  wall.  This  loss  to  the  wall  is  evident  from 
the  observation  of  an  oil  film  on  the  test  surface. 

IV . Optics  and  Measurement  Geometry 
A.  Position  Measurement 

The  present  experiments  utilized  the  laser-Doppler  velociroeter  in  the 
single  particle,  dual  forward-scatter  mode.  The  instrument  was  mounted  on 
a vibration  isolated,  two-axis  traverse  which  permitted  surveys  up  to  7S  cm 
in  the  streamwise  direction,  and  20  cm  in  the  direction  normal  to  the  test 
plate.  For  measurements  in  the  JPL  20-inch  wind  tunnel,  the  instrument  was 
mounted  on  top  of  the  wind  tunnel  test  section,  as  depicted  in  Figure  10. 

For  measurements  in  the  GALCIT  Merrill  wind  tunnel,  the  traverse  was  suspended 
from  the  ceiling  of  the  laboratory,  and  was  not  fixed  relative  to  the  test 


plate. 

The  vertical  position  was  sensed  by  means  of  a helipot  mounted  on  the 
vertical  traverse  and  was  repeatable  to  better  than  0.005  cm.  The  origin  for 
the  vertical  position  measurements  was  estimated  by  placing  a razor  blade  of 
known  width  normal  to  the  test  surface  at  the  location  of  the  focal  volume. 

The  origin  was  then  measured  by  determining  the  half  intensity  of  transmission 
as  a function  of  position  of  one  of  the  beams. 
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In  practice,  establishing  the  origin,  to  within  the  accuracy  warranted 


by  the  data,  by  direct  measurement  proved  difficult  for  a variety  of  reasons 


These  inherent  difficulties  in  measuring  the  position  of  the  wall,  relative 


to  the  LDV  traversing  mechanism  necessitated  establishing  the  origin  in  y 


using  the  data  for  the  mean-velocity  profile.  This  was  accomplished  by  a 


least-squares  fit  of  the  mean-velocity  data  to  the  law  of  the  wall  and  the 


law  of  the  wake  in  the  form 


and  where  y is  the  small  offset  of  the  position  of  the  wall  required  to 


establish  the  proper  origin 


For  this  purpose,  a least-squares  analysis  was  performed  to  determine 


the  parameters  u , II  and  y with  k = 0.41  and  c = 5.0  given,  and  6 was 


computed  as  a constraint  imposed  by  the  local  friction  law 


The  technique  for  this  analysis  is  outlined  in  Part  I of  this  work  by 


Collins,  Coles  and  Hicks  (1978)  in  conjunction  with  the  mean-flow  data 


presented  there.  The  results  of  this  procedure  lead  to  corrections  yQ 
of  the  order  of  the  focal  volume  diameter.  In  particular,  the  maximum 
correction  was  0.18  cm  for  one  case  of  the  high-speed  data  and  0.075  cm 
for  the  low-speed  flow  measurements. 

B.  Transmitting  Optics 

The  light  source  for  this  instrument  was  a Coherent  Radiation  Model 
52B,  4-watt  argon- ion  laser.  The  laser  was  operated  single  line  at 
0.5145  nm  and  used  an  oven  stabilized  etalon  to  provide  a single-mode  beam. 

It  was  mounted  on  the  tubular  support  structure  for  the  optics  as  shown  in 
Figure  10.  This  configuration  ensures  that  the  focal  volume  remains  unchanged 
as  the  assembly  is  traversed  and  also  that  the  relative  alignment  between  the 
transmitting  and  receiving  optics  is  preserved  to  within  a few  percent  of  the 
focal-volume  diameter. 

The  laser  output  beam  was  directed  into  the  vertical  support  tube  for 
the  transmitting  optics  by  means  of  two  separate  mirrors  in  order  to  preserve 
the  proper  polarization  of  the  beam.  The  transmitting  optics  cell  contains 
a fixed  mirror  to  direct  the  beam  into  the  horizontal  plane  maintaining  a 
vertical  polarization  vector.  The  beam  is  then  split  into  two  beams  of 
unequal  intensity  by  reflection  from  the  front  and  rear  surfaces  of  the  first 
of  two  multiply-coated  beam  splitters.  The  resulting  beam  pair  is  then  trans- 
mitted through  a second  multiply-coated  beam  splitter  which  divides  the  stronger 
beam  into  two  beams  of  equal  intensity.  This  arrangement  results  in  three 
parallel  beams  of  approximately  equal  intensity  which  form  a right  isosceles 
triangle  whose  base  is  nominally  parallel  to  the  test  surface.  The  three 
output  beams  are  then  focused  in  the  center  of  the  tunnel  by  a 1.0  meter  focal 
length,  antireflection  coated  lens,  to  a common  focal  volume  ~ 0.8  mm  in 


diameter.  The  optical  axis  intercepts  the  test  plate  from  below  at  an 
angle  of  ~ 1:100  in  order  to  permit  measurement  of  the  boundary  layer 
close  to  the  surface  of  the  polished  test  plate. 


I 


The  resulting  focal  volume  contains  three  independent  sets  of  virtual 
fringe  planes  correspondingly  perpendicular  to  the  u,  (u  + v) / J 2 and 
(u  - v)/y  2 velocity  vectors.  The  geometry  of  the  three  beam  pairs,  and 
actual  photographs  of  the  fringes  formed  by  each  beam  pair,  are  shown  in 
Figure  11.  By  selectively  blocking  one  of  the  three  beams,  any  one  of  the 
three  velocity  components  can  be  examined  without  a change  in  the  focal- 
volume  geometry.  This  method  we  consider  superior  to  the  more  common  one 
of  rotating  a single  beam  pair  in  the  desired  orientation  which,  in  our 
experience,  cannot  easily  be  done  with  sufficient  reproducibility  and 
without  slight  displacements  of  the  focal  volume. 

C.  Geometry  of  the  Intersection  Volume 

The  three  intersecting  beams,  taken  in  pairs,  define  three  planes, 
along  which  the  measurements  were  taken.  Even  though  considerable  care 
was  taken  to  orient  these  planes  at  0°  and  ± 45°  with  respect  to  the  wall, 
the  small  deviations  that  were  measured  were  taken  into  account  in  the  data 
analysis . 

The  subtended  angles  between  the  beam  pairs  were  measured  by  two 
different  methods  on  two  different  occasions  and  found  to  be  within  ~ 0.37. 
of  the  values: 
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Thus,  e2  °*  e3  define  a nearly  isosceles  triangle  as  illustrated 

in  Figure  11. 

As  a consequence  of  the  requirement  to  isolate  the  traversing  mechanism 
from  vibrations,  the  relative  orientation  between  the  hypotenuse  of  the 
triangle  and  the  wall  (angle  0,  Figure  12)  had  to  be  determined  by  a separate 
measurement  for  each  profile.  This  was  accomplished  by  permanently  mounting 
a second  razor  blade  to  one  side  of  the  test  plate  with  its  edge  aligned  with 
the  surface  of  the  plate.  The  relative  orientation  for  the  optics  could  then 
be  determined  by  measuring  the  offset  between  the  two  beams  forming  the  triangle 
hypotenuse.  This  angle  0 is  small  in  every  case  but  must  be  accurately  known 
to  estimate  any  quantities  that  involve  v,  the  velocity  component  normal  to 
the  wall. 

Using  the  geometry  as  defined  in  Figure  12,  we  have  from  the  law  of 


cosines. 


and 


cos  a 


cos  3 


cos  y 


(28) 


Assuming  now  that  the  velocity  vector  in  two  dimensions  can  be  defined 
by  its  components  in  the  three  directions  defined  in  Figure  12,  we  have 
that  the  three  measured  components  are  given  by 
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where 


b^  = cos0 

b21  * cos(y  + 0) 
b31  - cos (3  - 0) 
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sin0 

sin(v  + 0) 

-sin(3  ~ 0) 


(30) 


Equations  (29)  form  a redundant  set,  and  can  be  used  to  check  the  validity 
of  the  data.  In  particular,  since  sin(3  + v)  “ sino,  we  must  have 


(uj)  sincr  = (u2>  sin3  + (u^)  siny  . 
Recall  from  Equation  (1)  that 

2 sin0./2  6. 

V = uia'Tui  * 

so  that  the  consistency  Equation  (31)  becomes 


(31) 


(32) 


> 


<v + <v  • 


(33) 


Alternatively,  Equations  (29)  can  be  solved  in  pairs  to  obtain  two  in- 
dependent estimates  of  (u)  and  (v)  at  each  location. 
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Inverting  Equations  (34)  and  (35),  we  then  have 


(iO  sin(y  + 0)  - <u2)  sin0 
siny 


<u2)  cos (0)  - (ux>  cos(y  + 0) 


<uU> 


(u1>  Sin(p  - 0)  + <u3>  sin0 


<VII> 


(u3>  COS (P  - 0)  - <U3>  COS0 


The  final  values  for  (u)  and  (v)  were  then  estimated  by  averaging  the  two 
independent  estimates. 


<U>  - ~ (<ux>  + (un>) 

, X - <- 

<v>  = 2 (<VI>  + <VII>) 

* 

Using  Equations  (29),  we  can  also  relate  the  mean-square  fluctuations 


<VI>  + <VII> 


/2  il  tt 

of  the  measured  velocity  components,  to  (u ' ),  (v  ) and  (u  v ).  In 
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particular,  we  have 


/ / 2 2//2»  , 2 

<u  1 > “ bll  <u  > + b12 

z / 2v  . 2 - /2X  2 

<U  2 > " b21  <u  > + b22 
<u  3 > " b31  <U  > + b32 
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Note  that  while  Equations  (29)  are  redundant,  Equations  (39)  are  not, 
except  for  special  choices  of  p,  y and  0.  Inverting  Equations  (39), 
we  then  have 

<u'2>  « Cu  (u^2)  + C12  (u'22>  + (u^2) 

<v/2>  - C21  (u^2)  + C22  <u  '22>  + C23  <u'32> 

(u'v')  - C31  <u  /12>  + C32  (u'22)  + C33  (u'32>  J , 


(40) 


where  the  matrix  is  given  by. 


^ *\ 


Cjj  = sin(P  + v)  sin(p  - 0)  sin(y  + 0)/A 
C^2  = -sinP  sin(P  - 0)  sin0/A 


C^3  **  siny  sin(y  + 0)  sin0/A 


C21  = -sin(3  + y)  cos(p  - 0)  cos(y  + 0) /A 
C22  = sing  cos(p  - 0)  cos0/A 


C23  **  siny  cos(y  + 0)  cos0/A 


V > 


(41) 
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sin(P  + y)  sinO  - Y “ 20)/A  ^ 


C32  = ^ sinp  sin(P  - 20)/A 
C33  “ "2  sinY  sin^Y  + 20)/A 
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and  where 


A “ sing  • siny  • sin(&  + y)  • 


(42) 


D.  Receiving  Optics 

The  optical  axis  for  the  receiving  optics  is  aligned  at  an  angle  of 
~ 7.5°  with  respect  to  the  test-plate  surface,  thus  permitting  a full  view 
of  the  focal  volume  throughout  the  boundary  layer.  This  inclination  of 
both  the  transmitting  and  the  receiving  optics  with  respect  to  the  test 
surface  results  in  a reduction  of  the  spanwise  extent  of  the  focal  volume. 

The  collecting  lenses  consist  of  a pair  of  15.24  cm  diameter,  f/5  achromats 
which  have  been  anti-reflection  coated  and  are  mounted  back-to-back  in  a 
common  cell.  This  configuration  provides  a 1:1  imaging  system  that  was 
measured  to  be  approximately  twice  diffraction  limited.  The  collected 
light  cone  is  directed  by  a 10  cm  diameter  multi-layer  dielectric  coated 
mirror  into  a 0.05  cm  diameter  pinhole  which  spatially  filters  the  collected 
light  that  is  accepted  by  the  photomultiplier  assembly.  The  photomultiplier 
assembly  consists  of  a focusing  mount  with  a 20  mm  focal  length,  f/5  coated 
lens,  a 30A  wide  0.5145  pm  filter  and  magnetic  lensing  to  restrict  the  photo- 
cathode of  an  EMI  9813B  photomultiplier  tube.  The  photomultiplier  tube  output 
is  coupled  directly  to  a low-noise  fast  rise-time  preamplifier,  constructed 
inside  the  tube  housing  at  the  base  of  the  photomultiplier  tube.  The  photo- 
detection system  was  designed  with  a frequency  response  in  excess  of  50  MHz . 


The  effective  focal-volume  diameter  is  defined  by  the  pinhole  in  the 
receiving  optics  and  is  therefore  equal  to  ~ 0.05  cm. 

Oscilloscope  traces  of  raw  signals  from  the  output  of  the  preamplifier 
at  the  base  of  the  photomultiplier  are  included  in  Figure  13.  These  photo- 
graphs are  multiple  exposures  taken  of  particle  bursts  at  the  highest 
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recorded  velocities  of  >^~2.2  (Re^  ~ 4 x 10  ) at  the  edge  of  the  boundary 
layer.  The  frequency  of  the  u-component  signal  is  about  26  MHz.  The  fre- 
quency of  the  (u  + v) -component  signal  is  almost  exactly  one-half  the  ii- 
component  frequency.  This  feature,  which  will  be  discussed  later,  is  exploited 
by  the  processing  scheme.  Typical  signal-to-noise  ratios  of  15-20  db  were 
measured  at  these  velocities. 

V . Signal  Processing  and  Data  Acquisition 

The  output  of  the  photomultiplier,  coupled  through  the  integral  pre- 
amplifier, was  monitored  on  a HP  8553B/8552B  wave  analyzer  which  was  also 
used  to  align  the  optics  and  tune  the  entire  system  by  maximizing  the  signal- 
to-noise  ratio  of  the  Doppler  signal.  This  signal  was  filtered  by  a two-pole 
low-pass  filter  in  order  to  eliminate  the  substantial  shot  noise  above  the 
Doppler  frequencies.  The  pedestal,  resulting  from  the  dual-scatter  optics, 
is  not  removed  by  this  processing  scheme. 

A.  The  LDV  Processor 

The  filtered  output  signal  is  fed  into  the  LDV  processor  (Dimotakis  and 
Lang  (1977))  depicted  in  block  diagram  in  Figure  14.  The  various  components 
of  this  processor  are  described  in  the  following  paragraphs. 

1 . The  Analog  Input  Processor 

The  Analog  Input  Processor  defines  a reference  level  Vq  and  two  symmetric 
levels  VQ  - V^,  VQ  + about  VQ  as  shown  in  Figure  15a.  In  addition,  a fourth 


level  Vu  can  also  be  defined  which  must  be  higher  than  V + V. . The 
M 0 L 

amplitudes  V^,  and  V^,  as  well  as  the  polarity,  are  independently 
selectable.  The  filtered  Analog  INput  (AIN)  must  cross  V - V from  below 

vJ  L 

and  then  from  below,  in  that  order,  before  the  Digital  INput  (DIN) 

is  raised.  A subsequent  crossing  of  from  above  will  cause  DIN  to  go  low 
This  test  sequence  must  be  passed  by  the  analog  input  AIN  on  every  cycle, 

l 

otherwise  the  burst  is  rejected,  the  logic  is  reset,  and  the  processor  waits 

1 • 

for  the  next  burst.  If  the  analog  input  AIN  exceeds  at  any  time  during 
the  burst,  the  burst  is  rejected,  the  logic  is  reset,  and  the  processor  waits 

I 

for  the  next  burst  (cf.  Figure  15a).  Thus  the  amplitudes  V_,  V,  and  Vw 

U L M 

specify  the  range  of  particle  sizes  accepted  by  the  processor  for  the  mea- 
surement. The  frequency  of  the  signal  DIN  can  be  prescaled  (cf.  Figure  15b) 
to  produce  the  signal  DIN ' which  has  one  negative  slope  for  every  m negative 
slope  of  DIN.  The  constant  m is  front  panel  selectable  and  can  be  set  1.0 
1,  2,  4 or  8.  The  Analog  Input  Processor  is  designed  with  high-speed  ECL 
circuitry  and  stripline  technology  with  a maximum  Doppler  frequency  of  170  MHz. 

2 . The  Digital  Input  Processor 

The  prescaled  digital  input  DIN ' is  subsequently  processed  by  the  Digital 
Input  Processor.  The  negative  slope  of  DIN ' fires  a resettable  one-shot  T1 
of  duration  t^,  as  illustrated  in  Figure  16.  The  negative  slope  of  T1  fires 
a resettable  one-shot  of  duration  t^.  The  first  negative  slope  of  DIN ' also 
raises  a signal  G,  which  gates  a 100  MHz  crystal  clock  into  the  Flight  Time 
(FT)  counter,  as  illustrated  in  Figure  17.  The  next  negative  slope  of  DIN7, 
refires  the  one-shot  Tl,  resets  the  T2  one-shot  and  generates  a short  TRansfer 
Pulse  (TRP)  which  stores  the  FT  buffer  and  increments  the  cycle  counter  (cf. 

Figures  16  and  17).  If  the  T2  one-shot  times  out  before  it  is  reset,  an  error 
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flag  is  raised,  the  logic  is  reset  and  the  processor  waits  for  the  next 
burst.  As  a consequence,  the  logic  checks  that  every  period  of  the  pre- 
scaled digital  input  falls  within  predetermined  limits.  If  t is  the  period 
of  the  Doppler  signal,  we  must  have,  as  a consequence, 

T x < mT  < T x + T 2 , (43) 


where  m = 1,  2,  4,  8.  The  time  intervals  T ^ and  are  front  panel  selectable 
and  cover  the  range  60  nsec  < T^,  <3  msec.  The  maximum  Doppler  frequency 

that  can  be  handled  by  the  Analog  and  Digital  Processors  working  in  tandem 
is  given  by. 


<Vmax 


max(m) 

min^) 


8 

60  X 10 


130  MHz 


(44) 


3.  The  Timer 

The  100  MHz  crystal  clock  is  divided  by  10  to  form  a 10  MHz  time 

standard  which  is  used  as  the  clock  for  the  timer,  as  illustrated  in 

N 

Figure  17.  The  timer,  in  turn,  has  a t 10  prescaler  which  sets  the 
time  base.  This  increments  a free  running  32-bit  counter  which  can  be 
reset  if  required  but  is  normally  allowed  to  overflow.  The  gate  signal 
G is  inverted  to  form  G and  is  used  to  latch  the  contents  of  the  timer 
counter  into  the  timer  buffer,  as  illustrated  in  Figure  18.  The  purpose 
of  this  information  is  to  record  the  real  time  of  each  scattering  event, 
to  permit  the  reconstruction  of  the  velocity  fluctuations  in  time  (cf. 
Dimotakis  and  Brown  (1975)).  While  the  present  data  rate  was  often  high 
enough  (~  50  kHz)  to  warrant  an  examination  of  the  velocity  fluctuations 
in  time,  it  was  decided  not  to  include  such  data  in  the  present  report. 


... 
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Performance  and 


The  LDV  processor  operates  at  the  theoretical  performance  of  ± 1/2 
clock  count  (10  nsec)  in  determining  the  flight  time  of  the  i*^ 

scattering  particle.  The  overall  measurement  accuracy  is  of  course  limited 


by  a variety  of  other  factors  such  as  optical  precision,  s ignal -to-noise 
of  the  signal  at  the  processor  input,  etc.  In  all  the  high-speed  cases 
(h^~0.8),  the  finite  frequency  of  the  processor  clock  was  the  limiting 
factor  in  the  velocity  measurement  accuracy  of  a single  particle.  The 


highest  frequency  that  occurred  was  27  MHz  in  the  2.2  flow.  This  fre- 


quency was  measured  with  an  overall  accuracy  for  a single  particle  of  ~ 27.. 


A histogram  from  a 1,024  particle  record  taken  in  the  free  stream  2.2) 

of  the  flight  time  for  8 Doppler  cycles,  is  depicted  in  Figure  19.  It  can 
be  seen  that  ~ 2/3  of  all  the  counts  fall  in  one  10  nsec  bin.  At  the  lower 
velocities  the  limitations  arose  as  a result  of  the  combined  effect  of  all 
the  other  sources  of  error.  For  the  /^,~0.1  measurements  an  accuracy  of 
0.57.  for  a single  particle  was  achieved.  It  should  bo  emphasized  that  while 
the  accuracy  for  a single  particle  velocity  measurement  does  not  limit  the 
accuracy  with  which  mean  velocities  can  be  determined,  it  contaminates  the 
measurement  of  the  rms  fluctuation  levels. 

The  processor  is  capable  of  operating  in  a variety  of  modes.  IWo  of 


the  modes  that  were  used  in  this  work  will  be  discussed. 

In  the  first  mode,  a fixed  number  of  cycles  is  stored  in  the  Mode  Logic 
circuitry  and  compared  with  the  contents  of  the  Cycle  Counter  after  each  TRP 
[Milne  (cf.  Figures  16  and  17).  When  the  two  numbers  match,  the  Count  End 
)Mi)  flax  f"  raised  which  in  turn  lowers  the  CO  disabling  the  system 


reset  and  waits  for  the  next  burst  (dead  time  ~ 60  nsec).  All  the  data 


for  0.6  (JPL  20- inch  Wind  Tunnel)  were  recorded  in  this  mode.  Eight 
cycles  were  counted  for  the  beam  pair  at  0°  (u1  data)  with  the  prescaler 


m set  to  2.  Four  cycles  were  counted  for  the  ± 45 


u data)  beam 


pairs  with  the  prescaler  set  to  unity.  The  advantage  of  this  scheme  was 


that  the  digital  processor  saw  the  same  frequencies  and  numbers  of  cycles 


for  all  configurations.  This  occurs  for  the  present  measurements  because 


the  projection  of  the  velocity  vector  along  the  ± 45  directions  is  approxi 


This  uniformity  in  the  processing,  allows  all  the  signals  to  be 


processed  with  identical  settings  except  for  a change  in  the  prescaler 


from  1 to  2.  The  particle  velocity  is  then  determined  by  reading  the 


output  of  the  FT  counter  which  corresponds  to  the  flight  time  of  the 


particle  for  the  fixed  number  of  fringes  specified 


In  the  second  mode,  the  processor  requires  that  a fixed  number  of 


cycles  be  reached  or  exceeded  by  the  burst  and  records  the  flight  time 


and  number  of  cycles  of  the  entire  burst.  The  sequence  depicted  in 


Figure  15a  represents  such  a mode  where,  assuming  that  each  period  of 


Analog  Input  Processor.  The  flight  time  At  that  would  have  resulted  from 


such  a burst  is  indicated  on  the  figure  while  the  output  of  the  cycle 
counter  would  have  been  n **  3 (number  of  TRP  pulses  and  also  number  of 
fringe  plane  intervals).  Both  n and  At  would  be  recorded  in  this  mode 
and  the  velocity  would  be 


The  data  recorded  at  0.1  (GALCIT  Merrill  Wind  Tunnel)  were  acquired 
in  this  mode  requiring  at  least  four  fringes  for  the  ± 45°  data  (u^,  u^) 
and  at  least  eight  for  the  0°  data  (u^).  Thus,  in  both  modes,  the  same 
minimum  number  of  fringes  had  to  be  crossed. 

There  are  several  important  differences  between  the  two  modes. 

First,  using  the  facility  of  the  processor  to  recognize  the  end  of  the 
burst  while  retaining  the  previous  valid  measurement  allows  the  total 
flight  time  to  be  used,  instead  of  an  arbitrary  fraction,  in  determini.ig 
the  velocity  component  of  interest.  This  improves  the  measurement  accuracy 
in  an  obvious  way.  Second,  the  phase  errors  that  are  made  in  assigning 
equal  phase  to  the  crossings  of  the  reference  level  VQ  by  the  signal  are 
cancelled  when  the  entire  burst  is  used.  Third,  the  sampling  statistics 
are  different  between  the  two  modes.  In  the  first  mode  (fixed  number  of 
fringes),  all  we  knew  is  that  the  particle  crossed  at  least  that  many.  In 
the  second  mode  (number  of  fringes  unrestricted  provided  it  exceeds  a certain 
minimum)  we  know  that  the  particle  crossed  exactly  that  many.  The  sampling 
bias  is  different  for  these  two  cases,  as  will  be  discussed  in  Section  VI-B. 

B.  Output  Formatting 

The  (i)  flight  time,  (ii)  total  number  of  validated  cycles  and  (iii) 
the  real  time  are  encoded  by  the  processor  into  three  32-bit  words  as  binary 
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coded  decimals  (BDC) . These  data  serve  as  input  to  the  subsequent  format- 
ting electronics  which  selected  the  four  least  significant  decades  of  the 
flight  time  to  form  a sixteen-bit  word  and,  depending  on  the  fringe  mode, 
either  the  sixteen  bits  of  the  real  time  for  the  fixed-fringe  mode,  or  the 
number  of  fringe  crossings,  in  the  free-fringe  mode.  These  data  were  then 
combined  to  form  a single  32-bit  word. 

The  asynchronous  data  from  the  processor  were  recorded  in  dual  1024  X 
32-bit  buffers.  One  buffer  was  filled  asynchronously  with  the  processor 
output  while  the  second  was  clocked  synchronously,  as  a sequence  of  4 bits/ 
word  on  a Kennedy  Model  9100  digital  tape  deck.  Each  1024  X 32-bit  buffer 
constituted  a 4K-bit  record  on  tape.  Several  records  were  recorded  at  each 
station  depending  on  the  turbulence  level  and  average  data  rate. 


VI.  Data  Processing 

In  computing  the  Reynolds  stress  -pu  from  the  measured  velocity  data, 
it  is  necessary  to  estimate  both  the  first  and  second  moments  of  the  data, 
and  to  extract  relatively  small  differences  between  two  large  numbers.  This 
necessitates  a careful  accounting  of  all  random  and  systematic  errors  that 
result  from  spurious  measurements.  In  the  present  experiments,  such  errors 
are  caused  by  a finite  signal-to-noise  ratio,  and  appear  as  isolated  counts 
when  the  measurements  are  treated  in  histogram  form.  Errors  are  also  caused 
by  LDV  sampling  bias  and  by  both  the  finite  clock  period  and  by  noise  in  the 
processing  electronics  and  in  the  optics.  Errors  which  arise  as  a result  of 
the  finite  measurement  volume  have  been  found  to  only  be  important  in  the 
immediate  vicinity  of  the  wall,  and  have  not  been  corrected  in  the  present 
analysis.  The  correction  of  such  errors  would  result  in  a correction  of  at 
most  the  first  two  data  points  adjacent  to  the  wall,  for  which  we  have  the 
least  confidence,  and  would  not  affect  the  conclusions  presented. 
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A.  Histogram  Pruni 


The  data  from  the  present  experiments  were  recorded  as  multiple  records 
of  1024  individual  realizations  of  the  three  instantaneous  velocity  com- 
ponents u,  u + v,  u - v.  These  data  were  processed  by  forming  a separate 
histogram  for  each  record,  and  subsequently  performing  the  statistical 
computations  using  the  resulting  histograms.  For  the  high-speed  flow  experi- 
ments, the  data  were  recorded  in  the  fixed-fringe  mode  and  the  histograms 
were  formed  as  a function  of  the  integer  flight  time  At^.  The  low-speed 
measurements  were  recorded  in  the  free-fringe  mode  for  which  the  flight 
time  was  large  and  covered  a much  wider  range.  For  these  data,  the  histogram 
was  formed  as  a function  of  the  velocity  of  the  particles,  computed  as  an 
integer  percentage  of  the  maximum  velocity  that  the  processor  would  admit, 
i.e . , 

Ii  - 100  tx  n.  , (48) 


where  T ^ is  the  minimum  Doppler  period  and  n^  is  the  number  of  fringes 
crossed  by  the  i**1  particle  in  time  At^. 

The  pruning  was  performed  in  two  passes.  First,  data  in  any  bin  that 
contained  only  one  count,  and  did  not  have  neighboring  bins  with  more  than 
one  count,  were  discarded.  In  the  second  pass,  any  data  isolated  from  the 
main  body  of  the  histogram  by  at  least  one  zero  were  also  discarded.  If, 
as  a result  of  these  two  operations,  more  than  24  measurements  out  of  the 
total  of  1024  were  rejected,  the  entire  record  was  rejected.  This  scheme 
was  preferred  to  the  more  usual  one  of  rejecting  data  outside  a certain 
number  of  standard  deviations  from  the  mean.  By  way  of  example,  the  few 
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counts  in  bin  40  of  the  histogram  in  Figure  19  would  be  discarded  by  this 
scheme . 


The  histograms  which  result  from  the  data  in  the  present  investigations 
are  illustrated  in  Figures  20  and  21.  In  these  figures,  the  histograms  at 
several  locations  through  the  boundary  layer  are  superimposed  on  the  mean- 
velocity  profiles  for  that  case.  For  both  flows  represented  here,  the 
histogram  for  the  free-stream  flow  encompasses  only  four  bins.  Near  the 
wall,  the  histogram  is  spread  over  more  than  forty  bins. 

B.  Sampling  Bias 

The  complexity  of  processing  single  particle  laser-Doppler  velocity 
data  is  compounded  by  the  fact  that  the  fluid  velocity  is  sampled  in  a 
biased  manner.  This  problem,  first  pointed  out  by  McLaughlin  and  Tiederman 
(1973),  is  a consequence  of  the  fact  that  the  particle  flux  through  the 
measurement  volume,  and  the  resulting  measurement  probability  per  unit  time, 
is  higher  when  the  local  fluid  velocity  is  high  than  when  it  is  low.  The 
fact  that  the  measurement  probability  is  a function  of  the  measured  quantity 
gives  rise  to  the  sampling  bias  0(u),  where  u is  the  three-dimensional  velocity 
vector.  Consequently,  to  compute  the  expectation  value  of  a particular 
function  of  the  measurements  we  must  remove  this  bias,  i.e.. 


<1 


4 


(49) 


The  dependence  of  the  sampling  bias  function  3 on  the  local  velocity  vector 
u,  on  the  shape  of  the  measurement  volume,  and  on  the  minimum  number  of 
fringe  crossings  required  by  the  processing  electronics,  has  been  derived 
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by  Diraotakis  (1976).  For  data  recorded  in  the  fixed-fringe  mode,  if  the 
angle  between  the  intersecting  beams  is  small,  the  expression  for  0(u) 
simplifies  to 

« . ~ •>  2 ■ 

P(u;e)  “■  u^(l  + u„  /u^  ) [1  - e (1  + u„  /u^  )]  » (50) 

where  is  the  velocity  component  perpendicular  to  the  fringe  planes, 
is  the  velocity  component  parallel  to  the  fringe  planes  and  perpen- 
dicular to  the  beam  bisector,  and  e is  the  ratio  n . /n»  of  the  minimum 

ram  T 

number  of  fringe  crossings  required  by  the  processor  to  the  total  number 

in  the  focal  volume.  Note  that  recording  the  flight  time  for  a fixed 

number  of  fringes  in  fact  requires  that  at  least  that  many  were  crossed. 

Consequently,  the  fixed  number  of  fringes  is  equal  to  n^  in  the  context 

in  which  e was  defined  for  Equation  (50). 

Evaluating  f3(u;e)  for  the  u^  data  in  the  boundary  layer,  which  were 

acquired  using  the  fixed-fringe  mode  (0.6  ~ 2.2),  we  note  that 

2 2 2 2 

u */u  “ v /u  « 1 so  that,  in  this  case, 

* J. 

6(u;e)  ~u±*(l  - e2)  , (51) 

or, 

P(u;e)  « ux  • (52) 

In  going  from  Equation  (51)  to  Equation  (52),  an  additional  assump- 
tion has  been  made.  Even  though  n , is  fixed  by  the  signal  processing 

rain 

electronics,  the  total  number  of  fringes  crossed,  n^,  is  unfortunately  a 

complicated  function  of  the  particle  size.  This  is  a consequence  of  the 

fact  that  the  amplitude  levels  V.  and  V are  absolute.  A large  particle, 

U 1» 
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The  latter  form  is  preferable  in  avoiding  finite  precision  problems  in 
the  computer,  at  the  expense,  however,  of  substantially  more  computing 
time. 


In  evaluating  Equations  (54)  and  (55),  note  that  the  i^  measured 
velocity  component  is  given  by  Equation  (47),  or 


(58) 


where  n is  equal  to  the  number  of  fringes,  s is  the  fringe  spacing,  is 
the  clock  period  (10  nsec),  kj  is  the  (integer)  output  of  the  flight  time 
counter  and  x is  a Tandom  variable  representing  the  collective  uncertainty 
in  the  determination  of  the  flight  time  At  and  resulting  from  (i)  the 
finite  clock  period  , and  (ii)  the  uncertainty  with  which  the  equal 
phase  points  are  determined  by  the  comparators  from  the  signal  burst  as 
a result  of  finite  signal-to-noise  ratio.  The  1/2  is  added  to  k.^  to  remove 
the  systematic  truncation  error  of  the  digital  counter. 

Using  Equations  (54)  and  (55)  and  correcting  for  the  finite  variance 
of  the  random  variable  x,  we  have 


«■*>«--&)  XT  • 

<ul2)flx.d  “(f)  ( 2 J.,  3)  • 
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where 


<ki  +b~l 


(61) 


and 


2 

CT 


J x2p(x)dx 
J p(x)dx 


(62) 


is  the  variance  of  x. 

Unfortunately,  the  situation  is  less  well  defined  for  the  ± 45°  data 
recorded  in  the  fixed-fringe  mode.  In  that  case,  the  two  relevant  veloci- 
ties, u^  and  u^,  are  of  the  same  order  of  magnitude.  Nevertheless,  since 
only  one  velocity  component  was  measured  at  a time,  we  use  Equation  (52) 
for  the  bias  and  reduce  the  ± 45°  data  using  Equations  (59)  and  (60).  While 
this  is  clearly  not  correct,  it  is  better  than  ignoring  the  bias  altogether. 

Note  that  the  expression  for  the  sampling  bias  0 given  by  Equation  (52) 
is  not  appropriate  for  data  acquired  in  the  free-fringe  mode  because  this 
equation  was  derived  using  the  assumption  that  the  probability  of  measure- 
ment is  determined  by  the  fact  that  the  particle  has  crossed  at  least  the 
fixed  number  of  fringes.  In  the  free-fringe  mode,  the  measurement  of  the 
number  of  fringes  crossed  and  of  the  flight  time  in  units  of  corresponds 
to  a particle  crossing  of  exactly  n^  fringes.  The  probability  for  this  event 
can  be  determined  from  Equation  (50).  Since  0 (u;  n^/n^,)  is  proportional  to 
the  probability  of  crossing  at  least  n^.  fringes,  i.e., 

Pfixed  (^W  « P £n  < i^i  | u } , (63) 
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we  have  that 


p [n  = ni|  u}  *=  p [ti  < riil  u}  - p {n  C (r^  + 1)  | u } , 


where  p { } denotes  the  probability  per  unit  time  of  observing  the  event 
in  the  braces.  Consequently, 


Pfree^W  * Pf  ixed^W  ' Pf  ixed1^'*  (ni  + 1)/nT]  * (65) 


Pfree^— ;"i^nT^  “ ux"2(uJ-2  + u"2)  (ni  + 1/2)  » 


where,  in  the  context  of  Equation  (49),  the  assumption  has  been  made  that 
the  particle  size  is  uncorrelated  with  the  local  velocity.  In  the  context 
of  Equation  (49),  note  that 


whereas 


79  1 

(UjL  + U/,  ) « 7—  . , 


where  &t±  is  the  time  of  flight  for  the  measured  nt  fringes.  Substituting 
these  equations  in  Equation  (66) , we  have 

, (kt  + 1/2  + x)n.2 

sfr«  fe'W  “ „"+'i7i  ' 


where,  as  before,  kt  is  the  integer  output  of  the  flight  time  counter  and 
x is  a random  variable  representing  the  uncertainty  in  the  measurement  of 


the  flight  time.  Using  this  bias  function,  we  thus  have 


-ft)' Hr11)  • 


where 


K*,„  B £ 


(k.  + |)A(ni  + j) 


and  o is  the  variance  of  x.  Note  that  Equations  (70)  and  (71)  are  equally 
valid  for  the  0°  data  as  for  the  ± 45°  data.  The  data  recorded  at  ^=“0.1 
in  the  free-fringe  mode  were  analyzed  using  Equations  (70)  and  (71).  In 
order  to  compare  them  to  the  previous  data,  Equations  (54)  and  (55)  were 
also  used,  which  in  this  case  reduce  to 


where 


(+J(;iJerr“)  • 

Em 

i (kt  + j)£ 


Note  also  that  Equation  (69)  for  the  sampling  bias  in  the  free-fringe 

mode  is  in  disagreement  with  the  results  of  Hoesel  and  Rodi  (1977)  and 

Buchhave,  George  and  Lumley  (1979),  who  claim  that  a fit,  a (k.  + ^-)  . 

free  l l 2 

If  that  were  true.  Equations  (59)  and  (60)  would  be  valid  for  the  free- 

fringe  mode  as  well,  with  the  three  summations  as  defined  in  Equation  (61). 

2 

The  variance  o of  x was  estimated  from  the  mean-square  fluctuations 

/2 

in  the  free  stream  where  (u  ; ) =“  0. 

The  results  of  this  analysis  have  been  applied  to  the  data  obtained 
in  these  experiments.  The  bias  correction  is  most  important  in  the  region 
near  the  wall,  where  the  fluctuation  levels  are  highest.  These  calculations 
indicate  at  most  a 27.  correction  for  the  low-speed  flow  data  in  this  region, 
and  a 3.3%  correction  for  the  high-speed  flow  data. 

VII.  Results  and  Discussion 
A.  Data  Reduction 

Following  the  computation  of  the  proper  statistical  averages,  the  data 
were  reduced  to  a form  appropriate  for  boundary  layer  flows  by  utilizing  the 
scheme  outlined  in  Part  I of  this  report  by  Collins,  Coles  and  Hicks  (1978). 

For  this  purpose,  the  Van  Driest  (1955)  scaled  mean  streanwise  velocity  data 
were  fitted  to  the  law  of  the  wall  and  the  law  of  the  wake  using  a three- 
parameter  nonlinear  least-squares  fit  with  a single  constraint  equation  based 
on  the  local  friction  law.  The  profile  parameters  u^,  n and  6,  determined 
by  this  analysis,  were  then  used  to  infer  both  the  distribution  of  shear  stress, 
and  of  the  normal  velocity,  from  the  data  by  an  integration  of  the  equations 
for  the  conservation  of  mass  and  momentum  from  the  wall  to  the  free  stream, 
using  the  formulation  by  Spalding  (1961)  to  describe  the  flow  in  the  viscous 
sublayer.  The  details  of  these  computations  are  given  in  Part  1. 
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One  significant  difference  between  the  present  LDV  data  and  the  data 
obtained  using  standard  Pitot  instrumentation  (cf.  Part  I)  is  the  ability 


of  the  Pitot  to  obtain  valid  data  in  the  region  adjacent  to  the  wall. 

Because  of  noise  considerations,  the  present  data  are  limited  to  a distance 
of  approximately  one  focal-volume  diameter  from  the  wall.  Inis  fact  requires 
that  the  boundary-layer  scale  parameters  6 and  0 be  computed  directly  from 
the  law  of  the  wall  and  the  law  of  the  wake  representations  rather  than  from 
a direct  integration  of  the  data,  as  in  the  results  discussed  in  Part  I. 

The  scale  parameters  are  defined  as 
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dy  , (76) 
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In  terms  of  the  Van  Driest  (1955)  scaled  velocity  U, 
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and  U is  defined  by  the  Van  Driest  scaling 
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Thus, 
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The  integrals  involved  in  these  expressions  were  evaluated  in  terms 
of  the  integrals  P and  Q defined  in  Part  I.  These  integrals  are  functions 
only  of  the  mean-velocity  profile  parameters  u^,  H and  6,  and  are  evaluated 
in  terms  of  the  Spalding  formulation  for  the  boundary  layer. 

As  a result  of  the  differences  in  the  computational  procedure  for  the 
scale  parameters  6*  and  8,  these  quantities  differ  slightly  from  those 
computed  for  the  identical  flows  given  in  Part  I of  this  work.  Corrections 
for  these  discrepancies  have  been  included  where  comparisons  are  made  to 
the  previous  work. 

B.  Streamwise  Mean-Velocity  Profiles 

A comparison  between  the  mean-velocity  profiles  obtained  in  the  present 
experiments  and  those  obtained  from  the  Pitot  data  of  Part  I is  given  in 
Figures  22  and  23.  In  order  to  make  a consistent  comparison,  the  Pitot  data 
from  Part  I have  been  replotted  using  the  integral  scale  0 obtained  from  the 
LDV  experiments.  At  high  Mach  numbers,  the  results  exhibit  a small  discrep- 
ancy between  the  LDV  and  the  Pitot  data,  with  the  LDV  data  exhibiting  lower 
mean  velocities  near  the  wall.  The  cause  of  this  discrepancy  is  not  clear 
at  this  writing,  nor  is  it  clear  which  measurement  is  the  more  reliable. 

The  low-speed  flow  mean-velocity  profiles  have  been  computed  using  both  the 
exact  bias  correction,  according  to  Equation  (70)  and  using  the  harmonic 
mean  according  to  Equation  (59).  The  differences  are  small  as  are  the 
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differences  between  these  results  and  the  results  of  the  mean  value  of 
these  data  without  any  bias  correction.  The  comparisons  have  not  been 
plotted  here.  The  agreement  between  the  LDV  data  and  the  Pitot  data  is 
excellent  for  the  low-speed  results . 

The  data  obtained  from  the  present  experiments  have  been  plotted 
in  coordinates  appropriate  to  the  law  of  the  wall  in  Figures  24  and  25. 

The  solid  line  associated  with  each  data  set  has  been  computed  using  the 
formulation  by  Spalding  (1961),  and  illustrates  the  fit  obtained  in  the 
computation  of  u^,  II  and  6.  The  high-speed  flow  results  shown  in  these 
figures  illustrate  the  problem  of  obtaining  measurements  near  the  wall. 
These  data  exhibit  a departure  from  the  computed  profile  because  of  noise 
which  occurs  when  the  focal  volume  intercepts  the  wall.  This  fact  required 
that  the  data  included  in  the  least-squares  fit  be  limited  to  y+  ^ 200  for 
these  data,  as  was  the  case  for  the  Pitot  data  reported  in  Part  I. 

The  profile  parameters  u^  and  6 computed  from  the  least-squares 
procedure  agree  closely  with  those  parameters  computed  from  the  Pitot  data 
for  the  same  flows.  The  computations  of  the  wake  parameter  II,  however,  are 
considerably  higher  in  every  case.  The  self-consistency  of  the  LDV  results 
indicates  that  this  discrepancy  is  probably  not  the  consequence  of  data 
reduction  errors. 

C.  Streamwise  Velocity  Fluctuations 

The  fluctuations  of  the  streamwise  velocity  component,  normalized  by 
the  free-stream  velocity,  are  shewn  in  Figures  26  and  27  as  a function  of 
the  normal  coordinate.  These  results  are  compared  to  the  distribution  of 
the  u7  fluctuations  obtained  at  low  speeds  by  Klebanoff  (1954)  in  flow  at 
constant  pressure,  represented  by  the  solid  squares  in  these  figures.  For 


these  comparisons,  the  normal  coordinate  for  the  data  by  Klebanoff  have 
been  scaled  with  the  boundary-layer  scale  obtained  in  the  present  experi- 
ments . 

The  error  bars  in  each  of  the  subsequent  figures  have  been  computed 
as  the  standard  deviation  of  the  ensemble  of  records  used  in  computing 
the  mean  values.  In  each  case,  a few  data  points  for  which  the  error 
bars  exceed  the  mean  value  have  been  excluded.  From  these  results,  it 
is  clear  that,  near  the  free-stream,  the  velocity  fluctuations  are  properly 
computed  from  the  LDV  data.  As  the  wall  is  approached,  the  deviations 
between  the  present  results  and  those  of  Kelbanoff  are  within  the  confidence 
limits  in  the  present  data. 

It  can  be  concluded  from  these  data  that,  within  the  error  in  the 
present  experiments,  no  measurable  effects  of  compressibility  have  been 
found  on  the  mean-square  velocity  fluctuations.  The  data  for  Mach  numbers 
up  to  2.2  are  adequately  represented  by  the  data  of  Klebanoff.  The  differences 
exhibited  in  these  figures  between  the  high-speed  flow  data  and  the  low-speed 
flow  data  arise  from  the  differences  in  the  boundary- layer  Reynolds  number, 

Reg.  These  results  confirm  Morkovin's  (1961)  hypothesis  that  there  is  no 
essential  difference  in  the  dynamic  behavior  of  the  boundary  layer  at  constant 
pressure  for  Mach  numbers  up  to  4 or  5. 

D.  Reynolds  Stress  Profiles 

The  distribution  of  Reynolds  stress,  computed  as  - pu 'v is  shown  in 
Figures  28  and  29  for  the  data  of  these  experiments.  In  order  to  compare 
the  Reynolds  stress  deduced  from  the  experimental  data  with  the  expected 
distribution  through  the  boundary  layer,  a second  computation  based  on  the 
integral  formulation  outlined  in  Part  I is  displayed  as  the  solid  curve  in 
each  figure.  The  results  from  the  integral  formulation  coincide  with  the 
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results  obtained  using  the  Pitot  data,  with  minor  differences  arising 
from  the  difference  in  the  computed  value  of  the  wake  parameter  II. 

For  comparison  with  the  expected  behavior  at  lew  speeds,  the  data 
of  Klebanoff  (1954)  has  been  included  in  these  figures  as  the  dark  squares 
on  each  figure.  The  agreement  between  the  integral  formulation  and  the 
data  of  Klebanoff  indicates  that  the  integral  formulation,  based  on  the 
mean-velocity  data,  yields  the  expected  distribution  for  the  total  shear 
stress . 

The  dominant  feature  of  the  results  for  the  distribution  of  the 
Reynolds  stress,  as  measured  by  the  laser-Doppler  technique,  is  the  de- 
parture of  the  measured  Reynolds  stress  from  the  results  of  the  integral 
computation  in  the  region  y/S  ~ 3.  This  departure  represents  a major  error 
in  these  measurements  and  may  indicate  a fundamental  limitation  on  the 
applicability  of  the  laser-Doppler  technique  in  flows  of  this  type. 

Several  conclusions  can  be  reached  from  these  data.  The  anomaly 
represented  by  the  departure  of  the  Reynolds  stress  from  the  expected  value 
near  the  wall  has  been  reported  previously  by  Johnson  and  Rose  (1973),  by 
Yanta  and  Lee  (1974),  and  by  Abbiss  (1976).  An  explanation  of  this  phe- 
nomenon by  Sandborn  (1974)  assumed  that  the  phenomenon  is  related  to  the 
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of  the  mean-square  velocity  fluctuations  indicates  that  the  particles 
follow  the  flow  locally  throughout  the  boundary  layer. 

Similar  distributions  for  the  Reynolds  stress,  measured  using  hot- 
wire anemometer  techniques  by  Laderman  (1978),  have  recently  been  shown 
to  agree  with  the  present  results  given  by  the  integral  distributions 
shown  by  the  solid  curves  in  Figures  28  and  29. 

These  results  again  confirm  Morkovin's  hypothesis  (1961)  and  refute 
the  hypothesis  by  Sandborn  (1974)  that  density  fluctuations  may  be  re- 
sponsible for  major  changes  in  the  distribution  of  Reynolds  stress  at 
high  speeds.  Similar  conclusions  were  expressed  in  Part  I,  based  on 
the  results  of  the  Pitot  measurements  of  the  mean-velocity  profiles. 

E.  The  Problem  of  the  Normal  Velocity 

The  results  for  the  distribution  of  normal  velocity  are  illustrated 
in  Figure  30  for  the  high-speed  flow  at  low  Reynolds  number.  In  this 
figure,  the  solid  line  is  the  result  for  the  correct  distribution  of 
v/u,  computed  from  the  integral  representation  of  the  data  outlined  by 
Collins,  Coles  and  Hicks  (1978),  using  the  boundary- layer  parameters  , 
n and  6 derived  from  the  least-squares  analysis. 

The  errors  in  estimating  the  normal  velocity  from  the  present  data 
are  substantial  in  every  case.  The  principal  feature  of  the  data  is  that 
the  normal  velocity  is  large  and  negative  at  the  wall,  increasing  sharply 
in  magnitude  as  the  wall  is  approached.  The  disparity  between  the  large 
negative  value  and  the  error  bars  indicates  that  the  particles  are  being 
convected  toward  the  wall  in  this  region.  In  the  free  stream,  the  velocity 
is  positive  in  most  cases,  but  with  a magnitude  that  is  larger  than  the 
correct  free-stream  value  as  represented  by  the  integral  formulation. 


These  data  illustrate  the  difficulty  in  accurately  representing 
the  normal  velocity.  This  difficulty  arises  as  a consequence  of  the 
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fact  that  any  small  misalignment  of  the  optics  will  result  in  substantial 
errors  in  the  computed  value  of  the  normal  velocity,  which  has  a maximum 
value 


(v/u)  ~ 
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(84) 


where  H * 6*/0  is  the  boundary  layer  form  parameter.  In  addition,  because 
of  the  non-uniform  particle  distributions,  the  normal  velocity  components 
of  both  the  flow  near  the  wall  and  the  flow  in  the  vicinity  of  the  boundary- 
layer  edge  are  improperly  represented. 

The  fluctuations  in  the  normal  velocity  are  represented  for  these  same 
data  in  Figure  31.  These  data  have  approximately  the  value  of  the  data  of 
Klebanoff  (1954)  for  incompressible  flow,  represented  by  the  solid  squares 
in  the  figure.  However,  in  light  of  the  difficulties  in  measuring  the 
normal  velocity,  this  agreement  is  fortuitous.  The  magnitude  of  the  error 
bars  for  these  data  are  approximately  one-half  of  the  magnitude  of  the  data, 
and  indicate  a low  degree  of  confidence  in  the  results. 

VIII.  Conclusions 

The  present  experiments  lead  to  a number  of  conclusions  regarding  the 
applicability  of  the  laser-Doppler  technique  to  the  study  of  boundary- layer 
flows.  The  results  presented  for  the  behavior  of  the  streamwise  velocity 
components  indicate  the  importance  of  including  in  the  analysis  of  the  raw 
data  an  adequate  technique  for  pruning  the  initial  histograms  in  order  to 
remove  errors  in  the  data  introduced  by  noise.  These  results  further  indicate 
the  importance  of  including  proper  statistical  averaging  in  the  analysis  in 
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order  to  account  for  the  bias  introduced  by  sampling.  In  addition,  the 
effects  of  a finite  clock  period  and  noise  have  been  shown  to  be  important 
in  the  proper  determination  of  the  streanwise  velocity  fluctuations,  u /, 
near  the  boundary- layer  edge.  When  these  effects  are  properly  included 
in  the  data  analysis,  a direct  comparison  between  the  present  laser-Doppler 
measurements  and  measurements  in  the  same  flows,  presented  in  Part  1 of 
this  work,  indicates  good  agreement  for  the  streanwise  velocity  components 
u and  u/.  The  present  results,  when  properly  scaled  to  account  for  com- 
pressibility using  the  ideas  of  Van  Driest  (1955),  are  also  shwn  to  agree 
with  the  low-speed  data  presented  by  Kelbanoff  (1954).  This  observation 
further  corroborates  the  conclusions  expressed  in  Part  J of  this  work 
regarding  the  proper  role  of  compressibility  in  determining  the  behavior 
of  the  mean  velocity,  and  confirms  Morkovin's  hypothesis  (1961)  regarding 
the  role  of  density  fluctuations  in  the  boundary  layer  at  constant  pressure. 

The  present  measurements  of  the  Reynolds  stress,  -pu/v/,  and  for  the 
normal  velocity,  v,  are  in  substantial  disagreement  with  the  expected  results 
derived  from  the  mean  flow.  These  observations  seem  to  indicate  a funda- 
mental difficulty  with  the  application  of  the  laser-Doppler  technique  as 
implemented  in  these  experiments  to  the  measurement  of  the  detailed  structure 
of  flows  in  the  vicinity  of  a wall.  The  measured  Reynolds  stress  profiles, 
-pu'v',  deviate  from  the  expected  behavior  for  y/0  ~ 3,  indicating  an  error 
in  the  measurement  for  the  u 'v  ' correlation  in  the  region  adjacent  to  the 
wall.  This  error  arises  because  of  the  strong  correlation  between  the  normal 
velocity  component  and  the  particle  number  density  that  occurs  as  a result 
of  the  depletion  of  particles  near  the  wall.  These  results  indicate  that 
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the  flow  inside  the  viscous  sublayer  is  such  that  a particle  that  enters 

the  sublayer  has  a very  low  probability  of  leaving.  The  scaling  law  for 

particle  behavior  in  a turbulent  boundary  layer  (cf.  Equation  (13))  is  not 

applicable  in  the  viscous  sublayer.  A particle  which  enters  the  viscous 

sublayer  is  subjected  to  the  very  high  Lagrangian  frequencies  which  occur 

in  the  motion  of  the  longitudinal  vorticies,  whose  transverse  extent  is 
20v 

approximately  — , and  lagp  the  flow  in  this  region.  As  a consequence, 
a fluid  element  coming  from  the  wall  is  less  likely  to  carry  particles  than 
a fluid  element  moving  toward  the  wall.  This  results  in  an  under- 
representation of  the  positive  component  of  the  normal  velocity,  v,  in 
the  region  near  the  wall,  resulting  in  a negative  value  for  the  measured 
component  normal  to  the  wall. 

This  under-representation  of  the  motion  near  the  wall  also  results 
in  a lower  measured  value  for  the  Reynolds  stress,  -pa'v',  than  would  be 
expected  by  an  analysis  based  on  the  mean  flow  as  described  in  Part  I,  since 
the  upswelling  of  fluid  from  the  wall  is  associated  with  a lower  streanwise 
velocity.  This  bursting  motion  has  been  observed  by  Blackwelder  and  Kaplan 
(1976)  to  be  intermittent  and  quite  violent,  and  is  held  to  be  responsible 

/y°T  < \ 

for  a large  fraction  of  the  total  stress  near  the  wall  l — 1001.  The 

particle  transport  to  the  wall  via  the  viscous  sublayer  is  so  effective, 
even  at  the  lower  velocities  of  the  Merrill  wind  tunnel  experiments,  that 
the  wall  appears  as  an  infinite  sink  for  particles. 

We  conclude  from  these  observations  that  the  measured  disparity  between 
the  Reynolds  stress,  -pu'v7,  and  the  total  stress,  as  defined  by  the  inte- 
gration of  the  mean  flow,  arise  not  as  a consequence  of  compressibility, 
as  suggested  by  Sandborn  (1974),  but  as  a consequence  of  particle  depletion 
which  is  unique  to  boundary  layers  in  air.  This  is  a fundamental  limitation 
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of  Che  techniques  for  particle  seeding  which  must  be  solved  if  this 
technique  is  to  be  useful  for  detailed  measurements  near  the  wall. 

For  free>shear  flows,  because  of  the  absence  of  the  wall,  this  problem 
does  not  arise,  and  the  laser-Doppler  technique  will  give  an  accurate 
representation  of  the  mean  and  fluctuating  velocity  components  and  their 
correlations . 
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Figure  1.  The  Doppler  Frequency  as  a Function 
of  Flow  Velocity 


i 


5 U 


Figure  5.  The  Ratio  of  the  Boundary -Layer  Edge  Velocity  to 
the  Boundary-Layer  Thickness  in  Wall  Coordinates 
as  a Function  of  Reg. 
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Figure  7.  The  Distribution  of  Particle  Number  Density. 
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Figure  13. 

Laser-Doppler  Anemometer  Signal 
Turbulent  Boundary  Layer  - 2 
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Figure  15.  The  Analog  Processing  Technique. 
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Figure  17.  The  Digital  Input  Processor. 


Figure  18.  The  Timing  Circuit. 
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Figure  19.  A Typical  Histogram  from  the  Free-Stream 
Flow  at  M «*■  2.2. 
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Figure  20.  Representative  Histograms  from  the  High- 
Speed  Boundary  Layer. 
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Figure  26.  The  Streanwise  Velocity  Fluctuations 
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Figure  29.  The  Distribution  of  Reynolds  Stress 
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Nomenclature 


Symbol 

Equation 

Meaning 

c 

(22) 

constant  in  vail  law  (5.0) 

Cf 

(84) 

local  friction  coefficient 

d 

P 

(7) 

particle  diameter 

ki 

integer  output  of  the  flight  time 
counter 

n 

(21) 

particle  number  density 

ni 

number  of  observations  with  velocity 

Keg 

Reynolds  number  based  on  6 

6 

(2) 

fringe  spacing 

U,  V 

streamwise  and  normal  velocity  components 

ui 

(29) 

velocity  component  in  itk  direction 

UT 

(18) 

friction  velocity 

U 

JL 

(1) 

velocity  component  normal  to  the 
fringe  plane 

u> 

velocity  component  parallel  to  the 
fringe  plane 

x»  y 

streaowise  and  normal  coordinates 

er,  8,  V.  0 

Figure  12 

angles  defined  by  the  beam  planes 

P(u) 

(49) 

bias  function 

6 

(26) 

boundary-layer  thickness 

e 

(77) 

boundary- layer  momentum  thickness 

ei 

Figure  12 

beam  separation  angles 

K 

(22) 

Karman  constant  (0.41) 

X 

laser  wave  length  (0.5145  ^m) 

n 

viscosity 

V 

kinematic  viscosity 

nomenclature  (Cont.) 


Symbol 

Equation 

Meaning 

*0 

(1) 

Doppler  frequency 

n 

(22) 

strength  of  the  wake  component 

T 

C 

(5) 

processor  clock  period 

T 

(8) 

particle  response  time 

P 

T 

shear  stress  at  the  wall 

w 

Subscripts 

<>e 

edge  or  external  value 

pertaining  to  the  fluid 

C ) 

pertaining  to  the  particles 

P 

( )w 

wall  value 

( >' 

ms  fluctuations 

Superscripts 


value  made  dimensionless  with 
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Appendix  G 


Two-Point  LDV  Measurements  in  a Plane  Mixing  Layer 

# # **  * 

M.  M.  Koochesfahani  , C.  J.  Catherasoo  , P.  E.  Oimotakis  , M.  Gharib 

and  D.  B.  Lang* 

California  Institute  of  Technology,  Pasadena,  California 

Investigations  into  the  nature  of  the  large  structures  in  a two- 
dimensional  shear  layer  were  carried  out  using  Laser  Doppler  velocimetry 
in  the  GALCIT  Free  Surface  Water  Tunnel.  By  simultaneous  measurements 
of  velocity  at  two  points  outside  the  turbulent  region,  above  and  below  the 
shear  layer,  it  was  possible  to  measure  the  strength  (total  circulation) 
and  the  location  of  the  vorticity  center  of  the  large  structures.  It  was 
found  that  structures  not  in  the  process  of  pairing  convect  downstream 
with  the  center  of  their  cores  close  to  the  ray  y/x  along  which  the  mean 
velocity  Um  = ^(Ui  + Uj).  The  determined  value  of  the  mean  circulation 
is  consistent  with  the  independent  measurements  of  the  mean  spacing 
between  the  structures.  Results  indicate  that  if  the  large  structure 
vorticity  distribution  is  elliptical,  the  inclination  angle  of  its  axis  of 
symmetry  with  respect  to  the  flow  direction  is  small. 
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1 . Introduction 


In  a recent  experiment*,  the  induced  streamwise  velocity  fluctuations 
due  to  the  passage  of  vortical  structures  in  a two  dimensional  shear  layer 
were  measured  at  a point  just  outside  the  mixing  layer.  Such  a velocity 
record,  however,  (see  figure  1),  provides  ambiguous  information  about 
the  strengths  and  the  core  positions  of  the  vortical  structures,  since  the 
induced  velocity  due  to  the  passage  of  a strong  vortex  far  from  the  measur- 
ing station,  is  indistinguishable  from  that  due  to  a weaker  vortex  passing 
close  to  the  measuring  station. 

To  resolve  this  ambiguity,  simultaneous  measurements  at  two  points 
were  performed,  above  and  below  the  shear  layer,  to  allow  the  circula- 
tions (strength)  and  core  position  of  the  vortical  structures  to  be  determined. 


2.  Apparatus,  instrumentation  and  data  processing 
a.  Apparatus  and  the  resulting  flow  field 

The  flow  apparatus  and  the  resulting  shear  layer  have  been  docu- 
mented in  Ref.  1. 

From  the  measured  mean  velocity  profiles  (Ug  /Ux  “*  0.  19),  the 
vorticity  (maximum  slope)  thickness  6^,  normalized  with  the  distance 
downstream  from  the  splitter  plate  x,  was  found  to  be 


6 

“■  0.  10  . 
x 


(1) 


The  data  presented  in  this  paper  were  taken  at  c = 30  cm  with 
Uj  = 56  cm/s.  These  parameters  yield  a Reynolds  number,  defined  by 


AU6 


u> 


, of  1.  4 x 10^  at  the  measuring  location,  where  AU  is  the  velocity 


difference  across  the  mixing  layer.  Note  that  the  shear  layer  is  fully 
turbulent  at  this  Reynolds  number. 
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b.  L.DV  optics  and  electronics 

The  streamwise  component  of  velocity  was  measured  at  the  two 
points  by  a single-particle  Laser -Doppler  velocimeter  in  the  dual  scatter 
mode.  The  two  pairs  of  beams  were  generated  by  passing  the  laser  beam 
through  two  cube  beam-splitters  at  right  angles  (fig.  2).  The  two  focal 
volumes  (measuring  stations)  were  located  at  the  high-speed  and  low- speed 
edges  of  the  mixing  layer,  respectively  (see  fig.  3). 

The  outputs  of  the  two  photodetectors  were  time -multiplexed  at  the 
input  of  a single  LDV  processor.  The  processor  has  been  described  in 
Ref.  1. 


c.  Data  processing 

The  velocity  data  were  low-pass  filtered  digitally  according  to  the 
algorithm 


tyt_)  = y h(t  - t )U(t  )6t 
m m n n n 


n 


6t  = t - t , , 
n n n- 1 


(2a) 

(2b) 


where  h is  the  (normalized)  filter  function,  U the  unfiltered,  randomly 

sampled,  velocity  data  and  O’  is  the  filtered  velocity  data.  The  summation 

over  n is  such  that  -T,  s t - t s T,  where  2Tf  is  the  total  width 

f m n f f 

of  the  filter  h(t). 

A Gaussian,  low-pass  filter  was  used.  The  filtered  data  were 
sampled  uniformly.  See  fig.  4. 

3.  Data  Reduction 
a.  Model 

Recent  measurements  have  shown*'  ^ that  the  two- 

dimensional  shear  layer  is  characterized  by  an  array  of  large  vortical 
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structures,  which  are  convected  downstream  with  a velocity  U£, 
approximately  given  by  Uc  “■  |(Ut  + Uq  ).  Now,  in  a frame  moving 
with  Uc,  the  induced  velocities  u^,  u^,  at  the  points  (x,  b)  and 
(x,  -b)  respectively,  when  the  center  of  vorticity  of  a particular 
structure  is  at  (x,  c),  are  given  by 


_ r L_ 

UH  “ 2tr  b-c 


(3a) 


u 


L = 


r i_ 

2tt  b+c 


(3b) 


where  T is  the  total  circulation  of  the  vortical  structure.  Note  that 
F is  defined  to  be  positive  for  an  induced  velocity  in  the  clockwise 
sense.  Note  also  that  the  absolute  values  of  the  induced  velocities 
reach  a maximum  at  the  instant  the  center  of  vorticity  is  on  the  line 
joining  the  measuring  points.  In  this  model  the  effect  of  neighboring 
structures  is  neglected.  It  can  be  shown,  however,  that  the  resulting 
overestimation  of  the  circulation  T is  of  the  same  order  as  the  under- 
estimation of  r caused  by  neglecting  the  expected  higher  moments 
(quadrupole  moment)  of  the  vorticity  distribution.  The  effect  of  the 
quadrupole  moment  was  estimated  using  the  data  in  Fig.  6a  of  Ref.  2. 

Using  this  model,  an  isolated,  well-defined  local  maximum  in  the 
velocity  measured  on  the  high-speed  side,  associated  with  a corresponding 
well-defined  local  minimum  on  the  low-speed  side,  was  interpreted  as 
the  signature  of  the  passage  of  a single  vortical  structure  not  in  the 
process  of  pairing.  If  additional  significant  disturbances  were  discernible 
in  the  maximum -minimum  pair  (interpreted  as  resulting  from  the  presence 
of  additional  vortical  structures ),  the  pair  was  not  included  in  the  data. 
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By  way  of  example,  only  the  data  at  t = 0.65,  0.98,  1.33,  1.65,  2.35, 

2.  63  and  3.  35  seconds  in  figure  5 were  considered. 

The  maximum  induced  velocity  on  the  high-speed  side  and  the 
corresponding  minimum  on  the  low- speed  side  were  used  in  equations 
(3)  to  compute  the  position  >f  the  vortex  center  c and  total  circulation 
r of  each  vortical  structure.  The  resul’s  for  350  structures  are  pre- 
sented in  histogram  form  in  figures  6 and  7. 

b.  Cross-correlation  and  Auto- correlation  functicns 

The  correlation  functions  of  u ' = u„  - ut,  and  uT'  = uT  - u_  , 

H H H L»  Lj  L 

the  velocity  fluctuations  on  the  high-speed  and  low- speed  sides,  respectively, 

were  computed  from  the  filtered  data  by  means  of  the  equations 

M 

RHL(T)  = M ^ UH(tm)uL(tm  + T)  (4a) 

m=  1 

M 

“hH^  - T3  2 + T»  <4b> 

m=  1 
M 

RLL|T)  ‘ T3  2 UL(,m,UL(tm  + T>  <4c» 

m=  1 

where,  for  each  time  shift  T,  M is  the  largest  index  such  that 
t^  + T s T,  T being  the  record  length.  The  results,  averaged  over 
30  records,  are  plotted  in  figures  8 and  9.  The  smooth  curves  through 
the  data  points  have  been  drawn  to  aid  the  eye. 

4.  Results  and  Discussion 

a.  The  mean  core  position  (fig.  6)  is  slightly  below  (i.  e.  , on 
the  low-speed  side)  the  ray  c = 0 on  which  Um  = i(Ux  + TJ» ).  The 
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small  dispersion  in  the  histogram  indicates  that  structures  not  in  the 
process  of  pairing  convect  downstream  with  their  cores  closely  spaced 
about  a single  ray. 

The  mean  circulation  T (fig.  7),  normalized  by  the  product  A^AU, 
is  given  by 


UU 

(i) 


where  AU  = Ux  - Up  is  the  velocity  difference  across  the  shear  layer. 
Note  that  from  the  boundary  condition 

“■  AU  , (6) 

L 

and  equation  (5)  we  estimate  £,  the  spacing  between  structures,  to  be 

- 3.9  . (7) 

°uu 

This  result  is  consistent  with  earlier  independent  measurements  of  this 
ratio*’ It  should  be  mentioned  that,  strictly  speaking,  the  boundary 
condition  given  by  equation  (6)  should  be  written  instead  as  (T / 1)  = AU. 
The  use  of  the  ratio  of  the  mean  quantities  is  justified  a posteriori 
in  this  case,  since  the  distribution  of  T is  narrowly  peaked  about  its 
mean  value.  See  fig.  7. 

b.  Cross-correlation  and  Auto- cor  relation  results 


Results  of  other  investigators^'  ' suggest  that  the  vortical 
structures  in  a mixing  layer  are,  by  and  large,  cylinders  of  elliptical 
cross-section  whose  major  axes  may  be  inclined  with  respect  to  the 
streamwise  direction.  (See  fig.  3).  If  the  angle  of  inclination  is 
positive,  the  velocity  on  the  high-speed  side  will  reach  its  maximum 
before  the  low- speed  side  velocity  reaches  its  minimum.  This  would 
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cause  the  minimum  of  the  cross-correlation  function  of  the  high-speed 
and  low- speed  velocity  fluctuations  to  occur  at  a positive  time  shift. 

It  should  be  noted,  however,  that  the  minimum  of  the  cross-correlation 
function  has  occurred  at  a negative  time  shift.  See  fig.  9.  We  believe 
that  pairing  is  responsible  for  this  result. 

In  pairing,  neighboring  vortices  rotate  around  each  other  and 
coalesce  to  form  a larger  one.  Using  fig.  3 of  Ref.  5,  a schematic 
diagram  of  this  process  and  the  resulting  induced  velocity  fluctuations 
are  sketched  on  fig.  10.  It  can  be  seen  that,  in  the  process  of  pairing, 
a pair  of  vortices  taken  together  can  be  assigned  an  effective  negative 
angle  of  inclination. 

The  pairing  process  also  affects  the  auto-correlation  functions 
of  the  velocities  on  each  side  of  the  mixing  layer  (fig.  8).  The  auto- 
correlation function  of  the  high-speed  side  data  gives  a smaller  character- 
istic time,  T (twice  the  time  to  the  first  minimum),  than  that  for  the 
low-speed  side  (i.  e.  , T ).  This  becomes  clear  when  one  considers  the 

J-J 

fact  that,  in  the  process  of  pairing,  the  vortex  on  the  high-speed  side 
must  be  accelerating,  while  the  vortex  on  the  low- speed  side  is  decelerating. 
Note  that  this  is  consistent  with  the  results  of  Ref.  7.  It  should  also  be 
noted  that  the  data  support  the  relation 

thi.  “ i|TH  + V • |8) 

where  T T is  the  characteristic  time  determined  from  the  cross- 
correlation  function  (twice  the  time  interval  between  the  minimum  and 
maximum). 

If  we  use  the  characteristic  time  derived  from  the  cross- 

correlation function  as  a measure  of  the  mean  time  between  the  passage 
of  the  large  structures,  the  mean  spacing  would  then  be  given  by 
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1 THLUc  ' 


where,  from  the  cross-correlation  function,  we  have 

thlVx  “ °-  40  • 

Combining  equations  (1),  (9)  and  (10)  yields 


l/b 


uu 


4.  0 


(9) 


(10) 


(ID 


Note  that  the  independent  estimates  of  f/6^  as  given  by  (7)  and  (11)  are 
consistent. 

Direct  examination  of  the  velocity  records  shows  that  the  maximum 
velocity  on  the  high-speed  side  and  the  minimum  velocity  of  the  low- 
speed  side,  corresponding  to  the  passage  of  a single  structure,  occur 
at  the  same  time,  within  the  accuracy  of  the  present  measurements. 

Note  the  time  markers  on  fig.  5.  This  indicates  that  if  the  vorticity 
distribution  is  elliptical,  the  inclination  angle  is  small.  This  is  consistent 
with  the  findings  of  the  stability  analysis  by  Moore  and  Saffman  (Ref.  8). 

6.  Conclusion 

Simultaneous  records  of  the  streamwise  component  of  velocity  at 
the  high-speed  and  low- speed  edges  of  the  mixing  layer  were  used  to 
calculate  the  core  position  and  strengths  (circulation)  of  the  vortical 
structures.  It  was  found  that  an  isolated  structure  (not  in  the  process  of 
pairing)  travels  with  the  center  of  its  core  very  close  to  the  ~ y/x  on 
which  the  mean  velocity  is  Um  = i(Ux  + Ut ) and  has  a mean  circulation 
given  by  r/fi^AU  “■  3.9.  The  small  dispersion  about  the  mean  of  the 
distributions  of  core  position  and  circulation  is  noteworthy. 
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FIGURE  1.  VELOCITY  RECORD  TAKEN  AT  A MEAN  HIGH  SPEED 
FREE  STREAM  VELOCITY  OF  U.  = 56  CM/SEC 


as 


Sample  Size:  350 
Mean:  -0.025 


FIGURE  6.  DISTRIBUTION  OF  NORMALIZED  VORTEX  CORE  OFFSET 
DISTANCES 


FIGURE  7.  DISTRIBUTION  OF  NORMALIZED  VORTEX  STRENGTHS 
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The  turbulent  stuctures  in  the  flow  of  a jet  were  stuuied 
through  the  use  of  luser  induced  fluorescence.  The  fluip  to  be  issued 
from  the  jet  was  laoelled  with  a aye  that  fluoresced  resonantly  when 
illuminated  with  laser  light.  The  linearity  of  the  dye,  combined  with 
the  freedom  to  illuminate  the  flow  in  mar\y  ways,  makes  this  a powerful 
and  versatile  vecnnique.  dome  preliminary  observations  made  in  this 
fasnion  are  discussed. 
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The  study  of  turbulent  flows  has  been  curried  out  primarily 
by  qualitative  flow  visualization  techniques  ana  uy  quantitative  single 
point  measurements.  The  difficulties  encountered  in  correlating  the 
single  point  measurements  in  space  ana  time,  ana  interpreting  the  results 
if  the  correlations  are  carried  out,  huVe  limited  most  observations  to 
be  time  averaged  descriptions  of  the  flow,  buch  a description  neces- 
sarily contains  little  information  concerning  the  turbulent  structures 
tnemselves.  Consequently,  a detailed  description  of  the  nature  of 
turbulent  structures  has  not  yet  evolved  due  to  the  lack  of  adequate 
measurement  techniques. 

A detaileu  picture  of  the  structures  is  necessary  for  a workirg 

ki.owlege  of  the  piocesses  that  involve  turoulence.  buch  information 

would  contribute  to  the  understanding  of  energy  loss  due  to  turbulence 

and  especially  to  the  understanding  of  turbulent  dissipation  and  muting, 

as  for  chemical  reactions.  In  chemical  reactions,  such  as  combustion, 

the  reaction  occurs  at  the  interface  of  the  mixing  substances  and  this 

interface  is  defined  uy  the  shape  of  the  turuulent  structures. 

1,2,3  2 

Kotsovinos  and  List  have  constructed  a comprehensive 

description  of  the  properties  of  plane  (two  dimensional)  turbulent  jets, 
including  uouyancy  effects.  Their  description  includes  integral  prop- 
erties of  the  flow,  such  as  mass  flux  anu  momentum  flux,  and  fluctuation 
properties,  such  as  velocity  flueuations  anu  jet/ ambient  fluid  interface 
crossings. 

From  similarity  arguments,  which  are  supported  by  experiments 

> 

the  momentum  flux  through  aqy  plane  perpendicular  to  the  flow  is 
equal  to  a constant.  Thus:  (Fig.l) 
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where  J is  the  momentum  flux,  u is  tne  time  averagea  velocity  which 
has  its  maximum  on  the  jet  axis,  anu  b(x;  is  the  local  ^et  flow  wiath. 

U 

Dimensional  analysis  gives  (Lanaau  ana  Lifshita  ) 


(2)  b(x)~x 


so,  suostitution  into  (1)  gives: 
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This  is  not  a strictly  accurate  description  of  the  flow  in 

that  tiie  entrainment  flow  can  introuuce  a net  momentum  flux  that  invali- 

1 

dates  equation  (1),  as  Kotsovinos  noticed,  however  the  sign  of  the 
correction  due  to  the  entrained  flow  is  aependent  upon  the  jet's  envir- 
onment ana  since  the  correction  to  the  momentum  flux  is  ^17/b  at  an 
x/h  of  100,  tne  effect  can  generally  te  ignored.  The  change  in  ihe 
flow  due  to  the  nonconsuant  J is  that  the  growth  of  the  jet  width,  b(x), 
hecomes  nonlinear,  a decrease  in  the  momentum  flux  resulting  in  a 
great* r spreading  rate. 

5 

Uoertlor's  solutions  for  the  velocity  profile  of  the  turbulent 
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The  fluctuation  ana  intermittency  measurements  perform* -u  by 
Kotsovinos  inuicute  that  lor  the  two  dimensional  ,,et  no  ambient  fluid 
reaches  the  centerline  of  the  flow.  ihoy  inuicate  that  the  mixing  ana 
fluctuations  increase  on  each  sad**  of  the  centerline  before  going  to 


zero  again  in  the  ambient  fluid. 
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Profile  of  frequency  of  interface  crossing 
That  ambient  fluiu  does  roacn  the  jet  axis  wilL  l>c  born  out  uy  tills 
investigation,  previous  measurements  probably  being  limited  by  resol- 
ution considerations,  in  fuct,  part  of  this  study  was  directed  tow..rd 
the  structures  of  ambient  fluid  on  tno  axis  of  a three  dimensional  jet 
and  the  presence  of  ambient  fluid  on  the  two  dimensional  jet  axis  w.  s 
also  observed. 

What  is  missing  from  the  description  so  far  is  d.  ta  that  is 


correlated,  whether  over  space  at  a given  time  or  ovoi  time  .it  given 

<> 

spatial  cooruinates.  Perry  ana  him  have  developed  ..  technique  simll.,r 
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to  tiiot  used  in  this  investigation  that  enables  turn  to  make  such 
oosei-vations  of  the  entire  flow,  They  use  laser  illumination  of  a 
llow  la  celled  with  a smoke  tracer  to  visualize  a turmlent  flow  in  a 
gaseous  medium. 

The  technique  employed  in  this  investigation  also  uses  laser 
illumination,  but  here  the  jet  fluiu  is  laoelled  with  a fluorescent  dye. 
The  dye  is  Rhodamine  6-G,  which  is  used  in  (tye  lasers  and  fluoresces 
resonantly  when  illuminated  with  green  light.  The  fluorescence  is  at 
a longer  wavelength,  in  the  yellow,  and  is  excited  strongly  oy  the 
51kb>  & line  of  an  argon  ion  laser.  The  laser  liglt,  when  focused  to 
a tliin  pencil  or  sheet,  can  be  used  to  intersect  the  flow,  causing  the 
fluiu  that  was  issued  from  the  jet  to  fluoresce  while  it  simply  passes 
tiirougn  tne  amoient  fluiu.  The  fluorescence  is  unaffectea  by  the  flow 
(unlike  light  scattered  oy  smoke  particles  which  can  ue  rescatterea)  so 
that  observations  of  the  center  of  a flow  can  oe  maue.  nny  scattereu 
laser  light  is  easily  removed  by  an  appropriate  optical  filter. 

lip  to  near  saturation,  the  fluorescence  of  R6G  is  proportional 
to  the  incident  intensity  of  the  illuminating  beam,  bince  the  fluores- 
cence is  resonant  and  consequently  quite  strong,  it  is  not  necessary  to 
even  approach  the  point  of  saturation  in  practice,  so  the  fluorescent 
light  proviues  a quantitative  representation  of  the  jet  fluid  concen- 
tration in  the  flow. 

The  lasei  light,  as  descrioeu  by  gaussian  optics,  can  be 
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manipulated  to  proviue  illumination  that  is  sufficiently  precise  to 
resolve  Uie  finest  scale  of  turbulence  with  vitually  infinite  time  resol- 
ution. In  this  way,  the  technique  of  laser  induced  fluorescence  can  be 


useu  to  retrieve  nearly  ur^y  information  on  the  0et  fluiu  concentration 
fielu  as  a function  of  time,  ty  ooservmg  scattering  from  small  particles 
in  the  flow  for  a short  interval,  it  is  also  possible  to  characteiize 
the  velocity  field,  yielding  a substantial  amount  of  information  on  the 

flow. 


APfAhATUS 


both  two  and  three  dimensional  jets  were  studied  during  this 
investigation,  although  the  emphasis  was  on  the  flow  from  a 0.3" 
diameter  three  dimensional  jet.  This  jet  with  its  reservoir  of  aycd 
water  was  suspended  above  a glass  aquarium  (2'x  1'  x 1*-'  deep)  con- 
taining unqyed  water.  The  flow  from  the  jet  was  regulated  by  metering 
the  introduction  of  air  into  the  „et  reservoir.  A wide  range  of  jet 
exit  velocities  was  realized  by  tliis  arrangement,  making  Reynolds 
numbers  of  zero  to  several  thousand  accessible.  This  system  was  con- 
ceived and  designed  by  dimotakis.  a two  dimensional  Jet  with  a variable 
slit  width  (1-6  mm  x 23.$  cm)  designed  by  the  author  was  also  studied  in 
the  course  of  the  investigation. 

An  argon  ion  laser,  typically  outputting  several  hundred 
milliwatts,  was  the  illumination  source.  In  the  various  methods  of 
illumination  u^ed  here,  this  was  enough  light  intensity  when  used  in 
conjunction  with  a 0.6  mg/liter  concentration  of  i<6d.  This  concentration 
was  sufficiently  small  that  the  fluorescent  losses  oiu  not  attenuate  the 
illuminating  ueam  significantly  in  its  traverse  of  the  flow.  In  more 
quantitative  experiments,  it  m;iy  oe  necessary  to  consider  the  attenuation 


(6) 

out  a calibration  of  the  fluorescence  anil  the  calculation  of  the  losses 
is  not  difficult  to  carry  out. 

To  ooserve  a line  intersecting  the  flow,  a long  focal  length 
lens  focused  the  beam  to  a waist,  as  described  oy  gaussian  optics.  The 
waist  diameter  ana  the  length  over  which  the  beam  hau  not  appreciably 
expanded  permitted  viewing  of  the  smallest  structures  in  the  flows 
under  consideration  here.  To  observe  a plane  in  the  flow,  diverging 
cylindrical  lenses  wore  placed  in  the  ceam,  expanding  the  oeaiu  in  one 
dimension  perpendicular  uo  its  axis  of  propagation.  The  beam  was 
manipulated  with  a mirror  or  prisms  to  direct  it  through  the  flow  in 
the  desired  direction. 

The  resulting  fluorescence  was  recorded  in  two  different  ways 
for  these  two  observational  schemes,  une  method  consisted  of  illuminating 
the  flow  with  the  sheet  of  light  and  pliotographing  the  fluorescence 
with  a conventional  camera  focusta  on  the  plane  of  the  sheet.  The  second 
methou  used  a pencil  of  light  to  illuminate  the  flow  and  the  concentration 
pattern  along  this  line  was  ooserved  as  a function  of  lime.  The  data 
from  the  first  method  can  be  displayed  as  concentration  as  a function  of 
x arid  y at  a fixed  time.  The  second  metnous  results  can  be  displayed 
as  concentration  as  a function  of  x (or  y)  ana  t for  a fixed  y (or  x). 

To  collect  the  concentration  data  along  the  line  as  a function 
of  time,  a linear  optic  array  camera  was  used.  (fig.  1J  This  camera, 
designed  by  the*  author  with  bimotakis  for  tnis  use,  consists  of  a 
neticon  linear  array  of  photodioaes  and  its  associated  electronics  and 
an  optical  imaging  system,  a Nikon  photographic  lens  images  the  illum- 
inated line  on  the  array  of  oioues. 
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Although  there  exist  ooth  100  and  1000  element  arrays  that 
can  be  used  in  the  camera,  only  the  1000  element  array  (each  dioue  2J>  um 
square,  total  length  1")  was  useu  i'or  these  measurements.  With  an 
external  clock  and  power  supply,  tile  camera  outputs  sequentially  a 
signal  from  each  element  that  is  proportional  to  the  integrated  intensity 
that  element  received  since  the  last  sampling.  The  camera  could  be 
clocked  more  quickly  ( 200  samplings/sec)  than  the  fastest  fluctuations 
observed  in  the  flows,  permitting  observations  that  were  not  time  - 
resolution  limited. 

The  output  of  the  camera,  then,  is  a concentration  versus 
position  representation  of  the  flow  when  the  signal  is  viewed  as  a 
function  of  time  on  an  oscilloscope  for  a single  sampling  of  the  array. 
The  changes  in  tnis  functional  relationship  give  the  desired  time 
correlated  information.  To  get  a concentration  versus  position  and 
time  representation,  liorrison  uesigneo  an  optic  array  scanner  that 
synchronously  sweeps  an  oscilloscope  screen  and  mouulates  the  oeam 
intensity  with  the  output,  of  the  camera.  The  resulting  oscilloscope 
tube  image,  much  like  that  of  a radar  screen,  was  recorded  photograph- 
ically, giving  the  desired  concentration  versus  x (or  y)  and  t repre- 
sentation. 
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Industriol  water  passed  through  an  inline  10  um  filter  was 
the  fluid  usee  throughout  tliis  investigation.  Even  though  laser  light 
scattered  by  particles  could  be  optically  filtered,  it  was  deemed 
important  to  keep  the  water  as  clean  as  possible,  a Hoy a HMC  0(G) 
filter  was  used  to  preferentially  pass  the  fluorescent  light  while 
blocking  the  scattered  light. 

There  was  one  arrangement  where  the  particles  remaining  in 
the  water  were  useu  to  advantage.  (Fig.  2a.)  Without  any  dye  in  the 
jet  reservoir  and  with  a sheet  of  light  (xy  plane)  along  the  jet  axis, 
a velocity  field  representation  of  the  flow  was  recorded.  This  was 
accomplished  by  taking  a short  time  lapse  photograph  of  the  light 
scattered  by  the  uniformly  distributed  particles.  The  length  and 
direction  of  the  streaks  recorded  on  the  film  amount  to  a two  component 
(xy)  vector  velocity  field  for  the  flow. 

As  with  all  the  fluorescence  measurements,  these  velocity 

field  measurements  were  recorded  over  a range  of  Reynolds  numbers  where 

the  flow  was  turbulent  (500-1*000).  similarity  of  the  flows  was  preserved 

over  this  range  to  the  extent  that  it  was  investigated  in  this  study. 

That  is  to  say  the  growth  of  the  jet  width  and  growth  of  the  turbulent 

structures  was  similar  over  this  range.  The  importance  of  keeping  the 

densities  ( a. e.  temperatures)  of  the  jet  reservoir  iluid  and  the  ambient 

fluid  the  same  was  observed,  the  reason  being  that  the  effects  of  bouy- 

1 

anc.»  can  alter  the  flow  considerably  (Kotso vinos  ). 

Perhaps  the  most  intuitively  suggestive  measurement  is  the  sheet 
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illumination  of  the  flow,  'i'he  flow  was  illuminated  with  a sheet  of 
light  directed  upward  through  the  bottom  of  the  tank  (aquarium)  as  in 
the  velocity  measurement.  (Fig.  3a.)  However  the  scattered  light 
was  blocked  in  this  case  and  the  vet  fluid  was  dyea,  yeiloing  a con- 
centration field  representation  of  the  flow,  a conventional  photo- 
graphic camera  using  U00  AdA  Tri  X film  recorded  tnis  representation  at 
a fixed  time  value  sometime  after  the  flow  had  established  itself. 

Tnis  means  of  viewing  the  flow  was  intended  to  show  turbulent  structures 
of  all  scales,  their  variations  with  x/u  and  x/y,  and  the  spatial 
coherence  ana  ordering  of  the  structures,  (rig.  2b.) 

The  time  coherence  of  the  structures  is  also  of  interest, 
as  are  the  propagation  velocities  of  the  structures  (in  a structural 
sense  as  opposed  to  the  vector  velocity  of  points  vjithin  the  structures, 
as  was  previously  discussed).  Two  different  methods  were  employed  to 
study  time  variations  in  the  flow.  The  first  consisted  of  viewing  a 
line  intersecting  the  axis  of  the  jet  perpendicularly.  This  is  a 
concentration  versus  y anu  t representation  of  the  flow  at  a given  x 
value,  (rigs.  3 b.  and  U)  The  second  scheme  was  to  view  the  axis  of 
the  jet  as  a function  of  time;  an  xt  at  fixed  y (y*0)  representation. 

If  uhe  commonly  held  idea  that  no  ambient  fluid  reaches  the  jet  axis 
were  true,  this  would  have  been  a singulatl.,  uninteresting  measurement. 
The  opposite  turned  out  to  ue  the  case.  (Figs.  3 <-•.  ana  J>) 

oince  the  illuminating  ocam  was  made  as  small  as  possible  u> 
Observe  the  finest  structures  and  the  linear  array  of  photodiodes  is 
only  25  um  (1  mil)  wide,  the  alignment  of  the  system  was  not  a simple 
task,  (oee  APFhNhlX)  The  width  of  the  sheet  or  the  diameter  of  the 
pencil  can  be  calculated  from* 
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where  ltu  is  the  radius  oi'  the  laser  beam,  f is  converging  lens  focal 
length,  ana  r is  the  ui stance  along  the  oeam  axis  from  the  focal  point, 
being  x or  y in  the  two  cases,  do  for  51U5  X,  f a 1 m,  and  Rq  « 2ram, 
the  half  width  or  radius  of  the  illumination  is  0.082  mm. 

In  the  yt  method,  the  beam  was  brought  in  through  the  end  of 
the  tank,  intersecting  the  vertical  jet  axis  horizontally.  In  the  xt 
method,  the  beam  was  directed  up  through  the  bottom  of  the  tank.  In  all 
of  the  measurements,  the  flow  close  to  the  bottom  was  not  studied  because 
of  the  interruption  of  the  flow  due  to  the  bottom  and  its  effects  further 


upstream* 


dome  preliminary  observations  of  a two  dimensional  jet  designed 


anu  ouilt  for  this  investigation  were  also  carriea  out.  These  observa- 
tions consisted  of  the  instantaneous  concentration  fiela  in  a plane. 
When  the  plane  was  aligned  along  the  direction  of  the  flow  and  perpen- 
dicular to  the  slit,  a concentration  pattern  much  like  that  of  the  xy 
photographs  of  the  three  dimensional  jet  was  recorded.  When  the  plane 
was  aligned  again  along  the  direction  of  the  flow,  but  along  the  slit 
insteau,  a means  of  determining  the  two  dimensionality  of  the  flow  was 
realized.  ( Kigs.  6 and  7) 
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A comparison  of  the  concentration  field  ana  the  velocity 
lield  (figs.  2 a.  and  b. ) allow  a number  of  ooservations.  Where  the 
flow  is  still  laminar  near  the  jet  exit,  the  velocity  field  inuicau  s 
that  there  is  consicerable  entrainment  of  ambient  fluid  into  the  flow. 
(Xue  bright  spot  in  the  photograph  is  a reflection  of  the  jet  exit 
from  the  back  of  the  tank.)  apparently  the  simple  presence  of  entrain- 
ment aoes  not  necessitate  turbulence,  although  such  a laminar  flow  with 
so  considerable  an  entrainment  may  be  a metastable  state  that  eventually 

does  become  turbulent. 

7 

Townsena  aescribes  a model  for  entrainment  based  on  the 
diffusion  of  vorticity  into  the  ambient  fluid  through  the  flow's  bound- 
ary. The  ambient  fluid  is  considered  to  be  irrotational  outside  the 
boundaries  while  the  fluid  just  inside  the  boundary  is  gaining  vorticity 
by  a uniform  (on  the  scale  of  the  smallest  structures  present)  diffusion 
process.  Aqy  large  scale  disturbances  in  this  pattern,  that  might  trap 
irrotational  fluid  in  the  flow  for  instance,  are  expected  to  be  due  to 
normal  statistical  fluctuations. 

The  entrainment  and  the  acquisition  of  vorticity  are  not 
dil'fu sion-like  but  rather  have  some  definite  large  scale  characteristics 
in  this  flow  as  can  be  seen  in  rig.  2 a.  Large  scale  vorticity  extends 
deeply  into  the  ambient  fluid,  setting  up  a flow  that  must  have  some 
deterministic  nature.  Clearly  a moael  of  uniform  entrainment  does  not 
apply  in  the  ambient  fluid  just  outside  the  jet. 

The  instantaneous  concentration  fields  for  the  three  dimen- 
sional jet  (rig.  2b.)  and  the  two  dimensional  jet  (rig.  6)  are  quali- 
tatively much  the  same,  after  an  initial  laminar  legion,  the  flow 
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becomes  unstable  ;.nu  turbulent.  The  initial  instability  has  a form  that 
is  similar  in  these  two  cases  althogh  it  is  simpler  in  the  two  dimensional 
jet.  Because  of  the  relative  simplicity,  which  is  due  to  the  one  less 
degree  of  fre«.uom  in  two  dimensions,  the  x'ollowing  descriptions  are  of 
the  plane  ^et. 

There  are  two  distinctly  different  types  of  transition  to 
turbulence.  The  abrupt  change  between  the  types  occurs  at  a heynolos 
number  (^100)  where  the  flow  at  the  exit  shows  a stable  periodic 
vorticity.  The  transition  to  turbulence  for  a slow  jet  seems  to  come 
about  after  a developing  periodic  wander  attains  a critical  rotation, 
at  which  point  a vortex  is  created,  (tig.  10  a.)  In  the  case  where 
the  flow  already  has  a periodic  vortex  structure}  the  faster  flow;  the 
large  vortices  appear  where  the  wandering  begins,  (fig.  10  b, ) 

Apparently  the  flow  already  lias  sufficient  vorticity  to  be  turbulent 
ana  ary  asymmetry  in  the  flow  is  enough  to  trigger  the  large  vortices. 

Once  the  large  scale  vortices  appear  in  the  flow,  all  the 
chaotic  connotations  of  turbulence  are  realised,  a general  periodic 
wandering  trend  is  still  evident  in  the  flow  in  many  records,  but  its 
coherence  is  not  especially  great  - the  trend  sometimes  disappears. 


Compare,  for  example,  tig.  2b.,  which  is  fairly  typical,  to  tig.  9 in 
wliich  the  meandering  is  quite  pronounced.  But  beyond  any  shaaow  of  uoubt 


ambient  fluid  does  reach  the  jet  axis. 


It  should  be  noted  here  that  when  a wijietf^ sheet  of  light  was 
used  to  illuminate  the  flow,  the  ambient  fluid  on  the  0et  axis  was  not 


\ 


evident  at  small  x/Jj.  This  implies  that  the  structures  on  tne  axis  are 
fine  scale,  especially  closer  to  the  exit,  and  previous  measurements  may 
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have  not  found  them  due  to  a lack  of  resolution.  Kotsovinos  reports 
that  the  region  immune  to  ambient  contamination  narrows  with  x/D.  This 
is  consistent  with  the  similarity  result  that  the  scale  of  the  structures 
increases  with  x/d.  bon-equently  the  ability  to  resolve  the  structures 
improves  until  the  finest  structurs  are  observed  and  ambient  fluiu  is 
noticed  on  the  axis.  The  difficulties  in  observing  the  fine  scale  of 
the  structures  is  aggravated  uy  the  fact  that  the  velocity  is  at  a 
maximum  on  the  jet  axis  anu  so  time  resolution  is  important  as  well. 

In  checking  the  two  dimensionality  of  the  plane  jet,  striations 

in  the  flow  were  noticed  along  the  length  of  the  slit.  (Rig.  8)  Similar 

8 

phenomena  have  been  noticed  oel'ore,  in  particular  by  Konrad  in  a shear 
layer.  He  begins  to  explain  their  occurence  in  terms  of  a kayleigh 
instability,  which  comes  about  when  a structure  is  rotating  in  a sense 
opposite  to  that  of  the  more  general  flow.  This  model  turns  out  to  lit 
in  well  with  the  Keynolds  number  dependence  of  the  phenomenon,  a depend- 
ence that  appeared  to  be  present  to  some  degree  in  the  jet  as  well,  out 
by  polishing  the  contraction  surfaces  ol'  the  jet  and  taking  precautions 
to  avoid  their  contamination  with  particles  or  bubbles,  the  striations 
were  made  to  disappear.  This  effect  might  indicate  that  the  instability 
is  still  triggereu  by  imperfections  in  the  flow. 

Even  without  the  striations,  the  two  dimensionality  of  the  flow 
is  limited.  (Rig.  7)  The  initial  region  shows  good  two  dimensional 
behavior,  but  with  the  onset  of  large  scale  vorticity,  three  dimen.  ional 
structure  develops.  There  is  still  some  overall  two  dimensional  Leila vior 
seen  in  the  bands  of  turbulence,  but  this  has  the  same  sort  of  limited 
coherence  that  the  wandering  was  seen  to  have  had. 


The  incre-ose  of  the  turbulent  scale  with  x/h  as  given  by 
similarity  is  also  evident  in  rigs.  2 b.  and  6.  The  spreading  ol'  the 
jet,  aue  to  entrainment,  is  very  nearly  linear  as  previously  noted 
(b(x)fox).  Although  the  turoulence  is  not  uniform  within  its  boundaries 
it  doc  s have  linearly  spreading  edges.  That  the  effects  of  the  large 
vortices  extends  into  the  ambient  fluid  i.,  consistent  with  the  fact  that 
tile  edges  of  the  jet  are  not  maoe  up  of  uniform  structures.  The  concept 
of  a jet  boundary  is  still  valid,  but  the  rotational  influences  of  the 
flow  are  nut  reutricteu  to  oe  within  those  limits. 

it  nay  be  somewhat  sui’prising  that  the  size  of  structui'es 
uepends  oily  on  x/D,  while  there  is  considei’able  variation  in  the 
velocity  aci'oss  the  flow,  as  a function  of  the  y coordinate.  This 
result  implies  that  structui’es  in  different  environments,  moving  at 
different  velocities  change  their  scale  in  the  same  way.  Nevertheless, 
tnis  is  predicted  by  similarity  and  it  is  born  out  by  this  investigation 
Observing  yt  records  shows  that  structures  at  the  edges  of  the  flow  are 
elongated  in  the  t direction,  (rig.  h)  Thus  they  are  moving  more 
slowly,  yet  the  xy  records  show  that  in  spatial  coordinates  the  scale 
is  preserved  across  the  flow. 

The  remaining  thechnique  used  to  analyse  the  flow  was  the  use 
of  axial  illumination.  The  signal  was  inverted  for  ttii  s technique  so 
that  the  amoient  fluid  appears  as  oright  regions.  This  was  convenient 
otcause,  although  some  amoient  fluiu  cues  reach  the  axis,  the  greatest 
part  of  the  flow  w«s  issued  from  the  jut  and  a clearer  record  is  maoe 
with  fewer  oright  regions. 

The  average  velocity  along  thf  tliree  dimensional  jet's  axis 
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it 


goes  us  the  inverse  of  the  distance  along  that  axis.  On  an  xt  record, 
the  velocity  of  a structure  is  the  slope  of  its  path.  Thus  we  have: 

dx  1 2 

(6)  at  x or  x ~ t 

uo  the  paths  shoula  oe  parabolic,  in  an  average  sense,  That  this  is 
ruughly  correct  can  oe  seen  in  rig.  5*,  but  a larger  test  section  would 
be  necessary  to  verify  tins  observation. 

This  measurement  turned  out  to  be  particularly  revealing  in 
that  it  permitteu  the  observation  of  how  the  flow  grows  due  to  the 
interaction  of  the  turbulent  structures.  The  growth  of  the  turbulent 
scale  raises  the  question  of  the  mechanism  of  this  growth,  leather  than 
each  vortex  expanding  (which  might  wreak  havoc  with  the  velocity  depend- 
ence J the  flow  causes  vortices  to  combine,  becoming  part  of  larger 
structures.  This  can  be  seen  in  the  anti-branching  patterns  that  lie 
between  converging  parabolas.  The  complicated  structure  of  the  vortices 
sometimes  give  the  impression  of  anomolous  velocities;  a round  structure 
crossing  a lint  obliquely  can  have  botli  positive  ana  negative  apparent 
velocities.  (rig  11)  an  example  of  tliis  effect  is  recoroed  in  ilg.  £ 
with  the  nestled  arches  of  a spiraled  vortex  lx  clocks  from  the  left  ano 
0.75  blocks  from  the  bottom  of  the  figure. 

dome  insight  into  the  nature  of  the  turbulent  structures  is 
provideo  by  the  the  observation  that  their  trajectories  are  basically 
straight  lines  on  the  xt  record.  The  major  exception  to  this  rule  is 
when  several  structures  are  close  together,  in  ttie  process  of  combining. 
The  parabolas  are  made  up  of  line  segments  which  are  parabolic  only  in 
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an  average  sense.  This  provides  an  analog  to  a free  vortex  ring  which 
travels  with  a constant  velocity  through  the  ambient  fluid.  It  might 
oe  possiuie,  then,  to  formulate  an  atomistic  model  of  turbulence  in  which 
the  structures  behave  like  independent  vortex  rings  that  only  interact 
when  they  are  close  together.  This  woulo  not  be  a complete  description 
though,  because  it  coulu  not  explain  the  ptriouic  wandering  that  is 
observed  in  the  entire  flow. 


CUWCLUJ1UNS 


a mt-tnou  for  studying  the  structures  in  a turbulent  flow  has 
been  developed.  It  has  been  employed  to  make  detailed  qualitative 
measurements  of  the  flow  of  a turbulent  yet  and  has  the  potential  to 
permit  quantitative,  as  well  as  qualitative,  studies  of  many  flows. 

In  two  and  three  dimensional  jets,  ambient  fluid  was  observed 
to  reach  the  ^et  axis.  The  nature  of  the  onset  of  turbulence  was  studied 
and  it  appears  to  be  due  to  an  instability  caused  by  the  wandering  of  the 
flow.  The  growth  of  the  jet  was  seen  to  ue  accomplished  b„«  the  combina- 
tion of  turbulent  structures  as  the  0et  width  increases.  These  structures 
appear  to  travel  with  constant  velocity  like  free  vortex  rings  except 
while  in  the  process  of  combining  to  form  larger  structures. 
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Fig*  1 Experimental  set  up  and  definition  of  axes 
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hPPENuU 


The  linear  optical  system  is  an  important  part  ox’  the  laser 
induced  fluorescence  experiment.  The  illuminating  beam  from  tne  laser 
is  directed  through  the  fluid  flow,  exciting  ax^y  of  the  fluorescent  ctye 
in  its  path.  The  fluorescent  light  is  collected  by  a lens  and  focused 
upon  an  array  of  photodiodes  for  its  measurement.  (Fig.  Al)  The  illum- 
inating beam  is  a narrow  pencil  of  light  and  the  photooiooe  array  is  a 
line  of  small  elements.  This  could  conceivably  cause  some  difficulty  in 
trying  to  bring  the  image  of  the  segment  of  the  illuminated  flow  to  be 
observed  into  coincidence  with  the  diode  array.  Due  to  several  properties 
of  the  optical  system,  it  is  in  fact  difficult  to  determine  when  the 
optical  system  is  optimally  aligned. 

The  beam  from  the  laser  is  directed  through  a lens  that  x'ocuses 
the  beam  ai  the  center  of  the  line  segment  in  the  fluid  that  is  to  be 
observed.  The  converging/di verging  do am  can  be  described  accurately 
by  Gaussian  optics  - so  the  intensity  is  gaussian  in  the  distance  from 
the  axis  of  the  oeam  and  the  intensity  distribution  becomes  most  peaked 
where  the  beam  is  focused.  If  tnere  is  a constant  concentration  of 
fluorescent  dye  in  the  fluid  through  which  the  beam  is  passing,  the 
fluorescence  at  an,/  point  will  ot-  proportional  to  the  light  intonsity 
at  that  point. 

Gaussian  optics  gives  the  intensity  of  the  beam  as  a function 
of  position  as; 


a(1) 


l(x,y, z)  a Iq 


(30) 


where,  *(2) 


ana  wQ  is  the  waist,  the  smallest  radius  of  the  oeam,  given  by  Eqn.  (i>). 

The  fluorescing  dye  removes  some  intensity  from  the  ceam, 
resulting  in  an  exponentially  decreasing  total  intensity,  with  very  low 
aye  concentrations  these  losses  can  be  made  small  and  in  this  case  th*- 
amount  of  light  scattered  by  fluorescence  does  not  appreciai  ly  -le  : 
intensity  further  along  the  beam.  Thus  the  exponential  aecn  «.  • 
ignored  ana  the  fluorescence  becomes  proportional  to  l(x,y,i  . 

The  low  concentration  of  uye  provides  a nc«ns  o:  »*i«  i 
the  intensity  of  the  illuminating  beam  as  a functj.  : . 
fluorescence  can  be  measured  using  a collecting  1* ’ 
array  and  the  result  can  o«  compared  with  the 
see  if  the  array  it  aligned  to  be  coincio- id  . 
illuminated  segment,  oince  the  inunsit 
to  a high  degree  of  accuracy  uno  i>  ■ 
it  should  be  possible  to  c^lcul  i-  • • 
to  reach  the  array. 

To  a good  opp i • 
to  be  perfect.  That  is, 
point  in  the  f.  i 


alia  th*-  n 


uni 
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only  introduces  a constant  factor  proportional  to  the  aperature  size 
and  a cosine  factor  due  to  oblique  incidence  that  is  very  nearly  constant 
for  these  purposes  (within 

bue  to  the  changing  cross  section  of  the  bear!,  even  when  the 
axis  of  the  array  coinciues  (in  tiiree  dimensions)  with  the  image  of  the 
axis  of  the  illuminating  oeam,  the  light  intensity  that  falls  on  the 
photoaiode  array  will  have  a maximum.  This  maximum  occurs  where  the 
cross  section  is  most  peaked,  that  is  at  the  image  of  the  origin  where 
where  the  intensity  has  its  maximum  value.  If  the  dioue  array  is 
slightly  skewed  or  otherwise  misaligned,  it  is  still  likely  that  a peaked 
intensity  profit  will  be  observed,  but  it  will  have  a different  shape. 

In  practice  it  is  not  easy  to  tell  when  the  profilt corresponds  to 
optimal  alignment. 

To  calculate  how  the  collecting  optics  forms  an  image  in  the 
plane  of  the  array,  an  impulse  response  function  will  be  calculated. 

This  relates  the  amount  of  light  received  at  a point  in  the  array,  or 
image,  plane  to  that,  emitted  at  a point  in  the  illumination  field.  This 
function  can  then  oe  multiplied  by  the  intensity  actually  present  ana 
integrated  over  the  volume  where  there  is  a significant  intensity.  The 
impulse  response  function  will  depend  on  the  poj  nt  in  the  illumination 
field,  (x,y,z),  and  on  the  position  along  the  array,  y •.  The  integration 
of  the  product  of  the  impulse  response,  U(x,y,z;y«)  and  the  intensity, 

,)  will  yield  a value  for  the  signal,  S(y')  that  the  photodiodes 
will  measure. 

A point  source  emitting  equally  in  all  directions,  as  ooes  a 
point  in  the  dycu  fluid,  illuminates  an  apexature  of  radius  1 at  a 


1 


4 


2 2 

distance  uniform!  j if  f /^z0  & !•  'fhe  apcruture  of  the  lens  defines 


a cone  of  light  enitteu  fioin  u joint  that  is  then  focused  to  a point  in 
the  image  field,  (.-it'.  The  amount  of  light  from  this  point  that  is 

receiveu  at  a giv<  n detector  is  proportional  to  the  fraction  of  light  in 
the  converging  cone  that  is  inciui  nt  upon  the  detector.  The  fraction  of 
*-i{jht  collected  by  the  lens  that  is  incident  on  the  detector  is  given  by 
the  area  of  cone's  moss  section  that  ove/laps  the  detector  divided  uy 
the  total  area  of  the  cross  section,  (rig.  a.}) 

The  rauius  of  the  cross  section,  hc,  is  a function  of  the  z 
coordinate  of  the  emitting  point,  oince  any  point  in  the  xy  plane  is 
focused  u>  a point  in  the  plane  uei'ineu  by  z zQ  d0  the  lens  law  gives: 
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If  we  take  a point  not  in  the  xy  plane. 
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rVom  rig.  nU,  we  can  write. 
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If  hc  is  much  sinallei  than  the  size  of  a uetector,  the  cone 
of  light  can  be  consiuereu  focusea  and  the  response  of  the  system  will 
nave  a sharp  step  mien  the  focusea  point  moves  off  the  edge  of  tin 
detector,  fhe  Qetectors  are  2$  uw  square,  so  the  condition  for  a step  is. 


hc  fi  r 


<«  25  um 


riere  we  will  define,  as  is  conventional,  iV  3 f/2P  and  ii-z0/do,  wnere 
M is  the  magnification.  So,  for  ri»3  and  fi/~l.l*. 


z « 81*0  um 


The  question  now  arises  if  it  is  necessary  to  consider  z~81*0  um 
to  include  all  regions  where  there  is  appreciable  intensity  in  the 
integral  for  d(y').  we  will  have  to  consider  the  largest  values  of  3 
when  the  cross  section  is  the  least  peaked,  that  is  where  w(y)  is  large. 
With  a magnification  of  3,  the  detectors  will  view  a segment  7.62  cm 
long.  Thus  y can  take  on  values  between  -3.81  and  3.81  cm.  booking  at 
the  expression  for  w(y),  we  see, 


*(9)  wmaxa  W(ynex^  * 2^um  using  A*51U.£ 


nm,  wQc  2j>  um 


The  intensity  at  the  point  where  the  stop  condition  no  longei  holds 

2 2 

contains  the  term  exp(-2z  /w  ),  where  z 81*0  um  and  w 21*1*  um.  This 
term,  then,  introduces  a factor  of  !*■  x 10""^  indicating  that  there  is  no 
light  of  ar\y  significance  emitteu  where  it  is  not  valia  to  use  a response 
that  has  a sharp  step,  we  have  defined, 


mm*»y*** 


co 


*(10)  J(y') 


using  the  step  condition. 
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u(*>y>*;y')  i(x,y,z)  di 


-cr> 


h(11)  a(y 


'C 

-<*>  0 


*a/2  exp(-2(x^  z^j/w2) 


a/2  a/2 


w(yi‘ 


ox  dy  dz 


where  y0 ; rly  ' and  a-=i'l  2t  um  = 7i>  urn.  These  parameters  make  the  integration 
over  the  anti-image  of  tne  detector  that  lies  in  the  object  field.  We 
cun  integrate  z from  -co  to  co  since  there  is  no  appreciable  intensity 
outside  of  \z\<3U0  um  ana  its  inclusion  will  not  affect  the  results. 

The  function  w(y ) varies  slowly  as  y takes  on  values  from 
-3*tfl  to  3.01  cm,  separatee  by  76.2  um  as  y*  moves  from  detector  to 
detector,  do  w(y)  can  be  considered  constant  on  the  interval  (y0*a/2, 
y0-a/2)  and  equal  to  w(y0).  The  integral  over  2 introduces  a constant, 
wnich  can  be  incorporated  into  the  overall  constant,  and  cancels  a term 
w(y;  in  the  denominator.  The  intt gration  over  y alto  only  introouces  a 
constant.  We  now  have. 


pa/2  „ a£2 

A(12)  J(y')s  1,\  exp(-2x2/w2/  dx^I^S  exp(-2u £)  du 

i)-V2  w ^ 
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2.121 


To  account  for  the  index  of  refraction  of  the  watt  r through 
which  the  beam  passes,  it  is  convenient  to  define  y»»  y'/n*  y'/l.J3. 
the  measured  intensity  profile  is  slightly  narrower  than  this  calculated 
profile  (rig.  .1^)  due  to  some  misalignment. 
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Relative  intensity  vs.  array  position 
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'^The  more  important  results  and  developments  achieved  as  a direct  result 
of  this  contract  are  as  follows: 

A model  for  the  calculation  of  turbulent  flows  was  constructed  which 
incorporates  effects  of  Reynolds  stress  relaxation  with  the  generation  or 
destruction  of  turbulent  energy  by  rotation.  The  model  was  verified  by  com- 
parison with  experiment.  It  was  then  used  to  predict  the  properties  of 
isolated  turbulent  vortices  and  applied  to  the  turbulent  mixing  layer,  in 
particular  to  the  formation  and  interaction  of  the  oraani7Pd  structure  in  thf 
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mixing  layer.  An  alternative  approach  to  the  mixing  layer  in  which  it  is  model- 
led by  a rolling-up  vortex  sheet  between  fluids  of  different  density  has  also 
been  formulated. 

Calculations  of  strained  flame  elements  for  use  in  the  coherent  flame  model 
of  H2,  F combustion  showed  that  only  at  very  high  strain  rates  was  sufficient 
vibrational  nonequilibrium  produced  to  be  of  interest  for  a chemical  laser. 

This  result  is  reflected  In  the  fact  that  only  in  the  early  portion  of  the  jet, 
where  strain  rates  are  a maximum,  does  this  vibrational  nonequil ibrium  occur. 

Experimental  results  were  obtained  in  a new  facility  and  by  a new  technique 
for  the  amount  of  reaction  product  obtained  between  two  chemically  reacting 
aqueous  streams.  Results  at  high  and  low  Reynolds  number  showed  the  dependence 
of  the  mixing  on  Reynolds  number  and  Schmidt  number.  A transition  Reynolds 
number  region  was  found  and  the  nature  of  the  generation  of  small  scale  motions 
which  dominate  the  formation  of  reaction  product  was  explored. 

A new  facility  to  explore  energetic  reactions  with  substantial  heat  release 
has  been  designed,  is  being  built,  and  should  be  used  for  the  first  measurements 
in  the  summer  of  1979.  The  design  problems  were  formidable  but  the  facility  is 
unique  and  offers  the  prospect  of  obtaining  new  and  very  important  data  with 
which  to  compare  existing  models.  As  a result  of  the  design,  the  range  of 
problems  and  parameters  that  can  be  studied  in  the  facility  and  the  possibilitie 
for  the  development  and  application  of  new  instrumentation  and  diagnostics  are 
very  great. 

State-of-the-art,  very  high  speed  and  precise  laser  Doppler  velocity  meas- 
urements were  completed  under  this  contract.  The  acquired  experties  will  find 
application  in  the  turbulent  combustion  experiments.  It  has  already  been  ap- 
plied with  considerable  success  up  to  a Mach  number  of  2.2.  The  extension  of 
two-point  velocity  measurements  has  now  been  made  and  a multipoint  multichannel 
LDV  system  has  been  designed  and  is  nearing  completion. 

A laser  induced  fluorescence  technique  was  developed  and  successfully 
applied  to  the  observation  of  turbulent  mixing  in  a water  jet.  It  allowed  the 
direct  observation  of  a concentration  field  to  the  smallest  turbulent  scales. 
Both  direct  photographic  techniques  and  imaging  of  a line  on  a 1024  detector 
Recti  con  array  have  been  used  to  obtain  the  data.  Image  processing  will  give 
quantitative  measurements  of  the  concentration  field  and  potentially  an  inferred 
velocity  field. 
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